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Abstract 

I continue analysis of the Schrodinger operator with the strong degenerating mag- 
netic field, started in [Ivr6]. Now I consider 4-dimensional case, assuming that mag- 
netic field is generic degenerated and under certain conditions I derive spectral asymp- 
totics with the principal part x h~ 4 and the remainder estimate 0(fi~ 1 / 2 h~ 3 ) where 
H 1 is the intensity of the field and h <C 1 is the Plank constant; [ih < 1. 

These asymptotics can contain correction terms of magnitude // 5//4 /7~ 3 / 2 corre- 
sponding to the short periodic trajectories. 

Introduction 

0.0 Preface 

1 continue analysis of he Schrodinger operator with the strong degenerating magnetic field, 
started in [Ivr6] 

(0.1) A = \ QT Pjg Jk (x)P k - V) , Pj = Dj - fiVj 

j.k 

where g-* k , Vj, V are smooth real- valued functions of x 6 M 2 and (g jk ) is positive-definite 
matrix, < h ^ 1 is a Planck parameter and fi ^> 1 is a coupling parameter. I assume 
that A is a self-adjoint operator. 

*Work was partially supported by NSERC grant OGP0138277. 
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Now I consider 4-dimensional case, assuming that magnetic field is generic degenerated, 
i.e. of Martinet-Roussarie type [Ma, Rou] which will be described in details in subsection 
1.1. Here I just mention that in this case magnetic field F (which currently is considered 
as a closed 2- form) degenerates on the manifold T of dimension 3 and dimKerF = 2atZ. 
Further, dim(Ker F n TT) = 1 at T \ A and Ker F C TT at A where A C I is 1-dimensional 
manifold; furthermore, an angle between Ker F and TT. is exactly of magnitude dist(x, A). 
Finally, in an appropriate coordinates T = {xi = 0}, A = {xi = x 3 = x 4 = 0} and magnetic 
lines ^ G Ker F n TT are circles {xi = 0, x 2 = const, xf + xf = const}. 

My goal is to find asymptotics of J e(x, x, 0)^(x) dx with respect to h, /i, where e(x, y, t) 
is the Schwartz kernel of the spectral projector of A and ip is a smooth function supported 
in the vicinity of T. I assume that fih < const (otherwise e(x, y, 0) = 0(n~°°)). 

In the next paper I am gong to analyze the case of the generic non-degenerated field. 

0.1 Assumptions, notations and results 

So, magnetic field is characterized by Fj k = djV k — d^Vj. 

However, from the point of view of the classical and quantum dynamics and spectral 
asymptotics properties of (F J k ) = {g^){F lk ) are more important than those of (F jk ). Let 
i/Tx, ±/£ be eigenvalues of (F J k ), fj > 0. Then with the correct choice of notations fi x 
dist(x, T) and f 2 x 1. 

My first statement holds almost without any further assumptions: 

Theorem 0.1. Let F be of Martinet-Roussarie type and 
(0.2) V > e > 0. 

Let ip be supported in the small vicinity ofT. Then 

(0.3) | J (e(x, x, 0) - £ MW (x, 0))V(x) dx\ < C^ l / 2 h^ + C^ 2 h~ 2 

where 

(0.4) S MW {x,r) = 

(27r)-y/T 2 9 ( 2r+ V ~ ( 2m + l )^ hf i ~ {2n+l)fMhf 2 ) hh^fg 

(m,n)eZ+ 2 

is Magnetic Weyl Expression (see [IvrS] for the general case), g = det(g-' k )^ 1 . 
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Remark 0.2. Without any additional conditions remainder estimate (0.3) is sharp. An 
extra term 0(fi 2 h~ 2 ) in its right-hand expression is due to the possibility that for some 
pair (m, n) £ Z +2 function V — (2m + l)fihfi — (2n + l)fihf?_ identically vanishes in the 
vicinity of some point of E. However it is not as bad as it could be in the case of constant 
comeasurable f\, f-i in which case V — (2m+l)fihfi — (2n+l)fihf2 could vanish up to x 
pairs (m, ri) £ Z +2 and the remainder estimate would be the worst possible 0(fih~ 3 ). 

In the next few statements I am improving remainder estimate (0.3). 
Theorem 0.3. Let F be of Martinet-Roussarie type, conditions (0.2) and 

(°- 5 ) l(^) -(2n + l)/i/7| + |V L (-^)| <e ^ |detHess z (-^)| >e Vn £ Z+ 

'2 '2 '2 

be fulfilled. Let ip be supported in the small vicinity ofH. Then 

(0.6) | J (e(x, x, 0) - £ MW (x, 0))V(x) dx- J S^(x', r)^(x') dx'\ < C^ 2 h~ z 

where dx' is the standard density on Y. and£^„(x', r) is a correction term (see (4-51) which 
is 0( / u 5 / 4 /)- 3 / 2 ). 

Remark 0.4. (i) Correction terms and are of completely different nature; 

(ii) The first term in the left-hand expression of (0.5) is meaningful only in the case of the 
very strong magnetic field fih > e; so as fih < e condition (0.5) is equivalent to 

(0.7) |V z ( T )|<eo =^ |detHess z (-)| > e ; 

'2 '2 

(hi) One can weaken condition (0.5) to 

(0.8), 3n£Z+: |- - (2n + l)f 2fJl h\ + |V E (-)| < e => 

'2 '2 

Hess z (— ) has at least q eigenvalues with absolute values greater than e 

'2 

with q = 1,2,3 depending on the magnitude of \i\ as fih < e it is equivalent to 
(0.9), |V r (^)|<e =► 

'2 

Hess^ — ) has at least q eigenvalues with absolute values greater than e ; 

To 
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obviously (0.5), (0.7) are equivalent to (0.8) 3 , (0.9) 3 respectively; 

(iv) Under stronger condition 

(0.10) |(-) - {2n+l)fMh\ + |Vr(-)| > e Vn e Z+ 

the proof of theorem 0.3 is much easier; as fih <C 1 this condition is equivalent to 

(O.H) |V r (^)| >e . 

'2 

(v) In the very special cases (see Appendix A) one can derive remainder estimate (0.6) 
without condition (0.5) or even (0.8)i. Then the correction term could be larger, up to 
0(^1 A h^' 2 ). 

I will also prove 

Theorem 0.5. Let F is of Martinet-Roussarie type and condition (0.10) be fulfilled. Then 
as ip is supported in the small vicinity of T. H {V = 0} asymptotics (0.6) holds with 

«V,r) = o. 



0.2 Plan of the paper 

First of all, in section 1 I study the weak magnetic case /i < h~ s with sufficiently small 
exponent § > and derive remainder estimate 0( / u~ 1 / 2 /?~ 3 ) as the main part is given by 
the standard Weyl formula. I also analyze there the geometry of the degenerate magnetic 
field and the corresponding classical dynamics. 

Then, assuming that magnetic field is not weak, h~ 5 < n < ch^ 1 in section 2 I consider 
the different canonical forms of the operator in question; this canonical forms contain 
powers of yT 1 . There are universal canonical forms and also specific canonical forms as 
dist(x, T.) > /i~ 1/2 or as dist(x, A) > y 1/2 . 

Then in section 3 I derive remainder estimates 0(/i _1 / 2 /)~ 3 + fi 2 h~ 2 ) in the general 
case and 0(y~ 1 ^ 2 h^ 3 ) as /i > br 2 ^ and some non-degeneracy condition (depending on 
the magnitude of n) is fulfilled. The main part of this estimate is given by the standard 
intermediate formula 

(0.12) h- 1 f (F t ^ h - lr x T (t)VuQ;) dr 

J — oo 

where u is the Schwartz kernel of e' h ltA , Q = I and Xr(t) equals 1 as \t\ < |T and 
vanishes as |t| > T, T is rather arbitrary from interval [Tq, T^\ and the main fight is to 
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make 7~o as small and 7~i as large as possible. More precisely, the main part of asymptotics 
is given by the sum of expressions (0.12) with Q = making partition of unity and 
T= T L e[T {L)0 , T W ij. 

After this in section 4 I calculate (0.12) in more explicit way. 

Finally, Appendices A are devoted to the analysis of some special cases. 

1 Weak magnetic field 

1.1 Geometry of degenerating magnetic field 

So let Fjk(x) = djVk — dkVj be components of the matrix intensity of magnetic field; then 
uj = k Fj k dxj A dx k = d(J2k Vkdx k ) is the corresponding magnetic 2-form. According to 
[Ma], the generic form in dimension 4 never degenerates completely and has rank 2 on the 
submanifold E of codimension 1; more precisely, let 

(1.1) E = {x: rankw(x) = 2}; 
then 

(1.2) For a generic 2-form uj T. is the smooth manifold of codimension 1 and if ±ifi, 
±if 2 are eigenvalues of the corresponding matrix (Fjk), f 2 > f\ > ^ then fe(x) > e and 
/i(x) x dist(x, E). 

Let us consider Ker F(x) H 7~ X E in points of E. At each point it can be of dimension 1 
or 2; let 

(1.3) A = {x G E : dim (Ker F(x) n 7~ X E) = 2}; 

according to [Rou] 

(1.4) For a generic 2-form uj A is the smooth manifold of dimension 1; moreover in the 
appropriate coordinates 

(1.5) uj = dxi A dx 2 + x 3 dx 2 A dx 3 + d(x±x 3 + x 2 x 4 — ^-x|) A dx 4 . 



) This is a temporary notation, corresponding to Euclidean metrics g jk . 
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One can rewrite (1.5) as 

" + x!))d{x 2 + ^x 3 x 4 ) + (x 2 + x 2 )x x - ^(x 2 + x 2 ) 
after transformation x 2 i— > x 2 — |x 3 x 4 I get instead 



(1.6) uj = dx 1 A dx 2 — x 4 c6c! A dx 3 + x^dx\ A c/x 4 + x 3 c/x 2 A dx 3 + x 4 c/x 2 A dx A + 

2 ( x i - + *4 )) ^3 A c/x 4 = 
(*i — ^r 2 ) A c/x 2 + r 2 c/x! A dO + 2(x x - ^r 2 ) rdr AdO = c/((xi - ^r 2 )dx 2 + (x x - ^r 2 )r 2 ^ 

where x 3 = rcosO, x 4 = rsin 9. In contrast to the original Roussarie 2-form (1.5) this latter 
form is obviously (x 3 , x 4 )-circular symmetric. 
Then 

(1.7) T = { Xl = 0}, A = {x! = x 3 = x 4 = 0} 

and magnetic lines defined by 

dx 

(1.8) — G Ker F(x) fl T X T, x G I 

dt 

are helices {xi = 0,x 3 = rcos6',x 4 = rs\n9,x 2 = const — r 2 6/2] (with r = const) winging 
around A. 

Further, away from A one can rewrite 

(1.9) uj = 2x 1 c/x 1 A d6 + dy A dz, y = x x - ^r 2 , z = x 2 - 2y9 

and therefore there I have just a direct sum of two 2-dimensional magnetic fields; the first 
one is a generic degenerating field and second one is nondegenerate. 

1.2 Classical dynamics 

Classical dynamics described by the Hamiltonian 

(1.10) a(x,0 = \(£g ik {x)PjPk ~ V), PJ = - M 

j.k 



1 WEAK MAGNETIC FIELD 



7 



corresponding to (0.1), is much more complicated than the geometry because it also depends 
on the metrics (g-* k ) and the scalar potential V; however it appears that it depends mainly 
on f\, f 2 which are eigenvalues of — i{F J k ) 2 ) (where (F J k ) = {g-' l ){Fik)) and V and even on their 
ratios rather than these functions themselves. It also depends on the partition of energy 
(see below). 

1.2.1 Let us start from considering few examples. 

Example 1.1. As Fjk, gj and V are constant (which is not the case in our article) M 4 = 
Ki © K2 where Ky here and below are eigenspaces of —(F J k ) 2 corresponding to eigenvalues 
f?, dimKy = 2. Then the kinetic part of the Hamiltonian splits into 3° = + a\ where 
3° are Hamiltonians on Ky and 3° are movement integrals and the classical dynamics splits 
into two cyclotron movements along elliptical orbits in Ky with the angular velocities fifj] 
the sizes of these ellipses are x (/x^) _1 (ay ) 1 / 2 . 

Example 1.2. (i) Let in the frames of the previous example V be variable. Then 3y are no 
more movement integrals but they evolve with the average speeds 3 - 1 O^" 1 ) and in addition 
to the fast cyclotron movements there appear also slow drift movements along Ky with the 
speeds of x // _1 3y |Vk ; ( V/fj)\\ these drifts are orthogonal to V^V/fj). So these drifts 
depend on the instant energy partition. 

(ii) As (g-i k ) and (Fj k ) and thus fj become variable, the spaces Ky rotate with the average 
speeds 0(/i~ 1 ) but the previous statements remain true. 

Example 1.3. At points disjoint from A one can assume that magnetic field corresponds to 
uj = x\dx\ A dx 2 + cfcc3 A dx^, I used 2-form (1.9) and redefined coordinates x 1 = (x 2 , x 3 , x 4 ) in 
the obvious way. Assume that the metrics in these coordinates is Euclidean g^ k = 5jk and 
that V = V(x 1 , x 2 ). Then a 1 = a\-\V and 32 = 3° are movement integrals as in example 1.1 
and the movement splits into two movements again: one of them is the cyclotron movement 
in (x3,x 4 ) and another one is described by Hamiltonian a\ in coordinates (xi,X2). It was 
studied in details in [Ivr6]. 

Then there are outer and inner zones Z out and Z- mn respectively; in the outer zone Xi 
should be much larger than the cyclotron radius (associated with a{) which is 0([i~ 1 \xi\~ 1 ) 
and this is the ^out = {| x i| > Cji 1 / 2 }. However now a\ is not necessarily disjoint 

from (which was an assumption in [Ivr6]) and I will need to take it into account and 
redefine zones. 

2 ) In contrast to footnote 1) this is the final definition. 

3 ) While the speed of spatial movement is x 1 the drift is much slower and one can calculate all functions 
of x in the instant cyclotron centers. 
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Still I can conclude from [Ivr6] that if a trajectory (x(t),£(t)) on an energy level < c 
starts from the point (x(0),f(0)) with x(0) G 6(0, \) n {x : |xi| = 7 > C 2 ^ 1 / 2 }, then 
x(t) G 6(0, H {x : C _1 7 < |xi| < C7} as |t| < e^iry 2 . Further the drift speed is 
0(// -1 7~ 2 ); this estimate is sharp as long as a? is disjoint from 0. So again I have the 
cyclotron movement with angular velocities x fi\xi\ and \i and the drift movement mainly 
along x 2 . 

On the other hand, if a trajectory (x(t),£(t)) on an energy level < c starts from the 
point (x(0), f (0)) with x(0) € 6(0, §) n {x : |xi| < C/i" 1 / 2 } then x(t) G 6(0, \) n {x : |xi| < 
C 2 ^" 1 / 2 } as |t| < e and the averaged propagation speed is 0(1); this estimate is sharp for 
the "typical" (in the heuristic sense) trajectory. In this case there are no separate cyclotron 
and drift movements along but the cyclotron movement along K% remains. 

Example 1.4. Assume that in some coordinates (xx, ... , x 4 ) 2-form u is defined by (1.6) and 
also gi k = Sjk- Let V = const. I am going to show that the conclusion of the previous 
example remains true (with the drift directed mainly along magnetic lines). Now however 
the angle between Ki(x) and T x a is of magnitude r and therefore to prevent cyclotron 
movements from hitting T. it is sufficient to assume that C / u~ 1 7~ 1 r < 7 and to keep r of the 
same magnitude along trajectory it is sufficient to assume that r > C/i~ 1 7~ 1 . Alternatively 
I can assume that r < Cfi~ 1, -f~ 1 and then it will remain this way and the first inequality 
should be fulfilled with r = C/i _1 7 _1 . So I get zone 

(1.11) Z out = {C/i^S^ 1 < r < e^ 2 } U {r < C^ 1 ^ 1 , Cfi^j' 1 < e/i 7 2 } = 

{ 7 > Cmax^ 2 / 3 ,/:" 1 / 2 ^ 2 )}. 

1.2.2 Let instead of rather heuristic arguments above apply precise calculations. With uj 
defined by (1.6) and gi k = 5j k 

(1.12) a=i(£? + £+W(xi.r)). 

W( Xl , r) = (£2 - //(xx - l -r 2 )f + r 2 (ie - M*i " \r 2 )r 2 ) 2 = 

Then and £2 are movement integrals and the dynamics is restricted to the zone 3^ 2 ,§ fl = 
{(xx, r) : W(xi, r) < c}. While W is nonnegative and convex in the vicinity of (0, 0) it is 

a degenerate as £2 = £0 = 0: 1/1/ = /i 2 (l + r 2 )(xi — ^r 2 ) 2 and so particle seems to move 
far away if time is unbounded. So I will need to bound time (but this bound will be large 
enough to sustain a nice remainder estimate). 
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Now I can consider rigorously the general case. 

Let k a = (k™) be eigenvectors of (F J k ) = {g^')(Fik) corresponding to eigenvalues if a 
(a = 1, 2); then k^ are eigenvectors corresponding to eigenvalues —if a : 

(i.i3) E F i k « = frJC E F i k « = 



I normalize them so that 

(1.14) e sj^kK = °- E = °' E ®*^ k ? = 2< w- 

y,/c j,/c 
where the second equations in (1.13), (1.14) follow from the first ones. 

(1.15) One can select k,- G C°° satisfying (1.13), (1.14) in domain ft = 6(0, l)n{|xi < e} 4 >. 

Really, this statement is true for k2; therefore Kj G C°° and (F^) transforms Kj G C°° 
into itself and is skew-symmetric with respect to (g J ' k ). Then in orthonormal real base on 

Ki F = ( ^ qJ with /i G C°° and one can select ki G C°°. 

Let 

(1.16) z Q (x,e) = E^(*K(*-£), 4(*.o = E k - m w^( x ^). 

m m 

Pm = £m ~ AtV m (x), {py, P/<} = -yUFjfc 

where the second equation follows from the first one. Then since 

(1.17) \pj(x, £)| < Ci on the energy levels : a(x, £) < c } 
I arrive to 

(1.18) {Z Q , Z^} = {Zt, Z+} = 0, {Z Q , Z+} = 2/>f a 5 Q/3 mod 0(1). 



A > Provided I redefine temporarily f\ so that f\jx\ > 0. 
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1.2.3 Further, note that corrected xy. 

(1.19) xj^Xy-Zi-^PV J2 PJkF > 



kl 



satisfy 

(1.20) {x>,} = o(^V 2 ), 

(1.21) {x' j ,x' k } = ^F k + 0{^ 2 1 -% 



Also note that 



;i-22) a°(x, d = f i J>Vp/ = \ (W + |Z 2 | 2 ) 



and therefore 

(1.23) {a , Z Q } ee {a, ZJ = -/>f Q Z Q , {a , Z\} = {a, Z*} = ,>f a zj mod 0(1) 
and 

(1.24) {a°,6«} = {a, 6 Q } = mod 0(1), 6 Q = \Z a \ 2 . 

More precisely, (1.18) are fulfilled modulo linear forms 5 - 1 and therefore {a , b a } is a cubic 
formal) an( ^ | a j s a cubic form^) plus a linear forrn^). Moreover, these cubic terms 
should contain at least one factor Z\ or Z\ and also at least one factor Z> or Z|. 

(1.25) Then correcting b a (a = 1, 2) by cubic terms^) plus linear terms^), both multiplied 
by /i -1 , one can eliminate from {a, b a } all cubic terms but those with Z 1 Z 2 Z^ and Z\Z 2 Z\ 
and all linear terms but those with Z l5 Z\. Then these remaining terms produce 0(/i _1 ) 
error. 

Further, correcting in the zone Z oyit = {Cfi' 1 ^ 2 < |xi| < e} by cubic and linear terms of 
the above type but with coefficients of the type "x^ 1 x smooth function " one can eliminate 
the rest of the terms but now instead of 0(/i _1 ) error I get 0(/i -1 x{~ 2 ) error: 

(1.26) In the classical dynamics one can correct b a = \Z a \ 2 by 0(fi~ 1 x^ 1 ) so that their 
propagation speeds would be 0(/i~ 1 x 1 ~ 2 ). Therefore as long as x(t) e ft and |xi(t)| x 7, 
\t\ < T the following is fulfilled: 

\Z a (t)\ 2 = |Z Q (0)| 2 + Ofc-y- 1 + A* _1 7" 2 T-). 



5 ) With respect to (Z\, Z\, Z2, Z\ ) with smooth complex coefficients depending on x. 
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Similarly, consider {a, xj} which is a linear form K One can correct xj in the same way 
arriving to 

(1.27) In the classical dynamics one can correct xj by 0(yU~ 1 x 1 ~ 1 ) so that their propagation 
speeds would be O^^xf 2 ). Therefore as long as x(t) G 0. and |xi(t)| x 7, \t\ < T 

xj(t) = xj(o) + o{n~ x r l + ^t 2 t). 

Moreover, (1.24) and (1.25) with b a replaced by xj corrected yields that 
dx 

(1.28) — G Xi(x) + O^xf 1 ) for x, corrected modulo 0{fT V )• 

Still it is not yet what I want which is to calculate {a, xj} modulo 0(fi~ 1 x{ 1 ). However, let 
us consider first 

(1.29) b' 3 d ^ -Re(/>- 1 {4.^i}0M" 1 2'x)+^. <f>=h/*i 

and note that \Z\Z\, b' 3 } = (^(/i" 1 ). Also note that {a, 63} is a linear form of Z2, Z\ modulo 
0(fi~ 1 ) and therefore it could be corrected by a linear form of Z2, Z\ multiplied by so 
that after correction {a, b' 3 } = 0(1). Thus 

(1.30) In the classical dynamics one can correct b 3 = x\ by 0(//~ 1 ) so that its propagation 
speed would be 0(yU~ 1 ). Therefore x\(t) = x 2 (0) + O^yT 1 + pT x T\ and statements (1.26), 
(1.27) hold without precondition "as long as x(t) G ft and |xi(t)| x 7" which is fulfilled 
automatically for x(0) G 6(0, §) n {57 < |xi(0)| < 27}, T = e/^7 2 , 7 > 71 = C/i" 1/2 . 

(1.31) Also, trajectory originated in 6(0, |) n {|xx| < 7) remains in 6(0, 1) D {|xi| < C7} 
as 7 = 7i, T = e. 

Further, as long as |xi| < C7 and |{Zi,xi}| < cr, correction is 0(/i -1 (r + 7)) and 
therefore x 2 (t) = x 2 (0) + 0(fi~ 1 (r + 7) + /i" 1 "/"). Then I arrive to 

(1.32) As long as |{Zx,Xx}| < cr along trajectory, statements (1.30), (1.31) hold with 
71 = C/i~ 1 + C/i~ 1 / 2 r 1 / 2 . Moreover |{Zi, X!}(t)| < cr provided it was fulfilled (with constant 
c/2) at t = and 7 = e/i7 2 r, r > /i^ 1 / 2 + 7. 
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1.2.4 Now let us return to calculation of {a,xj} modulo 0(/i~ 1 ) where I will correct 
xj. Any error, larger than this, comes from eliminating linear terms containing Z\ or Z^, 
namely Re({Z|, Xj}Z\\ . This expression was eliminated by adding fi~ 1 f 1 ~ 1 Re(i{Z 1 ,xj}Zl) 
to xj which in turn generates an error 

(1.33) A *~ 1 Re/({a,C 1 {Zi, J& -}}4 + /i- 1 {4.J5}(/{a.^}-^i)). 

In the expression (1.33) the part containing factors Z 2 and Z\ in the combinations Z 2 , Z|, 
Z-j, Z 2 ' can be eliminated by the same way as above; the corresponding correction has an 
extra factor 0(/i _1 ) and the corresponding error just acquires factor 0(fx~ 1 x^ 1 ); therefore 
resulting error will be 0(/i~ 2 x 1 ~ 3 |Zi| + /i~ 2 x 1 -2 ) = 0(/i _1 + /i _1 7 _1 p) provided 

(1.34) \Zi\<cp, 7>c/i^ 1/2 . 

This leaves us with expression (which I do not call an error anymore) 

(1.35) i (-^VfH^J. *i} Re ('"{ z i. */} z i) + 0~ V 1 Re(«{Zi, x y })|Z 2 | 2 ) 
with the complex- valued coefficient a = a(x') provided 

(1.36) h = 0(x')x! + 0(x 2 ), x' = (x 2 , x 3 , x 4 ). 
Note that the first term in (1.35) is equal to 

(1.37) - ^xr 2 ^ 1 ■ Re({Zi, Xi}Zj) ■ Re(/{Z!, x y }Z/) = 

and one needs to calculate coefficients at Z 2 , |Zx| 2 as xx = only. Thus I need to calculate 
ki as Xx = 0. Let us decompose Fjk, gjk and k\ into powers of Xi; then F/^ki = /^ki implies 
that 

(1.38) F>? = 0, (F/ fc - + F? k k{ = 0, 
where Fj k corresponds to symplectic form (1.6) and therefore 



( 


1 


-x 4 


x 3 




-1 





X3 


x 4 




x 4 







2(x 1 -|(x| + x 4 2 ; 


) 


V"*3 


-x 4 


-2(x!-i(^+^); 
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Then (1.38), (1.39) imply for the "differentiation" part of Z\ 

(1.40) Zl, diff | Xl = = a([r l d e - rd 2 ) + p(rdx + <9 r )) 

with c- 1 < \a\ < c, c- 1 < Rei(3 < \(3\ < c. 

Then 

Re i({Z}, x 1 }Z 1 ) | xi=Q = 2\a\ 2 Re(i(3) (d e - r 2 d 2 ) 

is the differentiation along magnetic lines. Further, without any loss of the generality one 
can assume that a = 1. 

1.2.5 Consider first ip = Xi — r 2 /2 which is a regular function of x; note that Z\ij) = O(xi); 
therefore ip can be corrected by (^(/i" 1 ) so that ^ip = 0(/i _1 ) and therefore original ip 
is preserved modulo 0(/i _1 (T + 1)). Since I already know that xi is preserved modulo 
0(// _1 7~M 7 ~ + *)) I conclude that 

Proposition 1.5. Consider classical dynamics in ft = 6(0, 1) D {|xi| < Ce}, originated in 
a=B{0,±)n{\x 1 \<e}. 

(i) If in the original point |xi| = 7 with 7 > Cmaxf/i^r 1 , r 2 ) then both |xi(t)| and r(t) 
remain of the same magnitudes as T = e/ryr. 

(ii) If in the original point r 2 > 7 > C/i^ 1 / 2 /" 1 / 2 then both |xx(t)| and r(t) remain of the 
same magnitudes as T = eprfr . 

(in) If in the original point 7 > Cji~ 2 ^ and r < c / u~ 1 7~ 1 then |xi(t)| remains of the same 
magnitude while r(t) < 2c[i~ lr y~ 1 as T = e. 

(iv) If in the original point |xi| < 7 = c// _1 ' 2 r 1 ' 2 and r > C/i^ 1 / 3 then r(t) remains of the 
same magnitude while |xi(f)| < C7 as T = e. 

(v) If in the original point |xi| < 7 = c^T 2 ^ and r < cpT 1 ! 3 then r(t) < CpT 1 ^ and 
l x i(0l < Cj as T = e. 

On the other hand, at the energy level r one can rewrite the last term in (1.35) as 
/i~ 1 x ] ^ 1 Re(a{Z lt xj})(t + V) modulo C^/i" 1 + /i~ 1 7~ 1 p 2 ). Consider (1.38). Using the same 
technique as before one can eliminate 

fi-'x^r 1 ■ Re/(2{Z 1 t ,x 1 }{Z 1 t ,x / .}Z 2 ) 

by n~ 2 Xi 3 (fi~ 2 ■ Re^{Z/, xi}{Z|, XjjZ 2 ^ correction (which is 0(/i~ 2 7~ 3 rp 2 ) = 0(/i~ 1 7~ 1 p 2 )) 

as 7 > /i _1 / 2 r 1 / 2 adding 0(fi~ 2, y~ /[ r 2 p 2 + ji~ 2 ^~ 3 p 2 ) error; this error is less than e/i _1 7 _2 rp 2 
as 7 > C y u" 1/2 r 1 / 2 . 
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Also, any Z\ or Z\ unbalanced in the product could be treated in the same way as 
before. 

Then I arrive to 

Proposition 1.6. As 7 > Cp" 1 / 2 /" 1 / 2 the propagation speed with respect to 9 is 

(1.41) v(p; x, p" 1 ) = /i-V^x, p 2 , p) (p 2 + r~ 1 K 2 (x, p 2 , p) 7 ) + (p~ W^) 

with bounded and k disjoint from 0. 

I also need to consider dynamics of 6 2 = | -Z2 1 2 more precisely. Note that 

{|Z 2 | 2 ,a} = ^Re(({Z 2 ,Z 1 }Z 1 t + {Z 2 ,Zt}Z 1 -{Z 2 , \Z})zj) 

is a combination of cubic and linear terms (see footnote ^)) and the only terms where factor 
Zi is not compensated by Z\ and v. v. come from 

(1.42) {Z 2 , ZJ =a n Z 2 + a 12 Z 2 f + faZ x + ^Zj, 

(1.43) {Z 2 , Z^} =a 21 Z 2 + Q22Z} + faZi + &2Z+; 

these terms are equal to Re(cxZi)|Z 2 | 2 with a = a u + a 2 i- I need to calculate /3 as xx = 0. 
All other terms could be corrected by 0(p _1 ) leading to 0(p _1 ) error while this one is 
corrected by p _1 x 1 ^ 1 Re(//3Z 1 )|Z 2 | 2 leading to 0(p _1 Xj~ ) error; I will calculate this latter 
error more precisely below. Since 

«ii = (2/pf 2 )~ 1 {{Z 2 ,Z 1 },Zj}, 
«2i = (2/pf 2 )" 1 {{Z 2 ,Z 1 t },Z 2 t } 

calculated as Z\ = Z 2 = 0, I conclude that 

a = (2/pf 2 )- 1 (-{{Z 2 , Zi}, Zj} + {{Zj, ZJ, Z 2 } = -(/^{{Za, Zj}, Z x }. 

Therefore correction is 0(p -1 7 _1 p) and then the part of an error, with factors Z 2 and 
Z 2 balanced is 0(p~ 1 7~ 1 + p _1 7 ~ 2 pr). 

Moreover, the part of {|Z 2 | 2 , a} where Z 2 and Z 2 are balanced is k|Z 2 | 2 with 

K = (2p/ 2 )- 1 Re({/{Z 2 ,Z 2 t },Z 1 }| Zi=Z2=0 Z 1 t ) 
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since /'{Z 2 , Z 2 } is real; further, this expression is equal to 



(2f 2 r i {/2,i^ii 2 } + o(p 7 ) + o(/i- i 7- 1 ) 



(2f 2 )- 1 {f 2 , 3} - {2f 2 )- 1 {f 2 , |Z 2 | 2 } + 0(p 7 ) + OO^S" 1 ). 



Note that f 2 1 {f 2 ,\Z 2 \ 2 } has unbalanced Z 2 or Z\ and thus is eliminated by 0(fi 1 ) 
tion. So I arrive to an error — (2f 2 )^ 1 {f 2 , a}\Z 2 \ 2 . 



correc- 



Thus 



(1.44) (i) One can correct |Z 2 | 2 by 0(fi 1 + /i lr y 1 p) term so that the propagation speed 
after correction would be 0(/i _1 7" 1 + /i~ 1 7~ 2 pr). 

(ii) Moreover, one can correct f 2 _1 |Z 2 | 2 by (/i^ 1 + /i~ 1 7~ 1 p 2 ) term so that the propagation 
speed after correction would be 0(/i _1 + /i _1 7~ 2 p 2 ). 

I will also need the shorter term classical dynamics result: 

Proposition 1.7. Let us consider dynamics described in proposition 1.5(i)-(v). Then 
dist(x(t), (x(0)) > ept as \t\ < T x = e min(/i~ 1 7~ 1 , p). 

Remark 1.8. It effectively makes 7~i = e/i _1 7~ 1 because contribution to the remainder 
estimate of the subzone {|xi| x 7, \Z±\ < p, p < / u~ 1 7~ 1 } is 0(ph~ 3 x /i~ 2 7~ 2 x 7) = 
0(p _1 /7~ 3 7~ 1 ) which after summation over zone {7 > H~ 1/2 } results in 0(/i- 1/2 /r 3 ). 

1.3 Quantum dynamics. I. Outer zone 

In this subsection begin to deal with the quantum dynamics and the new parameter h 
appears; so one can now compare \i and h. Here I assume that magnetic field is rather 
weak; more precise requirements should vary from statement to statement but actually I 
do not care very much now and will just assume as needed that 



with small enough 6 > 0; I do not want to push it up so far. 

1.3.1 As I have mentioned there are few time scales and the shortest one is x /i" 1 . Let 
us deal with it first. Scaling x 1— > fix, t fit, h 1— > fih, [i \— > 1 I find ourselves with the 
standard Schrodinger operator with the propagator which is a standard Fourier integral 
operator. Returning to the original scale I get 



(1.45) 



C < n < h- s 
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Proposition 1.9. (i) As \t\ < 7~o = Cp 1 propagator e' h lAt is h-FIO corresponding 
to the classical dynamics <t> t = e tHa with Hamiltonian a(x,£). (ii) In particular, if ip £ 

Q°°(fl(x,e 7 )) 

(1.46) fnin-H-^A^ = e 2ni^h-^B^ 

with fz = 6(x) and the standard h-pseudodifferential operator 7 _1 B . 
Proof. Let us consider classical dynamics 4> t first; in (x, p) variables it is described by 
Hx- 

(1-47) ^/^E^MP/o 



dt 

;i.48) ^ = E + ^ E w*)p*p/ + a* -1 A/M 



(after rescaling t i— > with uniformly smooth /<*, /.*. Therefore in these coordinates 

t : (x, Z) i— > (x + p _1 X(x, p, t), V(x, p, t)). Thus ^ is a product of FIO corresponding 
to the symplectomorphism 4> t (a standard quantization, phase function is defined in the 
standard way and the symbol is just 1) and some PDO. 

In particular one can see easily that <t> t — / = 0(7) as t = 2nf 2 ^ 1 ; therefore I conclude 
that <t> t = e yHb with some symbol b = b(fi~ 1 x, Z). Quantizing it I get that (1.46) holds 
but with an extra PDO factor Q between FIO and i/j. However then one can perturb B by 
operator 7 _1 /?£>i so that 



e 2nifj,- 1 h- 1 f- 1 B q — e 2Trifj,- 1 h- 1 ff 1 (B+' r - 1 hB). 



I leave this easy exercise to the reader. □ 
Proposition 1.10. Let condition 

(1.49) V > e 
be fulfilled. Then 

(1.50) \F t -*h-WXT{t)r(i»l>Q)\ < Ch s Vr : |r| < e 

as p > K, T = f Ch\ log h\ < T < 2n(l — eo)p~ 1 ^~ 1 , 7 < e where ip is supported in B(x, 7), 
7 is calculated in x, Q is arbitrary h-pdo 6 ) . 

6 ) In our usual manner here and until the end of the paper x G C£°{[~ 1. — 5] U [5, 1]), X Q > °([— 1, 1]) 
are functions of [Brlvr] type, Xr{t) = x{t/T) etc. and T x (u) = u(x,x), fu = J u(x,x)dx. 
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Now let us consider in the outer zone the intermediate scale dynamics: 
Proposition 1.11. Consider x G Z out = f {7 > C/i -1 / 2 }. Then as \t\ < C/i~ 1 7~ 1 e'" ™" % 



h-FIO corresponding to the classical dynamics <t> t = e 3 w«£/i Hamiltonian a and 
(1.51) e* -1 ^ = e ih ' lt " A e ih ' ltlB ^ 

where t' = lixp^l^n, n = [(In)- 1 ? fit \, t" = t - t' . 

Proof. Formula (1.51) follows from (1.46) and the fact that |xx| retains its magnitude in 
the classical dynamics in the zone Z out . □ 

Proposition 1.12. (i) Let p < eh' 1 ] log/?!" 1 , x E Z out and Q be h-PDO with the symbol 
supported in {\Z\\ > p] with 



1.52) p > Cfi' 1 ^ 1 + C(/i/7 7 | log h\) 



1/2 



Then (1.50) holds as Cp 2 h\ log h\ < T < T\ d = e\p 1 7 1 . 

(ii) In particular, let 

(1.53) p > Cp- 1 ^ 1 + C(ph\ log h\f/ 2 . 

Then estimate (1.50) holds as T = Ch\ \ogh\ < T < eip' 1 ^ -1 . 

Proof. Let us temporarily direct K4 as a coordinate plane 7 -* {xi, x 2 } and consider p-admissible 
partition in (^1,^2)- Then one can see easily that the propagation along Ki has the speed 
x p. Further, condition (1.52) ensures that p retains its magnitude since the propagation 
speed of |Zi| 2 does not exceed Co/i~ 1 p 2 7 -2 + Cop" 1 and \Z\\ 2 is corrected by 0(pp^ 1 ). 
Furthermore, the last term in (1.52) ensures the logarithmic uncertainty principle. 

Then for time T the shift along Ki is exactly of magnitude pT and the logarithmic 
uncertainty principle pT x p > Ch\ log h\ is fulfilled as long as T > T' = Cp~ 2 h\ log h\. This 
implies (i). 

Now it follows from (i) and proposition 1.10 that estimate (1.50) holds as T = Ch\ log h\ < 
T < 7"i = eyU~ 1 7 ~ 1 provided V < e/i" 1 which means upgrade of condition (1.52) to 
(1.53). □ 



7 ) For T < cfi 1 7 1 and 7 > C\i x l 2 with C — C(e, c) our propagation is confined to B(x, ej). 
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1.3.2 On the other hand one can see easily that the contribution to the remainder estimate 
of the rather thin subzone where |xi| x 7 and condition (1.53) is violated does not exceed 
Cyu/j -3 (/i~ 2 7 -2 + fj,h\ log/j|)7; further, summation of this expression over Z out results in 
C/i -1 / 2 /?" 3 + Cp}hr 2 \ log /? I . Therefore I have proven 

Proposition 1.13. Contribution to the remainder estimate of the subzone of Z out where 

condition (1.52) is violated is O^/i -1 / 2 /? -3 + fi 2 h~ 2 \ \ogh\j. In particular this contribution 

is 0(p^ 2 h~ 3 ) as p< C{h\ \ogh\)- 2 / 5 . 

Actually, one can get rid off logarithmic factors here but under condition (1.45) it is not 
needed; this estimate (without logarithmic factor) is the best possible if no non- degeneracy 
condition is assumed. 

From now on I will consider only the main subzone of Z out where condition (1.53) is 
fulfilled and thus estimate (1.50) holds with T < T < Tj. I need to estimate its contribution 
to the remainder estimate. 

First of all one needs the quantum version of the results of the previous subsection. 

Proposition 1.14. Let condition (1-45) be fulfilled. Let function ip be properly supported 
in B(x,y), and Q be properly supported in {\Z\\ < p] with p > Qj/i -1 ^ 2 and 

(1.54) i > Cp^h\ \ogh\. 

Further, let 7 > CqP" 1 ^ 2 and T < T 2 d = e/i7 2 and let x" = (x 3 , x 4 ). Then 

(1.55) r = |x"| > r> = Co/i'S" 1 + Qi 
then estimate 

(1.56) |F t ^- lT xr(t)(l - il>i){wl> y Ql)\ < Ch s 
holds provided ipi = 1 in the domain 

(1.57) jx : C - X 7 < |xi| < C 7, C l r < \x"\ < C r, 

dist(x,x) > C 7'+ Cop' 1 !' 1 + C Q pT x {r + 7 ) 7 " 2 7}. 

(ii) On the other hand, if condition (1.55) is violated then estimate (1.56) holds provided 
tpi = 1 in the domain 

(1.58) {x : C _1 7 < |xi| < C 7, \x"\ < C r', 

dist(x,x) > C 7'+ C ^ 1 7 _1 + Q^ 1 ^' + l)T 2 T^. 
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Proof. Proof is the standard one (see e.g. [Brlvr]) based on the analysis of the symbol 
(1.59) w ^L( t±el/ -i(x(x,0-X(y,0))) 

where w is the same function as in [Brlvr], X(x, £) is one of the corrected symbols x 2 , 
x i — |l x "| 2 an d x j with j = 1, ... , 4 and i/ is the corresponding speed, namely p -1 , p -1 and 
pT x x[ 2 r respectively. 

Obviously, it is sufficient to consider T such that the last term in definitions of (1.57), (1.58) 
is dominant: fi~ 1 (r + / ~f) / ~f~ 2 T > 7' + yU _1 7 _1 . One can see easily that under this condition 
and (1.45) symbol (1.59) is quantizable; here also one can be more specific about exponent 
5 > in condition (1.45). □ 

Now I want to be more precise for smaller r, p: 

Proposition 1.15. In frames of proposition 1.14 as T < ep'y p estimate 

(1-60) \F t ^ T xr{t){l - QJM^yQy)] < Ch s 

holds provided Qi = I in the domain {c _1 p < \Z\\ < cp} and estimate (1.56) holds provided 
V>i = 1 in the domain 

(1.61) {x : Q _1 7 < |xi| < C 7, C^r < \x"\ < Qr, 

dist(x, x) > Q7' + Co/i-^" 1 + C p^~ 2 {rp 2 + 7) r| 

as r > r' , and in the domain 

(1.62) |x : Cq 1 -/ < |xi| < C 7, \x"\ < C r', 

dist(x.x) > Co7' + C /i~ 1 7^ 1 + C ^ 1 7~ 2 (r'p 2 + 7)^}- 

as r < r' . 

Proof. Proof is standard, based on the symbol (1.59) which is quantizable due to assumption 
(1.45); here again X is one of the corrected symbols \Zi\ 2 , x' and v is the corresponding 
upper bound for a speed, namely / u _1 7 _2 p, /i~ 1 7~ 2 (p 2 + 7) respectively. □ 

Note that under condition (1.45) r' x /i~ 1 7~ 1 and I want to eliminate the corresponding 
subzone from the future analysis. Really, let us note that the contribution of the subzone 
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{7 < < 27, \x"\ < r,\Zi\ < p} to the remainder estimate is 0( T^ 1 h~ 3 r 2 p 2 'y) = 
0(ph~ 3 r 2 p '1 ) (since I already eliminated subzone where condition (1.53) is violated) and 
plugging rp = p _1 7 _1 results in 0(p~ 1 h~ 3 ). Then summation with respect to 7 results in 
0(p~ 1 h~ 3 \ logp|). So, I conclude that 

(1.63) Contribution to the remainder estimate of Z out n {rp < C / u _1 7~ 1 } is o(p^ 1 ^ 2 h^ 3 ). 

1.3.3 From now on I restrict myself to the analysis of subzone Z out n {rp > Cp~ 1, -f~ 1 }. 

Now one can see that in Z out function v(p; x, p -1 ) defined by (1.41) has no more than 
one root p > Q/~ 1 r -1 7~ 1 ; I denote it by w = w(x, p^ 1 ). Otherwise (if such root does not 
exist) I set w = 0. One can see easily that 

(1.64) v{p;x,p' 1 ) x/i"V 2 (p 2 - w 2 ). 

If p — w x A then the propagation speed for corrected symbol 9 is v(p, .); then the linear 
shift for time e / u~ 1 7~ 1 is x r/i~ 2 7 _3 pA (where r appears because 9 is an angle). This shift 
is observable if the logarithmic uncertainty principle r/i~ 2 7~ 3 pA x A > Ch\ \ogh\ holds 8 -' 
i.e. 

(1.65) A > Cp(h\ log/i|) 1/ V 1/ V /2 ''" 1 - 

Since the propagation speed of symbol |Zi| 2 corrected is 0(/i _1 7~ 2 p 2 ), the magnitude of A 
is preserved during the time interval T" = e/i7 2 p _1 A which is larger than C/i~ 1 7~ 1 as 

(1.66) A > CpT 2 p-^- 3 . 

This condition is stronger than (1.65) as r > p~ x and (1.45) is fulfilled with small enough 
5 > 0. 

To justify this assertion properly one needs to operate with the quantizable symbols 
and one can see easily that this is the case under assumptions (1.45) with small enough 
5 > and (1.65). Then the contribution to the remainder estimate of the corresponding 
subzone (where magnitudes of 7, p, r and A are fixed) does not exceed Cr 2 h~ 3 p , -fA/ T" = 

Qi^/rVpV 1 . 

Let us note that as magnitudes of 7, r, p are fixed and A ranges from C/i" 2 p~ 1 7 -3 
dictated by (1.66) to ep, the number of such elements does not exceed Clog p. So the 
summation with respect to this zone results in Cp~ 1 /?~ 3 r 2 p 2 7~ 1 log p. Further, as A > ep, 

8 ' Now A is the scale in p. 
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there are no more than C such elements and therefore the contribution to the remainder 
estimate of the part of the zone 

(1-67) y^ P ,r={\l<\xx\<2 1 , ^p<\Z 1 \<2p, \r<\x"\<2r) 

where condition (1.66) is fulfilled does not exceed Cp,~ 1 h~ 3 r 2 p 2/ ~f~ 1 logp as well. Summation 
with respect to p with p < Cipi/r) 1 / 2 results in Cp~ 1 h~ 3 r\ogp and then the summation 
with respect to all r, 7 results in /i _1 /j~ 3 | logp| 2 = o(p^ 1 ^ 2 h^ 3 ). So, contribution to the 
remainder estimate of the part of the zone Z out where r > fi~ 1/ -f~ 1 , p 2 < C7 and condition 
(1.66) is fulfilled is 0{p^l 2 h^). 

The same arguments one can apply as p 2 > C7/V but in this case A = ep and con- 
dition (1.64) is fulfilled automatically, there will be no factor log/i and the summation of 
C/i~ 1 /)~ 3 r 2 p 2 7~ 1 with respect to p, r, 7 results in CpT x l 2 h~ 3 . 

Therefore I arrive to 

(1.68) Contribution to the remainder estimate of the subzone of Z out where condition (1.66) 
is fulfilled is 0(/i- 1 / 2 / ) -3). 

Finally, contribution to the remainder estimate of the subzone y jiP with 7 > CpT 1 ! 2 
where condition (1.66) is violated does not exceed 

C(/T 3 p 7 x p-W 3 x Tf 1 ) = Cp-^h- 3 

since T\ = ep" 1 7~ 1 and summation with respect to 7 results again in 0(p _1 / 2 /)~ 3 ). 
Combining with (1.68) I obtain immediately 

Proposition 1.16. Under condition (1-44) the contribution to the remainder estimate of 
the zone Z out is 0(p~ 1 / 2 /7~ 3 ) . 

1.4 Quantum dynamics. II. Inner zone 

Now I need to analyze the contribution to the remainder estimate of the inner zone Z m „ = 
{\xi\ < 70 = C/i" 1 / 2 } under condition (1.45). Both classical and quantum dynamics are 
confined to this zone as T < e and it starts in B(0, |); propagation speed with respect to x 
does not exceed Q>. Note that the correction procedure still works as long as one eliminates 
only unbalanced factors Z2, Z\. 
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1.4.1 Contribution to the remainder estimate of subzone {\Z%\ < p, |xi| < 7} does not 
exceed CT -1 h~ 3 p 2 r 2 7 = C/i/?~ 3 7p 2 with T = T = eji~ l ] plugging 7 = cpT 1 ! 2 and p = 
cpT x l 2 I get 0(p~ 1 / 2 h~ 3 ) and therefore one needs to consider only the contribution to the 
remainder estimate of the subzone ^inn n{ | ^1 1 > CyU -1 / 2 }; but then proposition 1.12 remains 
true here and one can increase T = Tq to T = Ti = e/i -1 / 29 ). 

Then the contribution to the remainder estimate of zone {|Zx| < p, \x±\ < 7, \x"\ < r} 
does not exceed CT^ 1 h~ 3 p 2 r 2 j = Cp}l 2 h~ 3 ^r 2 p 2 with T = e/i -1 / 2 ; in particular, the total 
contribution of the whole inner zone Z- mn is 0(h~ 3 ) and one needs to recover factor p" 1 ^ 2 . 

To recover this factor one needs just to increase T — 7~i = epT 1 ! 2 to T = T 2 = e. 
Before doing this just note that the contributions to the remainder estimates of subzones 
Z inn n {\Z X \ < cp- 1 / 4 }, Z ]nn n {|x"| < cp- 1 '*} and Z ]nn n {\Z X \ ■ \x"\ < cp~ l l*\ \ogfi\~ 1 } are 
0(p>~ 1 / 2 h~ 3 ) and therefore 

(1.69) One needs to consider only the contribution to the remainder estimate of the subzone 

^inn H {|X"| > Cll T X '\ \Z X \ > Cl fi-y\ \Z X \ ■ \X"\ > C^-V^ \0gp\~ 1 }. 

I remind that the classical dynamics which starts at {|xi| < 7, |x"| = r > Ci/i" 1 / 4 } 
remains in {|xi| < C7, ||x"| — r < Ci/i" 1 / 2 /"" 1 } as 7" = e; actually it could be larger as 
pr < 1. 

Repeating arguments of the proof of proposition 1.14 one can justify it easily for quan- 
tum propagation as well: 

Proposition 1.17. Under condition (1-45) with small enough 5 > the quantum dynam- 
ics which starts at {Ix^ < 7, \x"\ = r > Ci/i" 1 / 4 } remains in {\xi\ < C7, ||x"| — r < 
Ci/i" 1 / 2 ^ 1 } asT = e. 

Now let us consider again function |x 2 ; then 

(1.70) ^{a,x 2 } = Re({Z 1 t ,x 1 }x 1 Z 1 + {Z 2 f , x^Za) 
and thus should be corrected by 

(1.71) /i- 1 Re/^.XiKf'xiZj + {Z 2 t ,x 1 }^ 1 x 1 Z 2 ) 
which is 0(yU _1 r) as r > p~ x l 2 and leads to the new error 

(1.72) p^ 1 Re /({a, {Zj, x 1 }ff 1 x 1 }Z 1 + {Zj, x^f ^({a, Z x } + />Zi)+ 

{a, {Z 2 t ,x 1 }/ 2 - 1 x 1 }Z 2 + {Zj.x^^XiHa, Z 2 } + />Z 2 )). 

9 ) As in footnote ®) one must take T\ — e min(/i~ 1//2 , p) but contribution to the remainder estimate of 
subzone Z mn n {|Zi| < Ai~ 1/2 } is 0{\ih~ 3 = CK/i- 1 / 2 /)- 3 ). 
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This new error is 0(/i 1 ) and thus taking T = er one can see that the oscillation of x\ 
would be 0(p~ 1 r) which leads to 

Proposition 1.18. (i) The classical dynamics which starts at {{x^ = 7, \x"\ = r > 
Ci/i" 1 / 4 } with CoyU -1 ^ 2 /" 1 / 2 < 7 < Co/i" 1 ^ 2 remains in {C _1 7 < |xi| < C7} as T = er. 

(ii) The classical dynamics which starts at {\x\\ < 7, \x"\ = r > Cip^ 1 ^} with 7 = 
Co/i _1 ^ 2 r 1 / 2 remains in {|xx| < C7} as T = er. 

(Hi) Under condition (1-45) with small enough 5 > statements (i),(H) remain true for a 
quantum dynamics as well. 

So, one can see that the dynamics in the zones described in proposition 1.18(i), (ii) 
are different with the dynamics in the former resembling the dynamics in Z out . Then the 
arguments of the previous subsection work perfectly in the zone described in proposition 
1.18(i); tedious but easy 10 '' details I leave to the reader: 

Proposition 1.19. Under condition (1-45) with sufficiently small 6 > the contribution 
to the remainder estimate of the zone {Cop^^-^lx"] 1 / 2 < |xx| < Cq/i^ 1 ^ 2 } does not exceed 

Thus I am left with the zone 

(1-73) {Ci^ 1/6 < |x"| < c, |xi| < d/i-^lx"! 172 }. 

Again the dynamics in the zones D^ 7 , p , r with pr < ep^ 2 and pr > ep'y 2 are very different 
with the former one more similar to the dynamics in the outer zone (as long as it remains 
there). 

Similarly the propagation speed of the symbol \Zi\ 2 corrected by 0(/i~ 1 7~ 1 r|Z 1 | + /i" 1 ) 
does not exceed C(yU _1 7~ 2 r|Zi| + /i -1 ) and therefore the propagation speed of |Zi| does 
not exceed C{uT x ^~ 2 r + /i -1 / 2 ) and thus magnitude of \Z\\ x p is preserved on the time 
interval T = e(/i~ 1 7" 2 r + p^ 1 ^ 2 ) 1 p as long as 7 > C(prp) 1 ^ 2 . As a result I arrive to 

Proposition 1.20. (%) The classical dynamics which starts at {\x\\ = 7, \x"\ = r > 
Ci/i -1 / 4 , \Z\\ = p] with CqP' 1 / 2 r 1 / 2 p 1 ! 2 < 7 < Cop' 1 ^ 2 remains in {C _1 p < |Zi| < Cp] as 

T = e min [r, pqrp~ x , p^p^ 1 J . 

10 ) Because one needs to worry only about eliminating factor p 1 / 2 while in the previous subsection the 
offending factor was p 3 ^ 2 . One can just take A = Cp^ 1 ^ 2 . Also one can assume that |xi| > 7 = cipr 3 ^ 5 ; 
otherwise \x"\ < r — C2p~ 1 ^ 5 and the contribution to the remainder estimate of this zone would be 
0(Ai/7- 3 7 2 r 2 ) = 0(n- 1/2 h- 3 ). 

But then p > pT 1 ^ -1 in the zone of interest and one can finally upgrade 7~i = e/i -1 / 2 to T\ — ep^ 1 ^ 1 
and the crude estimate to 0(ph~ 3 p 2 j 2 r 2 ). 
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(ii) The classical dynamics which starts at {\x\\ < 7, \x"\ = r > C\p, 1 ^ 4 , \Z\\ = p}, with 
7 = Co/i~ 1 / 2 r 1 / 2 p 1 / 2 remains in {C~ 1 p< \Z\\ < Cp] as T = e min [r, p^rp" 1 , /i~ 1 p _1 j. 

(Hi) Under condition (1-45) with small enough S > statements (i),(ii) remain true for a 
quantum dynamics as well. 

1.4.2 Then again arguments of the previous subsection work perfectly in the zone de- 
scribed in proposition 1.20(i); tedious but easy details (see footnote ^) with the obvious 
modifications.) I leave to the reader: 

Proposition 1.21. Under condition (1-44) with sufficiently small 6 > the contribution to 
the remainder estimate of the zone {Cqh~ 1 I 2 \x"\ 1 I 2 \Zi\ 1 / 2 < |xi| < Cop^ 1 ^ 2 } does not exceed 

c^i 2 h-\ 

Therefore I am left with the true inner zone 

(1.74) Zl^llx^C^lx'I^IZil 172 }- 
Similarly to (1.69) 

(1.75) One needs to consider only the contribution to the remainder estimate of the subzone 
ZL H {|x"| > c^-Ve, \Z X \ > d^-Ve, |Z X | ■ |x"| > c x yr^\ \ogp\~ 1 }. 

Let us consider classical dynamics in this zone first; I am interested in the time interval 
T < er. There 

{a, Z x } = {Zj, Z 1 }Z 1 + {Zl Z 1 }Z 2 + {Z 2 , Z t }Zl - {V, Z x } 
and then I correct Z\ to 

(1.76) Ci = Zi + ip^f^iA, Z 1 }Z 2 - i f i- 1 f 2 1 {Z 2 , Z y }Z\ 
satisfying 

{a, Ci} = {Zl Z 1 }Z 1 - ifi'p- 1 ({Z 2 , {Zj, Z 1 }} - {Zl {Z 2 , Z,}}) \Z 2 \ 2 - {V, ZJ = 

{zl z 1 \z 1 + //r V 1 ! 7 !. I 7 ?. z 2 t }}|z 2 | 2 - { v, Zi} = 

{Z 1 t ,Z 1 }Z 1 -/ 2 ~ 1 {Zi,/ : 2}|Z 2 | 2 -{\/,Z 1 } mod 0{p~ 1 ' 2 ). 
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Since I do not mind to change the speed to one of the same magnitude and I consider only 
level 0, I can assume with no loss of the generality that 

(1.77) h = 1; 

then finally {a, 0} = {ClCi - V, &} mod O^ 1 ' 2 ). 

Also corrected functions xj satisfy {a, xj} = {CxCi ~~ V, xj} mod 0(/i~ 1 / 2 ) and I arrive 
to the following conclusion: 

(1.78) Under condition (1.77) the classical dynamics in zone {|xi| < Cji" 1 ^ 2 } in variables 
xj, Z\, Z 2 is described modulo 0(/i _1 + /i~ 1 / 2 T) by the solution to the short system 

(1-79) ^F =1 2 {Z ^ Zl - V ' Zl} ' 

Hy- 1 

(1.80) dt = 2 {Z ^ Zl ~ V ' Xj} ' 

(1.81) {Z x , Z 2 j = O(^), {Z lt Zj} = O^" 1 ). 

Further, for each given initial point one can replace all functions modulo (xx — ), (x 2 — 

Or 2 ) so that I get two-dimensional system in variables (r, 9) described by the Hamiltonian 
^{Z\Zi — V) with Z\, Z\ satisfying commutator relation 

(1.82) {Z 1 ,Zl} = ^?(r--r) + 0(^ 1 )) 

which is exactly system studied extensively in [Ivr6] with magnetic field 0(r — r) and 
coupling constant Jx = fi?. Condition (1.74) means exactly that I am in the inner zone 
{|r — r\ < Cjl~ 1 / 2 \Zi\ 1 / 2 } for this system; since (in contrast to [Ivr6]) \Z\\ is not disjoint 
from anymore one needs to define inner zone in this way. 
Then the drift speed with respect to 9r is 

(1.83) x (77- k*(j)- 1/2 pjpr 

where constant k* « 0.66 is defined in [Ivr6] and 77 is the corrected symbol /ix 2 /2 which 
leads to its propagation speed given by (1.72) multiplied by /x. 

One can see easily that the fourth term there \Z\, xi}f 2 ~ 1 xi ({a, Z 2 } + //JZ2) is 0(/i -1 / 2 ). 
Further, the third term {a, {Z\, xx}/ 2 ~ 1 x 1 }Z 2 could be rewritten modulo 0(/i~ 1 / 2 ) and terms 
with unbalanced Z 2 , Z\ as |{Z 2 , x 1 }| 2 / 2 _1 |Z 2 | 2 and disappears after applying Re /'. 

So I am left with just two terms 

(1.84) Re /({a, {Z*, Xi}^ 1 }^ + {Zj.Xi^fta, Zx} + ,>Zi)) 

where = f\l x \ and one can replace a by Z^Zx — V (other terms have unbalanced Z 2 , Z|). 
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1.4.3 So I am exactly in [Ivr6] situation. I remind that in [Ivr6] symbol Zi could be 
reduced to 

(1.85) Z x = e'^W(6 + ia(x)(& - V 2 (x))) 

with real- valued a, V2 such that a = 1 + 0(xi), V 2 = x 2 /2 + 0(x 3 ). 

One can see easily that in this case that corrected px 2 would be £2 + 0(p,~ 1 ^ 2 ); I want 
to remind that in [Ivr6] exactly deviation of £2 from /r*0 _1//2 p was used as the measure of 
the drift with p = V 1 / 2 . Note that the potential V now should be replaced by (V — £2) 
with E2 = |Zf|. 

Therefore due to the logarithmic uncertainty principle the violation of periodicity is 
observable after the first turn as 

(1.86) (77 - /c*0- 1/2 p)prp~ 1/2 x (77 - /c*0~ 1/2 p) > Ch\ \ogh\; 

under assumption (1.45) the latter condition is the automatic corollary of 

(1.87) A = \ V - /c*0" 1/2 p| > Cyr^r^p- 1 ' 2 . 

Under this condition the drift speed for the "short" system is greater than Ceip -1 / 2 which 
is the larger than the error in the drift speed cp -1 / 2 (see (1.78). 

In this case Ti = erA and the rest is easy but extra logarithmic factor would appear; 
however considering the right direction one can take Ty = erp 1 / 2 A 1 / 2 and following [Ivr6] 
with the standard justification on the quantum level one can prove easily that 

(1.88) The contribution of the zone where condition (1.87) is fulfilled to the remainder 
estimate is 0(p~ 1 / 2 /? -3 ). 

On the other hand, contribution to the remainder estimate of zone where this condition 
is violated would not exceed 

p-^h-'r^p-^rdrdp") = 0(p^ 2 h- 3 ) 

as well. 

So, the main result of this section is proven: 

Proposition 1.22. Under condition (1-45) with small enough 5 > the remainder estimate 
is 0(p~ 1//2 /)~ 3 ) while the principal part is given by the standard Weyl formula (which is equal 
modulo the remainder estimate to magnetic Weyl formula) . 

This statement implies trivially 

Corollary 1.23. Under conditions (0.2) and (1-45) estimate (0.3) holds. 
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2 Strong magnetic field. Canonical forms 

Now I want to consider the main case 

(2.1) h-^fi^Ch- 1 

with an arbitrarily small exponent 5 > 0. In this section I consider different canonical forms 
of the Magnetic Schrodinger Operator in question (depending on the zone). 

2.1 Precanonical Form 
2.1.1 Let 

(2.2) (Mx,£) = ^*)- 1/2 Z 2 (x,£) 
Then 

(2.3) * -/> _1 {L/ 2 , U\} = 2 mod O k 

with k = 1 where O k = d?*(/i -1 Zi, [i~ x Z\, /i~ 1 Z 2 , /i~ 1 Zj) is the space of the sum of polyno- 
mials containing monoms of order k or higher with respect to (fi~ 1 Z 1 , {i~ x Z\, fi~ 1 Z 2 , /i~ 1 Zj) 
with the coefficients, smoothly depending on x. 

Note that correcting U 2 by U 2 O k one trades (2.3)* by (2.3)* + i and therefore after an 
appropriate correction (2.3) holds with arbitrarily large k = M. 

So, correcting U 2 : 

(2.4) U 2 ^U 2 + Yl a klpq {x)U*Ul'Z?Zl q v 1 - k -'-»-o 

k+l+p+q>2 
k>l 

one one can achieve 

(2.5) -ifi^{U 2 , U\} = 2 mod 0(^ M ) 

with arbitrarily large M. 

Note that {U 2 , Z{\ = {U 2 , Z{\ = modulo 0\ both originally and after correction (2.4). 
I claim that 

Proposition 2.1. Correcting 

(2.6) Z X ^Z 1+ a> klpq (x)U*U±'Z 1 p Zl<i^ k -<-i>-<<, 

k+l+p+q>2 
k+l>l 

(2.7) xj ^xj+Y a] klpq (x)U k Ul'Z?Zl V^'^- 

k+l>l 
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one can arrange 

(2.8) {U 2 ,Xj} EE {Lfe.Zi} EE {U 2 ,Z\} EE 0, 

(2.9) {*4.*/} = {^2- *i} = {^4- 4} = mod 0(/i- M ). 

Proof. Really, correcting in the way described (but with the sum over k > 1 rather than 
k + / > 1) I can achieve (2.9). However, so far corrected symbol Xj is not necessarily 
real-valued and corrected symbols Z\ and Z\ are not necessarily complex conjugate. 

Then due to (2.9), (2.5) and Poisson identity {U±,{U 2 ,y}} = with y = Xj, Z x , Z\. 
Therefore 

l+p+q>l 

{u 2 ,zi}= Yl /3UW^ t, zfz;v- , - p - £? - 

l+p+q>l 

{U 2 ,xj}= Y %, pq (x)U± l Z?Zl^- l - p - q 

l+p+q>l 

where x, Z\ and Z\ are already corrected symbols. Since at this moment Xj are not neces- 
sarily real-valued one needs to plug them formally into functions and use a Taylor decom- 
position with respect to [i~ x Z\, [i~ x Z 2 , [i~ x \J\, yT x \J 2 . 

Then one can correct xj, Z\ and Z[ according to formulae (2. 7), (2. 8) with summation 
k = so that (2.8) holds. 

Note that instead of Z-y and Z\ one could consider ReZi and Im Z\ and correct them 
deriving (2.8)-(2.9). At this moment corrected symbols xj, ReZi and ImZi are not neces- 
sarily real-valued. However U 2 and \j\ are truly complex conjugate and then (2.8)-(2.9) 
for corrected symbols Xj, ReZ 1; Im Z\ imply the same equalities for the real parts of them. 
Let us replace then (corrected) symbols xj, ReZi and \mZi by their real parts. After this 
(corrected) symbols xj become real- valued and (corrected) symbols Z\ = Re Z\ + / Im Z\ and 
Z| = Re Zx — /' Im Z\ become complex conjugate. □ 

2.1.2 One can rewrite operator in question in these new variables as 

(2.10) 2A ee f 2 {x)UlU 2 + Z\Z X - V(x)+ 

Y b klpq (x)U 2 k Ul l Z?Zl^ 2 - k -'-r-« + Y b> klpq (x)UZul'Z?Z}<^ k -i-i>-« 

k+l+p+q>3 k+l+p+q>3 
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modulo lower order terms where the first line is the "main part" and the second line is the 
"perturbation"; here the first part of the "perturbation" comes from ^2 + Z^Zi and the 
second part comes from — V. 

Proposition 2.2. For arbitrarily large M there exist e = e(M) > and a real valued symbol 

(2.11) c= e kPq n,(x)uZul'z?zl q f,- 1 - k - l ->-<- m 

k+l>l 

such that in the strip {|xi| < e} the following equality holds (so far only for the principal 
symbols of operators in question): 

(2.12) A* ^ e-'^-^Ae^- 1 ^ = ±(f 2 (x*)U*Uf + Zfzf - V(x#)) + 

£ ) ( u* up) k z# p zp y -^-p-,-2 m 

2k+p+q+2m>3 

where e~ lfl lfl lc , e ,At lh l£ are (formal) fi^h-FIOs and 

(2.13) xf = e-^ lb - lc xje ! ^ lh - lc , U* = e^^^e^ 1 ^, 

Z* = e- i »- ih - lc Z x e i *~ lh ~ lc 

are linked to the original symbols xj, U2, Z\ by formulae similar to (2.4), {2.6), (2.7) but 
containing also factors /i" 1 exactly on the same role as U2 or il\ so that k + I in the 
condition of summation is replaced by k + / + m ll > . 
All equalities here hold modulo 0(/i~ M ). 

Proof. Proof is standard (see [Ivr5]). The problem of resonance is avoided because I remove 
only terms b if (x)U 2 k U\ 'Zf Z\ q with k ^ I which leads to the denominator 

- f)h + (P ~ q)fi) 

n ) More precisely 

(2.4)* U 2 ^U* = U 2 + J2 aki Pqm (x)U*Ul'ZPzl q ^- k -l-P-<>- m , 

k+l+p+q+m>2 

(2.6) * Z 1 ~Z* = Z 1 + J2 a> klpqm (x)U*Ul<Z>Zl<^- k ->-P-'-'» l 

k+l+p+q+m>2 

(2.7) * XJ ^ X # =XJ+ J2 a% lpqm (x)U*Ut l Z?Zl''n- k - l ->>-«-">. 

k+l+p+q+m>l 
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in the correction term; this denominator is of magnitude fi(k — I) since k ^ I, f 2 x 1 and 
f x x |xi| < e(M). □ 

Now one can upgrade all the above arguments to operators so that equalities (2. 5), (2. 8) 
and (2.9) hold for commutators rather than for Poisson brackets and (2.12) is replaced by 

(2.12)* A* = e- i ' t ~ lh ~ lc Ae i ' t ~ lh ~ l£ = - [f 2 w U^ •U 2 + Z* • Z 1 - V w ^j + 

E B * Pqm s(U2 • U*)'Z[Z* ^2-2 k -2 m - P - q -s h s 

2k+p+q+2m+2s>3 

where a w is the Weyl quantization of symbol a and one should remember that the corrected 
symbol x is not just a function of coordinates but symbols and therefore the functions of 
x are replaced by the symbols of / u~ 1 /?-PDOs. Here and below • means the symmetrized 
product: K • L = \(KL + LK). So, I arrive to 

Proposition 2.3. For arbitrarily large M and arbitrarily small 5 > there exist e = 
e(M, 5) > and a real valued symbol 

(2.11)* C= £ i kpqms {x)U^Ul'ZPZl q ^ x - k -'-^- m - s h s 

k+l+p+q+m+2s>l 

such that under condition (2.1) in the strip {\xi\ < e} formula (2.12)* holds where e~' M 1/1 ^ , 
e iti h c are fo_pio s an d x # t \jjf \zf are full symbols still defined by (2.13) and they 
are linked to the original symbols xj, Z 1; U 2 by formulae similar to (2. 4)*, (2. 6)*, (2. 7)* but 
containing also factors (/i _1 /7) of the double value of U 2 or Ll\ or /i -1 so that k + / + m in 
the condition of summation is replaced by k + / + m + 2s l2 \ 
All equalities here and below hold modulo 0(h M ). 

12 ) More precisely 

(2.4)** U 2 ^U* = U 2 + ]T a k i pqm5 {x)U 2 k Ui'Z?Zi q ii 1 - k -i-r-*- m - 5 h s , 

k+l+p+q+m+2s>2 

(2.6) ** Z 1 ^Z* = Z 1 + ]T a' klpqmm (x)U^Ul'z^^- k -'-^- m - s h s , 

k+l+p+q+m+2s>2 

(2.7) ** xj ^xf=xj + Yl a^^U'ul'Z^Z^trX- 1 -^-^^, 

k+l+m+2s>l 
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2.1.3 Note that due to (2.5) for Uf there exists /z^-FIO T such that 

(2.14) FT = TF = I and FU*T = (hD 2 - i/ix 2 ). 
Then 

(2.15) F{ U* • U* ] )F = {h 2 D\ + ^ 2 x 2 2 ) 

and also due to (2.8), (2.9) ^ X FZ*T, Fx* J 7 and FB kpqms T are fi^h-PDOs with 

symbols depending only on x',£' (x' = (xi,X3,X4)). Let us redefine Z* re = f fi -1 F Z* J 7 , 

B k Pqm s re =^* B k Pq ms^ and A# re =FA#F, then also automatically Z* ]re = )i~ l FZ*F is an 
adjoint operator. I remind that f 2 and B™ pqms are operators with the real-valued symbols 
and A* is again defined by (2.12)*. 

Now one can decompose u into series 



(2.16) u(x,y,t)= ^ 1 ^n\x',y',t)T n {x 2 )T n \y 2 ), T n (x 2 ) d ^ v n {^ 2 h~ l l 2 x 2 ) 

n,n'& 



where v n are (real-valued and orthonormal) Hermite functions and then replace operator 
(h 2 D 2 + /i 2 xf) by (2n + l)/i/), thus reducing operator /A* to the family of 3-dimensional 
^/j-PDOs 



(2.17) A„ d = f \ (f*(2n + l)»h + (Z* . Zf) _ V #) + 



2k+p+q+2m+2s>l 

again the first line is the main part of operator and f&, V* are transformed f 2 and V. 
Remark 2.4. (i) In this and below formulae one actually needs to consider 

(2.18) n < C Q /{nh) + C\ \ogh\. 

Really, I am interested in the domain in the phase space where {a < c} or, equivalently, 
{\U 2 \ < c, \Zx\ < c}. Note that condition {\U 2 \ < g} does not contradict to the logarithmic 
uncertainty principle as long as q 2 > Cfih\ log h\ and thus one can take an upper bound for 
a as g 2 = c + Cjih\ \ogh\ which implies (2.18). This estimate perfectly suits my purposes 
as long as 

(2.19) h~ s < fi < eh' 1 ] \ogh\~ 1 . 
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(ii) From the operator point of view which one needs to apply only as 

(2.20) e/T 1 ! log/)]" 1 < n < e/T 1 
in the final analysis I need to consider only n < Cq/(^iH). 

Remark 2.5. (i) Obviously ff, V# and yU _1 Z* are /i" 1 /)-PDOs with symbols /yoV|/ 2 , VoV|/ 2 
and fi~ 1 Zi o \|/ 2 respectively where : M 6 3 (x', £') — > x G K 4 is some smooth map. 
Actually I am interested in the symbols of ff and V* only modulo symbols of fi~ 1 Z-f, 
li~ x Zf* and therefore actually I am interested only in the map T. — > R 4 . However at this 
moment I do not have a natural parametrization of Z. 

(ii) Note that ff is not necessarily xi anymore. Still since VReZi, VlmZi, VReZ2, 
VI1T1Z2 and V/i were linearly independent 13 -* this linear independence statement is true 
for VReZ*, VlmZ* and Vf x * 14 ^ as well. Therefore after /i -1 /)-FIO transformation (in 
(x', /i~ 1 /?D / )) one can achieve ff = X\. 

So, I arrive to 

Proposition 2.6. Let condition (2.1) be fulfilled. Then by means of [i' 1 h -transform, de- 
composition (2.16) and one more fi" 1 h-transform one can reduce operator to the family of 
3-dimensional operators (2.17) with n satisfying (2.18) and with 

(2.21) '> _1 { Z i # - Z i #t ) = f t = xi mod (Z x , Z\). 
2.2 Canonical form away from A 

2.2.1 Let us consider first the canonical form in the subdomain, disjoint from A which 
means that 



(2.22) \{Zx,fi}\ >e - 

Let us reduce (2.17) to a canonical form (multiplied by an elliptic operator). Note that 
the multiplication of Zf by symbol (3 implies the multiplication of Zf by ft and therefore 
the multiplication of ff by \(3\ 2 modulo Oi and finally the multiplication of {Zf, ff} by 
\f3\ 2 f3 modulo 0\ as well where here and below O m means the sum of monoms of the type 
«//'/</W / ~ m x{Z 1 * J Zf^ k with m > j + k + / and smooth coefficients a^. Also to maintain 
(2.12) and (2.17) one needs to multiply operators A and A n (and thus ff, V* and all the 
perturbation terms) by |/3| 2 . 

13 ) All the linear independence statements are uniform. 

14 ' Where in the former case it was V x ,j and in the latter one V X ',f; I remind that x' — (x1.x3.x4). 
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Then picking up 

(3 = {z*,f*y.\{z*j*}\-^ 

one can achieve 

(2.23) {Z*,f*} = 1 + a 

with a G O m and m — 1. Consequently, multiplying Zf by (1 + (3) with a G O m one 
changes {Zf , ff } by {2(3 + /J^){Zf , fi} modulo O m +i and picking /3 = -c^ — |a one can 
achieve (2.23) with a G (5^ with arbitrarily large M (and after asymptotic summation with 
M = oo). Unfortunately elements of (in contrast to elements of Om) are not necessarily 
negligible unless |xi| < h s . 

However, it is easy to see that multiplying Zf by + with (3 G O m nOoo one changes 
{Z*,f*} by {Zf ,0} with 

= UP + 0*) + //^{Zf , ~ i^{Zf\ (3} 

modulo O m r\Ooo. Obviously there exists (3 G OiflOoo eliminating error a in (2.23) modulo 
0\ Pi Ooq. One can see easily that the only restriction to G O m _i fl (5oo is that it must be 
real valued. 

Therefore step by step one can make Re a G Om H Ooo'. {Re Zi, fi} = 1 mod Om H Ooo 
with arbitrarily large M. But then setting ff = //i~ 1 {Zf, Zf *} = 2 / u~ 1 {Re Zf , Im Zf } one 

gets 

{ Re Zf, {ImZf.fi}} = {lm Z* , {Re Zf, fi}} G M -i n 

and since {Im Zf, fi} G O m _i as fi = and {ReZf , ff} ~ 1 I derive that {Im Zf , f x } G 
Om-i n Ooo. Then 

(2.24) -ifjL-^Zf, Z* ] } = Iff, {Zf, f*} = 1 mod 0{^ M ) 
and then after an appropriate /i _1 /?-FIO transformation 

(2.25) Z* = hDx + i{hD 3 - ^ixf) 

(again I am using a linear independence of V Re Zf , V Im Zf and Vfi). However, one can 
replace operators (3{x' , fi~ 1 hD') by their decompositions into powers of fi~ 1 hDi = |(Zf + 
Zf^) and fi~ 1 hD 3 — = ^/i _1 (Zf — Zf thus arriving to decomposition (2.26) below 
with all the operators of the /3{xi,x", fi~ 1 hD/\; /i" 1 /?) and the following statement is proven: 
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Proposition 2.7. Let conditions (2.1) and (2.22) be fulfilled. Then by means of n~ x h- 
transform, decomposition (2.16) and one more /i -1 h -transform one can reduce original op- 
erator A to the family of 3-dimensional operators 

(2.26) A n = \a((Z*.Z^)-W*+ 

E B ^s x ((2/7+ l)nh) k Z*<>Zp V- 2 *- 2 ™-^) 

2/c+p+g+2m+2s>3 

where symbol a = \{Z 1 , fi}| 2//3 is bounded and disjoint from 0, Zf given by (2.25) and 

(2.27) \Nf is an operator with the symbol a _1 (V — /^(2n + l)fih) o U' is £/ie smooth 
diffeomorphism R 4 — > IR 4 , wit/i transforming Span((0, 1, 0; 0), (0, 0, 0; 1)) into K 2 , 
Span ((1, 0, 0, 0), (0, 0, 1, 0)) into Kx and, in particular, (0,0, 1;0) into some element of Kx 
o/(0, *,*,*) /orm. | det OM/ 1 1 = f 2 . 

Again the first line in (2.26) is the main part. 

2.2.2 Now let us consider domain where condition (2.22) is violated. Note first that 
{(x,£) : fi = Zi = {Z\, fi} = 0} was an involutive manifold of codimension 5 in IR 4 x M 4 
with 

(2.28) V/i, VReZi, VlmZi, V{ Re Z x , h} and V{ I m Z x , linearly independent on A' 13 ). 
Also note that 

{Z x , xi} = x 3 + /x 4 , {Zx, x 3 } = x 3 + /x 4 , {Zx, x 4 } = x 4 - /x 3 as xi = 

and 

(2.29) {a*, P*} = {a, 0} + /> -1 ^ -1 ({Z 2 f , «}{Z 2 , /?} - {Z 2 , «}{Z 2 f , /?}) + 0(/^ 2 ). 
However {Z 2 ,xj} are actually arbitrary and therefore while we know that 

(2.30) A* = {f* = Z* = {Zf, f*} = 0} 
is a manifold of codimension 5 with 

(2.31) Vff, V ReZf, V Im Zf, V{ReZf, f x } and V{lm Zf, f*} linearly independent on 
A # 13),14). 
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but its symplectic structure is not fixed and there is no canonical form uniformly near 
Instead let us consider point z (= (0, 0) £ M 6 for simplicity of notations) with 

(2.32) \fi\<l = er 2 , \{Z 1 ,f 1 }\ = r with r > fi 5 ~ 1/3 , 

and reduce an operator to a canonical form in its 5(7, r)-vicinity with respect to (xi; X2, x 3 , X4), 
leaving domain > e|{Zi, h}\ 2 } for a later. 

Then after reduction of the previous subsection condition (2.32) remains valid for Z* and 
therefore now multiplication by a symbol (3 = {Zf, ff}^ ■ \{Zf , f^}! -4 / 3 /- 1 / 3 provides 
a modified equality (2.23), namely, 

(2.33) {Z*,f*} = r + a. 

However, (3 is an uniformly smooth symbol; more precisely (3 = /3(r _1 x', r -1 £') with uni- 
formly smooth symbol [3, which means that after original rescaling x' 1— > rx' 15 ) (3 = 
(3(x, n -1 hr~ 2 D') and x^ 1 a is of the same type. 

Continuing as before 16 -' one can achieve a modified equality (2.24), namely, 

(2.34) -i^r^Zf, Z* ] } = 2ff, {Zf, ff} = r mod 0(// _s ) 
where now /i _1 r~ 1 Z 1 is /i _1 /?r~ 2 -PDO and then modified equality (2.25) 

(2.35) Z* = (rh,D 1 + h*D A ) + i(h*D 3 - -fine 2 ) 

where extra terms /?*D 4 and /)*D 3 appear because (original) V Re Z\ and V Im Z\ are linearly 
independent, /?* = hr^ 1 . 

Then after another rescaling xi i— > rxi, D 1 r~ 1 D 1 (2.35) becomes 

Z* ee (/j,Di + /j*D 4 ) + /(/7,D 3 - ^fir 3 x 2 ) 
and after transformation (xi, x 3 , x 4 ; Di, D 3 , D 4 ) 1— > (xi, x 3 , X4 + Xi; Dx — D 4 , D 3 , D 4 ) becomes 

(2.36) Zf = /?*Di + /'(/?* D 3 - ^x 2 ), h* = hr\ ^ = fir 3 

which is exactly equality (2.25) with h and fi replaced by /)* and /i* respectively. 

15 ) Since symplectic structure of A# is not defined one cannot confine this non-smoothness to a couple of 
coordinates. 

16 ' With effective semiclassical parameter fi^ 1 hr^ 2 rather than /i -1 /?. 
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Meanwhile, in the original operators of type f3(r x' , jX ^hr^D'^/i^hr 2 ) (with the 
smooth symbol (3) were transformed subsequently into f3(x' , /i" 1 hr" 1 D'\ fi~ 1 hr~ 2 ), then into 
(3{rxi, x", yU _1 /7r _3 D 1 , hr~ 2 D h '; yU _1 /7r~ 2 ) and finally into 

/3(rxi, x", fi'Hr^Di, fi^hr 2 D 3 , fi^hr^D^, fi^hr 3 ). 

So, I arrive to (2.17)-like decomposition with Z[ given by (2.36) and a, \N n and B™ pqms 
of the above type. However, one can replace such operators by their decompositions into 
powers of f i- 1 hr~ 3 D 1 = \^- x r' 2 (Z* + Zf + ) and ^hr~ 2 D z - \r 2 x 2 = ^/i" 1 r'^Zf - Z# f ) 
thus arriving to decomposition (2.38) below with all the operators of the simplified type 

(2.37) /3{rx lt x", ^Hr 2 D 4 ; fi^hr 3 ); 

factors /i _1 /7r -4 appear because one needs to commute Zf, Z* between themselves and 
with operators of (2.37)-type. 
So, I finally arrive to 

Proposition 2.8. Let conditions (2.1) and (2.32) be fulfilled. Then by means of [T x h-FIO 
transform, decomposition (2.14) and the series of rescalings and yr x hr^ 2 -FIO transform 
one can reduce described above original operator A to the family of 3-dimensional operators 

(2.38) An*±a#((Z*.Z#t)-W*+ 

B^ kpqms X ((2/7 + \)lih) k Z* P Zf «y-2*-2m-p-„-s^ r -2p-2«,-4m-4s\ 

2k+p+q+2m+2s>3 

where a* is transformed operator with symbol a = \{Z\, fi}| 2 / 3 r -2 / 3 which is bounded and 
disjoint from and Zf given by (2.35). 

Again the first line in (2.38) is the main part. 

Remark 2.9. (i) Here operators cr*, Wf and B™ pqms are of the (2.37)-type with uniformly 
smooth symbol (3. Therefore these symbols are quantizable provided under condition (2.32) 
for sure; 

(ii) I leave for a section 3 the the calculation of the symbol of W and the discussion of the 
corresponding diffeomorphism. 



Figure 1: Reduction was done in the colored zone. 
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2.2.3 Let us look what happens with the cut-off symbol ip as a result of these transfor- 
mations. After transformations of subsubsections 2.1.1, 2.1.2 after decomposition (2.16) it 
becomes a matrix operator 

k,l 

with basis-shifting (in /- 2 (R)) operators (2(n + l)^/?) -1 / 2 ^ and (2n/x/7) _1 / 2 L/ 2 - 

Then after transformations of subsubsection 2.2.2 these matrix operator becomes 

(2.39) <Pu pqm ^^'zf p Zt^ q ^~ k ~'~^ q ~ S ^r- 2 ^^ 2n ^ 

k,l ,p,q,m,s 

with the scalar main part of the form (2.37). 

2.3 Canonical form in the strictly outer zone 
2.3.1 Now one can properly redefine the first part of the outer zone 

(2.40) Z outJ = | Q/- 1/2 dist(x, A) 1/2 < dist(x, Z) < edist(x, A) 2 | 

which in the notations dist(x, T.) = |xi| and dist(x, A) x r is {yii" 1 / 2 /' 1 / 2 < |xi| < er 2 }. 
Note that after rescaling of the previous section xi i— > xir~ 2 , /i \— > = fir 3 it becomes 
{|*i | > C[u 1/2 }. 

In this subsubsection I restrict myself to the part of strictly outer zone 

(2.41) Z* ut l = | Qi 3 ^ 1/2 dist(x, A) 1/2 < dist(x, Z) < edist(x, A) 2 } 

which in the same notations is | / u~ 1 / 2 r 1 / 2 < |xi| < er 2 } and after becomes {|xi| > Cfu S 1 ^ 2 } 
(may be with a different 5 > 0) and provides a nice quantization of all symbols below. 

Then since /i* > h~ s one can apply the canonical form of the Schrodinger operator 
with the strong magnetic field. Let us consider an operator obtained after reduction of the 
previous subsection; there \xi\ x 7* = / -/r~ 2 if originally |xi| x 7. Notation 7 = r 2 was 
temporary. 

Applying rescaling (x!,x 3 ) 1— > (x^7* < -*37*) one arrives to transformed (2.36) of the same 
type but with /)* and /1* replaced by /?*7~ 1 = hr^~ x and /i*7 2 = /i7 2 r~ 1 respectively: 

(2.42) Z* = n-fr^KDi + i(hD 3 + ^x 2 )) 
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with |xi| x 1 and 

(2.43) h = fi-Hr 2 ^ 3 < c^h^r 1 ' 2 < ch 1 ' 2 . 

Also n- l hr~ 2 D z h+ ^h-y^D^. 

Let us consider the part where xi > 0; another part is separated and treated exactly in 
the same manner. Introducing 

(2.44) (A = x~ 1/2 (hDt + i(hD 3 + ^x 2 )) 
one can see easily that 

(2.45) |*A| <ci/ -1 , i/ = /17V 1 > h- 5 , 

(2.46) {U lt Ul} = l mod L/ 1( L/f 

where (2.45) follows from |Z*| < c and (2.46) is understood in the sense of ^-symbols. 

Therefore there exists h-FlO transformation T' of (x 1 , x 3 , HDi, hD 3 ) thus not affecting 
(X4, Li~ 1 hr~ 2 D/\) so that 

(2.47) (A i-> if = hD l + Dfi mod 2 (L/i, U\)\ 
then operator (2.38) is transformed into 

(2.48) -aa x [u 2 U* • Up - 1/1/**+ 

£ B^x((2#*l),,/,)>i#)>^ 

2/c+p+q+2m+2s+2/>3 

with 1/1/** the result of transformation of xf 1 !/!/* . 

Here all the operators <x, a%, 1/1/** and B™ pqmsj of the type 

P(x lt x 3 , x 4 ; ftDi, 7iD 3 , yT x hr 2 D A ). 

Note that 7m/ = hr^" 1 <C /i 1 / 2 /?r 1 / 2 < /7 1//2 in the zone in question. One can decompose all 
such operators into power series with respect to x 1; fiD^ and rewrite (2.48) with all operators 
of the type 

(2.49) /3(x3,x 4 ;nD3,/i- 1 /?r- 2 D 4 ), h = fi^hr 2 ^ 3 
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and with (z/i/*) P , {yUf^ q replaced by [uUf) P+P v p ', (vUf^ q+q v q ' respectively since 
v > h' 5 . 

Repeating arguments of the previous subsections one can transform this operator into 

(2.50) -o-o-i (y 2 Ut • - W**+ 

E B41 x (^ 2 Uf • Uf) q u 2 - 2 ^'h'+ 

2<j+2/+2/>3 

E B™ kqjmsl x ((2n + l)fih) k (u 2 U* . Uf ] ) q ^- 2k - 2m - 2q - 5 h 5 r- Aq - Am - As u- 2 ^) 

2fc+2q+2m+2s+2/>3 

where the second line comes from decomposing of 1/1/** and because it is rescaled before 

(2.51) 1^,1 < C 1 r\ 

Now one can apply decomposition with respect to h~ l ' 4 v p (h~ 1 ' 2 Xi): 

(2.52) u nn ,(x',y',t)= ^ 1 »n P n> P >(x",y",t)v p (h~ 1 / 2 x 1 )v p ,(h- 1 / 2 y 1 ) 

p,p'ez+ 

with x" = (x 3 ,x 4 ). Then the operator (2.50) becomes 

(2.53) l -aa x {[2p + l)/i/> 7 - W#*+ 

E B% x ((2p + l)^h 1 ) q v 2 - 2 ^h'+ 

2<7+2/+2j>3 

E B ^s, X ((2/7+ l)A*/l)*((2p + l)p/7 7 ) V 2 - 2 ^ 2 ^^^- 4 ^^ 2 ^'). 

2k+2(7+2m+2s+2/>3 

since i/ 2 /i = p/77 and 2k + 2q + 2m + 2s + 2j + 21 > 3 is equivalent to /c + qr + /77 + s+_/ + / > 2. 

In the decomposition (2.53) the role of m, s, _/ and / is just to bound a magnitude 
and indicate dependence on p and h since dependence on r, 7 (but not magnitude) is not 
important. However v = \x ■ 7 2 r~ 1 < fir 2 , h = p -1 /? ■ r 2 7 -3 > p _1 /jr -4 ; therefore one can 
join terms with s > or m > to the terms with s = and m = without changing s + /, 
m + j thus getting the simplified correction term 

E fi 2U/ x (( 2 " + l)/^)*((2p + l)^7)V 2 -^^ s r-\- 2 ^' 

/(+<7+j+/>2 

Therefore I arrive to the following 
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Proposition 2.10. Let condition (2.1) be fulfilled. Let us consider zone Z* utl , described 
by (240). 

Then by means of fi~ 1 hr~ 2 -FIO transforms, decomposition (2.16), the "special" h-FIO 
transform and decomposition (2.52) one can reduce the original operator A to the family of 
2-dimensional operators 

(2.54) A„ p = ~(T(Ji{H pn + 

E ((2n + l) / i/ ) )^((2p + l) / i/, 7 )V 2 " 2 ^ 2 ^ s ^ 4 V- 2 ^') 

k+q+j+l>2 

with n < Co/(/i/)), p < Co/(filr-f) where all operators are of (2. 4 9) -type and a, o\ are 
bounded and disjoint from 0. 

Remark 2.11. (i) While in the canonical form (2.38) only one cyclotron movement was 
separated, in the canonical form (2.54) both of them are separated leaving only pure drift 
movement which will be studied in the next section; 

(ii) I leave for a section 3 the discussion of the actual elements where reduction was done, 
corresponding diffeomorphism and the calculation of the symbol of 1/1/; 

(hi) Cut-off symbol ip which after previous transformations became (2.39) now becomes 

(2.55) J] <f> kipqsl X L // L /t/c^#p^#t^-/-^2 m -p- q -s^s r -2p-2 q -4 m -4s zy -2i-/-^/ 
k,i,p,q,s,l 

with (2.49)-type operators (f>ki pq si- 

2.3.2 Now my goal is to achieve a similar reduction in the second part of the outer zone 

(2.56) Z* out ll = |Cmax(dist(x,A) 2 ,/ <5 - 2/3 ) < dist(x, E) < e}. 

Then the previous subsection 2.2 becomes irrelevant; however I can start from (2.17). 
Again with no loss of the generality one can assume that = x-y. Let us consider a point 
z = (x , £ ) G M 6 and an element (2.55) with 7 = |x{| and r = 7 1 / 2 (exactly as in the 
previous subsection it was defined the other way. Without any loss of the generality one 
can assume that z = (7, 0, ... , 0). 

Consider Taylor decomposition (with a remainder) of /i _1 Z* with respect to x', fi~ 1 hD' 
and r (in r-vicinity of all of them are 0(r)); then due to {Zf , Z^} = O(r), {Zi, xi} = 
0(r) with no loss of the generality one can assume that yU _1 Z* = /i _1 /?D 3 + ifi~ 1 hD 4 + 
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0(r 2 ). Note that p^/iDi could be either as square (or higher degree) or with a linear 
or quadratic factor with respect to other variables. Also note that the quadratic terms in 
0(r 2 ) containing exclusively x 3 , x 4 must commute with D 3 + /D 4 (otherwise {Zf, Zf^} ~ 2xi 
would not be possible) and therefore must be of the form c(x 3 + /x 4 ) 2 . However one could 
remove such terms by /i^/j-PDO affecting only x // , / u~ 1 /?D" and leaving x 1 , fi~ 1 hD 1 intact. 
So, terms in the decomposition must contain either one of the factors x 1; fi~ 1 hD$, fi^^-hD^, 
or two factors /i~ 1 /)D 1 , or three factors x 3 ,x 4 ; also x 1 must be accompanied by 0(r)-type 
factor. 

Then scaling x" \— > rx', xi > r 2 xi, D" \— > r~ 1 D", D\ i— > r~ 2 D\ and dividing by r 3 
one gets completely legitimate p _1 /)r~ 4 -PDO where each factor /i~ 1 /7D" can accommodate 
division by r 3 and /i~ 1 /7D 1 can accommodate division by r 2 but there either at least two 
such factors, or there is an extra O(r) type factor. This operator, with the semiclassical 
parameter h = /i^/ry" 2 satisfies again {Zf , Z^} ~ 2x 1 but must be considered near point 
(1, 0, ... , 0) and then by the standard arguments one can reduce operator to the canonical 
form (2.54) with all the operators of the type 

(2.57) (3{x 3 ,x A ;hD 3 ,hD 4 ), h = fi~ 1 h-f~ 2 . 

Note that as r = 7 1 / 2 this is consistent with the definition of h in the zone Z out j but form 

(2.57) is more general then (2.49). 
So, I arrive to 

Proposition 2.12. Let condition (2.1) be fulfilled. Then one can reduce the original op- 
erator A to the family (2.54) °f 2- dimensional operators with n < Co/(fih), p < Co/(^h'~f) 
where all operators are of (2. 49) -type and a, a 1 are bounded and disjoint from 0. 

Remark 2.13. (i) All statements (i)-(iii) of remark 2.11 remain valid. 

2.3.3 I would like to finish this section by the following 
Remark 2.14. In 2D-case canonical form would be [Ivr6] 

(2.58) A p = aa 1 ((2p + 1)/i/j 7 - W*#+ 

£ Bl x ((2p + l)^ 7 ) ' x p 2 - 2 '- 2 ^/7V 4i - 3t ) 

2/+2/+2t>3 

with 

(2.59) W## = P(T l *i, // _1 fry- 2 P3). 

There instead of 2-parameter family of 2-dimensional PDOs I had only 1-parametric family 
of 1-dimensional PDOs. 
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3 Estimates 

In this section I derive remainder estimates with the main part in the standard but rather 
implicit form which is the sum of expressions (0.12) 

(3-1) h~ l J2 f (F t ^rXTAt)YuQl y ) dr 

J — oo 

with Q L partition of unity and an appropriate T t where as in my previous papers Xr(t) = 
x{t/T) and x is supported in [—1, 1] and equal 1 at [— ~, \], while the remainder estimate 
could depend on non- degeneracy condition. 

I am going to consider different zones and apply different approaches and use different 
canonical forms in each of them; I start from the strictly outer zone. 



7 




Figure 2: My first target is strictly outer zone. 



3.1 Estimates in the strictly outer zone. I 

I will start from the most massive Strictly Outer Zone I Z* ut , described in (2.40) with an 
arbitrarily small exponent 5q > and sufficiently small constant e > 0. 

3.1.1 After rescalings and transformations described in the previous section the original 
operator somewhere there was reduced to the family of two-dimensional (^-3, ^-PDOs 
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A pn {x", h 3 D 3 , h 4 D 4 ) given by (2.54) with 

(3.2) h 3 = fi^hr 2 ^ 3 < i?/ 2 h, h A = fi^hr 2 < ^ 1/3 h; 
one needs to consider only magnetic numbers 

(3.3) n < Coi^h)- 1 , p < C^/x/ry) -1 as dist(x, E) x 7; 

these two restrictions mean that energy levels do not exceed e. 
These final operators are defined in 6(0, 1). Therefore 

(3.4) A contribution of each (pair, element) 17 ) to to the asymptotics is 0(^3 = 
0(/i 2 /j- 2 7 3 ) 18) - 

However, one can see easily that this element before the final reduction was the "slab" 
6(0, 1) fl {|xi — xi| + |x 3 — x 3 | < 57*} with 7* = 7r~ 2 , |xi| x 7* where here 6(0, 1) is an 
image of the original element 6(y', r) n Z out j. Therefore the total contribution of the given 
pair (p, n) over subzone 

Z( rn) ={7 < dist(x, E) < 27, r < dist(x, A) < 2r} 

with r > e7 1 / 2 to the asymptotics is 0(/i 2 /?~ 2 7 3 x r 2 7 _1 ) = 0(yU 2 /)~ 2 7 2 r 2 ) and the total 
contribution of subzone 2(r, 7 ) to the asymptotics is 0(/i 2 /? _2 7 2 r 2 x (/z/?7) _1 x (/i/j) -1 ) = 
0(/?~ 4 7/" 2 ) where the second and the third factors are numbers of permitted indices p and 
n respectively. 

Then the total contribution of zone Z* ut i to asymptotics is 0(/?~ 4 ) 19 ^. 

The same approach works in Z* ut n : here one has the family of two-dimensional /z 3 -PDO 
A P n(x", h 3 D 3 , h 3 Da) with h 3 = /i" 1 /^" 2 . Therefore a contribution of each (pair, element) to 
to the asymptotics is 0(h^ 2 ) = 0(yU 2 /?~ 2 7 4 ). However there was nondiscriminatory rescaling 
(x 3 , x 4 , £ 3 , £ 4 ) 1 — > (r _1 x 3 , r _1 x 4 , r _1 ^ 3 , r^ 1 ^) while the number of such balls should be 0(r~ 2 ) 
since codim^ A = 2. 

Therefore the total contribution of the given pair (p, n) over subzone 

Z [rn) ={7 < dist(x, E) < 27, dist(x, A) < 2r} 

with r = 7 1 / 2 to the asymptotics is 0( / u 2 /7 _2 7 3 ), the total contribution of subzone ^( r , 7 ) 
to the asymptotics is 0(h~ A j 2 ) and the total contribution of zone Z* ut t to asymptotics is 
0(/r 4 ). ' 

Now I need to analyze the contribution of such elements to the remainder estimate. 

17 > Where here and below a pair is (p, n) and an element means the final element 6(0, 1). 

18 ) And for many pairs it is actually of this magnitude. 

19 ) And is actually of this magnitude under condition (0.2) as [i < eft -1 . These results are consistent with 
what one can get without reduction. 
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3.1.2 Let us consider some index pair (p, n) and some final 6(0, 1) element. Assume first 
that after all these transformations operator A pn is regular on this (pair, element) in the 
sense of (3. 5), (3. 6) and (3.8) 

(3.5) Va\^x"rAn\ < C J Va : < \a\ < K, 
with rj a = 1 as a = (*, 0, *, 0) and rj a = 'yr~ 2 otherwise and I satisfying 

(3.6) £>h- Ko . 

Then the propagation speed with respect to (xs,hsD 3 ) is 0{lh-i/ti) = O^fi' 1 r 2 7~ 3 ) and 
the propagation speed with respect to (x 4 , /l 4 D 4 ) is 0{lh^jh x r 2 /^) = 0(£h 3 /h) as well. 
Therefore the dynamics starting in 6(0, |) is retained in 6(0, 1) as \ t\ < 7~i with 

(3.7) T x = ehr 1 ^ 1 x fi^r 2 ?- 1 ; 

I remind that ^> h 4 . In the same time under respective assumptions 
(3-8) | V^Apnl + (\V x ^A pn \ x I C = ir' 2 

one can pick up To from the uncertainty principle 20 ) ihjh^To > hjh~ 5 ' 21 ) or equivalently 

(3.9) T = Ch x - S 'r x . 

with an arbitrarily small exponent 5' > 0. 

Then the contribution of this (pair, element) to the remainder estimate does not exceed 

(3.10) CfcT 1 /^ 1 x T Tf 1 x h^h^' = /i/r^'r 2 

where Ch^h^ 1 x T estimates sup| T | <e \F t ^ h -i T XT{t)f~ (uip) | and factor Tf 1 is due to the 
standard Tauberian approach. 

So, in comparison with the estimate of the contribution to the asymptotics (3.4) the 
estimate of the contribution to the asymptotics picked an extra factor h^h^ 5 with an 
arbitrarily small exponent 5' > 0; actually one can reduce 5' to but it is not needed 
here. 

Let us pick up (p, element), then break this element into C-subelements with respect to 
(x 4 , £ 4 ); then on each such subelement operator A pn is elliptic 22 ) \A pn \ > e£ for all indices n 

20%> I don't need the logarithmic uncertainty principle in this place and use the simpler version of the 
microlocal uncertainty principle. 
21 ) Alternatively one could scale x 4 i— > (ft 4 /fi,3) 1 / 2 x 4 , £ 4 i— > (fru/fr?,) 1 ! 2 ^!, ft 4 i— > S3. 

22 - ) See next subsubscction for justification; to prove (3.14) and (3.15) which are main results of this 
subsubsection one does not need these arguments. 
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but C(£/(/i/?) + l). Therefore the total contribution to the remainder the pair (p, element) 
combination such that conditions (3. 5), (3. 6) and (3.8) are fulfilled with t < I does not 
exceed 

(3.11) Qz/T^V x (iifih)- 1 + l) = C (£ + fth) h~ 2 - s ' r 2 . 
Then summation over Z( r i7 ) results in 

(3.12) C/rV 2 x fi-'hr 2 ^ 3 x h~ s ' {£ + /i/?) = C(I + ^h)^ 1 h~ 3 ~ 5 ' r 4 ^ 2 
and therefore 

(3.13) The total contribution to the remainder of all (pair, element) combinations residing 
in Z( rn \ C Z* ut i and satisfying conditions (3. 5), (3. 6) and (3.8) with £ < £ does not exceed 
(3.12)'. 

It implies 

(3.14) The total contribution to the remainder of all (pair, element) combinations residing 
in the far outer zone Z* ut i R {dist(x, T.) > /i^ 1 / 4 /)^ dist(x, A) 2 } and satisfying conditions 
(3. 5), (3. 6) and (3.8) does not exceed Q/" 1 / 2 /? -3 . 

The same arguments work for Zu,i) with r = 7 1 / 2 resulting in the estimate C /x" 1 h~ 3 ~ 8 ' of 
its contribution to the remainder; summation over Z* ut n results in expression 
Therefore 

(3.15) The total contribution to the remainder of all (pair, element) combinations residing 
in Z* ut u and satisfying conditions (3. 5), (3. 6) and (3.8) with ( = 1 23 ^ does not exceed 

I consider the exceptional regular elements (not covered by (3.14) in the next subsubsection. 
The rest of the analysis of this subsection is devoted to detecting and analysis of the irregular 
elements. 



' One should pick up here £1 = 1 because hj, — h^. 
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3.1.3 Note that 

(3.16) H pn = ^(V -(2n+l)fihf 2 -(2p + l) f ihf 1 )^ o V 

with uniformly smooth map W which does not depend on V; this exonerates elliptic argu- 
ments of the previous subsubsection and in particular (3.13) since perturbation part satisfies 
\d n B pn \ <C fih; one can consider n and p as continuous parameters. 

Proposition 3.1. (i) After reduction inZ* utl symbols Woty, f\o^ ', f 2 oty satisfy conditions 
(3.5), (3.6) with £ = ^r~ l and indicated above rj a and (; furthermore 

(3.17) |V X3 , 6 /ioV|/|x 7 ; 

(ii) After reduction in Z* utll symbols W o v|/ ; fi o v|/ ; f 2 o v|/ satisfy conditions (3.5), (3.6) 
with I = 7/" ; furthermore 

(3.18) \V x »,z»fi o v|/| x 7. 

Proof. Proof follows easily from the reductions. More precise results will be needed and 
proven later. □ 

Therefore 

(3.19) As p > C /(fihr) symbols H pq satisfy conditions (3. 5), (3. 6) and (3.8 with £ = fitryp > 
Ql/ r. 

On the other hand one can prove easily from the reduction that 

(3.20) Perturbation symbols B pn satisfy conditions (3.5) and (3.6) with I replaced by £' = 
/>VV(p + l) 2 + ^h^r 2 < /i/7 7 (p + 1). 

Note that if I restrict myself only by p < Co/(/i/?r) (condition to be explained later) 
instead of p < Cof(fih'y), transition from (3.11) to (3.12) picks up an extra factor 7/r and 
therefore 

(3.21) The total contribution to the remainder of all (pair, element) combinations with p < 
Co/(fihr), residing in Z* ut , and satisfying conditions (3. 5), (3. 6) and (3.8) does not exceed 
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Now I want to cover indices p > Co/(^hr) and elements in Z* ut , which are not in 
the far outer zone and prove that their total contribution to the remainder also does not 
exceed Cp -1 / 2 /? -3 . To achieve this I plan gain a factor 7r~ 2 in the estimates of the previous 
subsubsection, starting from (3.10) by increasing T\ given by (3.7) to 

(3.22) 7 2 = er/i, 7 2 = e/i 7 2 . 

Note that (3.7) appeared because it is a time for which propagation is confined to 6(0, 1) 
in the final reduction, for larger t shift with respect to (x 3 ,£ 3 ) would be too large while 
shift with respect to (x 4 , £ 4 ) would be still less than e as \ t\ < T 2 . 

However considering propagation after intermediate reduction (i.e. reduction of subsec- 
tion 2.2) the shift with respect to x 3 would be of magnitude //7 1 7~ 2 (/i/)7p)| t\ which would 
be observable but less than 1 as 7~i < \t\ < 7~ 2 ; this propagation just moves to the adjacent 
"slabs" in the intermediate element 6(0, 1) corresponding to the different final elements. 

Then the contribution of such (pair, element) to the remainder estimate is given by a 
modified expression (3.10) 

(3.23) Ch^ 1 ^ 1 x ToT"^ 1 x p/r^V -1 x /i/r 1 -^/^)- 1 
which leads to modified expression (3.12) 

(3.24) p/T^V" 2 x (fih)" 1 x (/i/ry) -1 ] \ogh\ x r 2 7 " ! 

where the second, the third and the fourth factors estimate respectively the number of 
indices n, the sum of (p/?7p) _1 with respect to p and the number of the final elements with 
given (r, 7 ). 

Then the sum with respect to 7 could be estimated by by the same expression with 
7 = /i -1 / 2 ^/" 1 / 2 , which is C/i~ 1 / 2-,5 /j~ 3-5 '/' 3 / 2 | log h\; finally summation with respect to r 
results in 0(/i _1 / 2 /)~ 3 ) as 5' > is small enough.. 

Therefore I arrive to 

(3.25) The total contribution to the remainder of all (pair, element) combinations with p > 
Co/(fxhr), residing in Z* ut t does not exceed Cp" 1 / 2 /?" 3 . 

which together with (3.21 concludes the analysis of regular elements; their total contribution 
to the remainder is 0(p~ 1 / 2 /)^ 3 ). 
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3.1.4 Therefore in what follows I am interested only in non- regular elements, with p < 
Co/(/i/?r); then H pn with n violating ellipticity condition satisfies (3. 5), (3. 6) with £ = •yr~ 1 

in z Lt,i and with ^ = 7 r_1 > Va = (1, 1, 1, 1) in 2 * ut „. 

First consider elements in Z* ut t with fixed p < Co/(fJ>hr) and n. After rescaling 

(3.26) x 4 i-> 7r~ 2 x 4 , £ 4 i-> 7/"~ 2 £ 4 ft 4 i— > ft 4 d = /x _1 /?7 _2 r 2 

symbol 7~ 1 r/-/ pf , becomes uniformly smooth and therefore one can introduce a scaling func- 
tion 

(3.27) £ = £ pn (x", £") = er 7 - 1 | V x » e »^ pn | + 4, 4 = C^ /2 /T 5 "; 

obviously | Wl\ < \. Let us introduce ^-admissible partition of 6(0, 1). Then the symbols 
are still quantizable after rescaling (x", £") i— > £~ 1 (x" , i— >• £~ 2 hz, ft! A i— > £~ 2 ft! A . 

Further, as £ > 2£ , the propagation speed does not exceed 

(3.28) v x f j,- 1 >y- 3 r 2 \VAp„\ X ^ n ~ 2 £ 

in the corresponding partition sub element 24 ^ (in the coordinates before rescaling of the 
subelements) and therefore the dynamics started in B(z, £(z)) remains in B(z, 2£(z)) as 
\t\ < Ti with 

(3.29) 71 = e£v^ = efir 1 ^ 2 . 

On the other hand, the shift v\t\ is observable as \t\ > T with 7" defined from the 
uncertainty principle 

(3.30) vT x i > h z h~ 5 ' . 

Really, if | V X3i £ 3 .4. pn | x £7 then the propagation speed with respect to (x 3 ,£ 3 ) is exactly 
of magnitude v and the uncertainty principle is due to (3.30). Otherwise |V X4 ^ 4 ^4 P n| x £7 
and the propagation speed with respect to (x 4 , £ 4 ) is exactly of magnitude 1/7 and the 
uncertainty principle is 7 vT x £ > h 2 h~ 5 'r- 1 which is due to (3.30) as well. 
So one can pick up 

(3.31) t = n^'v^r 1 x h 1 - 5 ' 1 - 1 rr 2 . 

Therefore the contribution of this partition subelement to the remainder does not exceed 

(3.32) Cf^%-H 4 x 7~o Tf 1 = Q^-i-Vr-^ 2 

24 ) I call this partition elements "subelements" because actually 6(0, 1) itself is a rescaled partition element 
by itself. 
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while in comparison its contribution to the main part of asymptotics does not exceed 

On the other hand, the contributions of the subelement with £ < £q to both the main 
part and the remainder do not exceed 

(3.33) Ch^h'fH 4 = Cfih-^r^Il 

which is less than (3.32) (with a larger but still arbitrarily small exponent 5'). 

Let us notice that for given index p and ^-element operator A pn is elliptic as \n— n\fih > 
CqE 2 ^/^ 1 for some n; therefore ellipticity fails for no more than (Co£ 2 'yr~ 1 (fih)~ 1 + l) indices 
n. 

Then the contribution to the remainder of all the indices n (while index p and t- 
subelement remain fixed) does not exceed 

(3.34) Qi/T 1 - V^ 2 x (Qfiifihr)- 1 + l) . 

As 

(3.35) t^l^e^h 1 ' 2 ^- 1 ' 2 

this expression is x Ch~ 2 ~ 5 ' 7 3 r~ l £ A and the sum over ^-subelements in 6(0, 1) with £ sat- 
isfying (3.35) results in Ch~ 2 ~ 5 7 3 r~ 1 . 

Then the sum over all such subelements residing in >2( r , 7 ) does not exceed 

(3.36) C/T^Vr- 1 x r 2 7 " 2 x r 2 7 " 1 = OT^V 

where the second factor in the left-hand expression are the number of elements which 
appeared after rescaling (3.26) and the corresponding partition I made in the beginning of 
subsubsection and the third factor is the number of the "slabs" corresponding to the final 
elements after the second transformation in the "intermediate" ball 6(0, 1) obtained after 
the first transformation. 

Then the sum over the whole Z* utj does not exceed Ch~ 2 ~ s \ \ogh\. So far index p 
remains fixed; the sum over p < Co/(p/j) does not exceed C [i" 1 h~ 3 ~ 5 ' \ \ogh\ = 0(fi~ 1 ^ 2 h~ 3 . 
Thus I arrive to 

(3.37) The total contribution to the remainder of all (pair, subelement) combinations in zone 
Z* ut j with £ > £1 does not exceed Cfi^^h^ 3 . 
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Therefore only subelements with £ < l\ remain to be considered. 

For such subelements expression (3.34) becomes Cp/?" 1 " 5 7 2 ^ 2 ; however since without 
non- degeneracy condition the number of such subelements in 6(0, 1) could be x £~ A , the 
sum would become singular. To overcome this obstacle let us redefine the scaling function 
(3.27) in the manner more suitable for the analysis of elements with i < i\. 

1 /2 

(3.38) t=t{x'\i") = 6xn\n( 1 - 2 r 2 \V x , l ^A pn \ 2 + H~ 1 \A pn \) +£ Q , 

with 4 = ?h /2 h~ 5 ' = iT^h^'r^r. 

Note that as the minimum in the right-hand expression is achieved for (p, n) = (p, n)(x" , £"), 
for every index p inequality I > £ p = £* + e|p — p\jihr holds as the n = n(p, (x",£")) is 
selected to break ellipticity condition; this index n is unique. 

Then the contribution of all subelements with £* x A in 6(0, 1) (obtained after rescaling 
(3.26)) to the remainder (after summation with respect to n and p = p) does not exceed 
the sum with respect to p of expressions (3.34): 

(3.39) Cp/r^V J2i x + e \p - pI^)" 2 ; 

Pt^P 

I included in this sum only p with |p— p| > 1 leaving special ((n, p),subelement) combinations 
for a special consideration. Expression (3.39) obviously does not exceed Cfi~ 1 h~ 3 ~ 5 7 2 r~ 2 
and therefore the contribution of all such (pair.subelement) combinations residing in -Z( r , 7 ) 
does not exceed 

Cfi-'h'^'r 2 x r 2 7 " 2 x r 2 7 ~ 1 = C fi' 1 /r 3 - 5 ' r 2 ^ 1 

with the same origin of the second and the third factors in the left-hand expression as in 
(3.36). 

Finally after summation over Z* ut s one gets 0(p- 1/2 /7~ 3 ). Therefore I arrive to 

(3.40) The total contribution to the remainder of all ((p, n),subelement) combinations in 
zone Z* ut j with £ <£\ and (p, n) (p, n) does not exceed Cp" 1 / 2 /?" 3 . 

3.1.5 So, in what follows I need to consider only special ((n, p),subelement) combinations 
with £* < £\. Exactly for these combinations non- degeneracy condition becomes crucial. 
In the general case however A pn with (p, n) = (p, n) could be "flat" and then there would 
be no difference between estimates of the contribution of such subelement to the remainder 
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and to the main part of the asymptotics which would be Ch 3 1 ft! 4 1 x p 2 /j 2 y>r 2 . Then the 



with the same origin of the second and the third factors in the left-hand expression as in 



Finally after summation over Z* ut t one gets 0(p l ^ 2 h 3 ). Therefore I arrive to 

(3.42) The total contribution to the remainder of all ((p, n),subelement) combinations in 
zone Z* ut i with £ < l\ and (p, n) = (p, ft) does not exceed Cp 2 /? -2 . 

Combining with (3. 37), (3. 40) I arrive to 

Proposition 3.2. Under condition (2.1) the total contribution to the remainder of the zone 
Z* ut i does not exceed Cjji 2 h~ 2 + Cp" 1 / 2 /?" 3 . 

Remark 3.3. All the arguments leading to proposition 3.2 are applicable in Z* ut u with the 
following minor modifications and simplifications: 

(i) operators A pn are ^3-PDOs from the very beginning (so there is no need in (3.26) 
rescaling; there are no "slabs" etc, everything is homogeneous; 

(ii) A factor r 2 in estimates translates into 7 which makes life easier; in particular there 
is no need in the special propagation analysis leading to (3.25). 

So I arrive to 

(i) 0(p -1 / 2 /7 -3 ) estimate for the contribution into the remainder of all ((p, n),subelement) 
combinations in zone Z* ut u with either I > £ 1 or £ < £ 1 and (p, n) 7^ (p, n) and 

(ii) 0(p 2 /?~ 2 7 3 ) estimate for the contribution into the remainder of all ((p, n),subelement) 
combinations in zone Z* ut n D {dist(x, Y.) < 7} with £ <£\ and (p, n) = (p, h). 

In particular: 

Proposition 3.4. Under condition (2.1) the total contribution to the remainder of the zone 
Z* ut a does not exceed Cjji 2 h~ 2 + Cp -1 / 2 /?" 3 . 




(3.36). 
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3.1.6 Estimate (0.3) 0(/j?h~ 2 + /i" 1 / 2 /?" 3 ) is the best possible in the general case (see 
Appendix A. 3) and it coincides with the best possible O^ji^^h^) as /i < h~ 2 ^ 5 ; so \i > 
h~ 2 / 5 until the end of this subsection. 

However, one can do better under some non-degeneracy conditions. Namely, as before 
one can estimate the contribution of the special pair ((p, n),subelement) to the remainder by 
Cfih-^'^f as£>£ (see (3.32)) and by C/j, 2 r~ 2 y 5 (see (3.33)) where now £ = t defined 
by (3.38). Therefore the contribution to the remainder of all the special subelements residing 
in the final element 6(0, 1) (obtained after rescaling/partition (3.26)) does not exceed 

(3.43) C/i/T^Vr- 2 / t~ 2 dx"dC + C^-VV / dx"d£" 

J{h<e<h} J{t<£ } 

where I remind £ = ^h 1 / 2 ' 6 ' V 3 / 2 r. l x = (j, 1 / 2 h 1 ' 2 rV 2 y-V 2 . 

Then the total contribution to the remainder of all the special subelements with U'-image 
residing in -Z( r , 7 ) C Z* ut , to the remainder does not exceed 



Cu/r^V'" 2 / - - t~ 2 \6etDV\- 1 dx 

?o < t o M/- 1 < 4} n Z {rn) 



Cu, 2 h- 2 r*j 5 / _ IdetDvl/l-^x; 



Since 

(3.44) | det OM/ 1 x 7 V 4 
one can rewrite this expression as 

(3.45) C/i/T 1 - 5 ' / 1 ~ 2 r 2 t- 2 dx + 

J{£ <£*o V" 1 < £,} n Z 

Cfi 2 h~ 2 [ _ ydx 

J{£* o vi/- 1 <£ }nz 

With Z = Z( r ,-y)- 

I leave to the reader the similar analysis in Z* ut n leading to the same estimate for the 
contribution of Z^.-y) C Z* ut u with r = 7 1 / 2 : 

Proposition 3.5. Under condition (2.1) 

(i) The total contribution of the special subelements with -image residing in Z( rn ) C Z* utl 
does not exceed (3.45) with Z = 2( r , 7 ); 
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(ii) The total contribution of the special subelements with '-image residing in 2( r , 7 ) C Z* ut u 
does not exceed (3.45) with Z = Z( rn y, where r = 7 1 / 2 . 

(Hi) The total contribution of all the special subelements (with -image residing in Z* ut 
does not exceed (3.45) with Z = Z* ut , 

Note that all elements of DU' do not exceed C'jr -1 . Combining with (3.45) and (3.16) 
one can conclude that all elements of DU' -1 do not exceed C7~ 1 r. Therefore if £* was 
defined by (3.38) with H pn instead of A pn then the following inequality would hold 

(3.46) for 1 > L= e\V x (V/f 2 )\ + min \V r (V/f 2 ) - (2n + l)^h\. 

n 

Then under condition (0.10) fxl even if i* are defined by A pn and the first term in (3.44) 
does not exceed Cfih^ 1 " 6 7~ 1 r 4 while the second term vanishes and the total contribution 
of Z* ut ! to the remainder does not exceed C/i 3 / 2 /?" 1 . 
Therefore I arrive to 

Proposition 3.6. Under conditions (2.1) and (0.10) the total contribution of Z* ut to the 



remainder does not exceed Cfx x l 2 h 



-3 



Furthermore, under assumption (3.46) and condition (0.8)^ one can see easily that (3.45) 
with Z = Z* ut does not exceed C/i 3//2 /? _1 plus 



(3.47)„ C/i 2 /7~ 2 / j£ q dx 

which is 



J q 1 ~ x d 1 ■ r'dr, J q = h~ 2& " fh^^hr^-y 2 ^ 2 



with integral taken over Q = j//'^ 1 / 2 / -1 / 2 < 7 < r 2 < 1}. 

Obviously (3.47)2 does not exceed J 2 calculated as 7 = 11 s _1 / 2 , r = 1 which is 0(yU _1 / 2 /7~ 3 . 
Further, (3.47)i does not exceed J\ calculated as 7 = r = 1 which is C/j?/ 2 h~ 3 / 2 ~ 5 " with 
arbitrarily small 8" > 0. 

So, one arrives to the following statement (still under (3.46) 

Proposition 3.7. Let condition (2.1) be fulfilled. Then 

(i) Under condition (0.9) 2 the total contribution to the remainder of the zone Z* ut with does 
not exceed C [x~ x l 2 h~ 3 ; 

(ii) Under condition (0.9)\ the total contribution to the remainder estimate the zone Z* ut 
does not exceed C /j 3 ^ 2 /j~ 3 / 2 ~ 5 + C/i -1 / 2 /? -3 with arbitrarily small 5" > 0. 
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Proof. To get rid of assumption (3.26) one needs to take in account perturbation term B 
in A pn ; then £* o v|/ _1 > L — \V^B pn o V|/ 1 j . Thus it would be enough to replace £ by 
io + |V r S pn o V- 1 ! in (3.47). 

There are two leading terms in £> pn ; the first one is yU~ 2 r 2 7~ 4 ((2p + l)/i/?7) 2 ; since 
P < Co/(fihr) this term does not exceed C/x _2 7~ 2 . The second term is 0(fi~~ 1 hr 2/ -f~ 3 ) but 
looking its origin one can notice that it must contain | V^( V/^)! 2 ; then 

|VzH pn o vj/- 1 ! < /.x = CyT^r. 
Plugging £ = /i _2 7~ 3 r into (3.47)i results in 0(yU 1//2 /?~ 2 ). □ 

Remark 3.8. Note that in the extra terms 0(/i 2 /?~ 2 ) and 0(/i 3//2 /) _3 / 2_<5 ") appearing in the 
general case and under condition (0.8)^ the main contributor is far outer zone where 7 ~ 1 
and r Ri 1 while in the sharp estimate 0( / u~ 1 / 2 /7 -3 ) the main contributor is zone where 
7 « /i -1 / 2 and r ~ 1; 

(ii) Actually one can greatly improve remainder estimate in (0.8)i case using more refined 
partition-rescaling technique like in the proof of proposition 4.5. I believe that it can be 
even brought to 0(fi~ 1 / 2 h~ 3 ). However I think that really interesting are only general 
(q = 0) and generic (q = 3) cases. 

3.2 Estimates in the near outer zone 

After deriving remainder estimates in Z* ut I need to consider the remaining part of Z out 
(which will be the near outer zone 

(3.48) Z' out = Z out n {dist(x, E) < /1- 1 / 2 /)- 5 } 

and the inner zone. However I also introduce and start from the inner core where I am 
able to apply pretty non-sophisticated approach and still derive proper estimates, just to 
keep r more disjoint from to be able to use canonical form (2.38). 

3.2.1 In this subsection my goal is to derive remainder estimate 0(fj>~ 1 / 2 h~ 3 + fi 2 h~ 2 ). 
Therefore my approach here would be pretty unsophisticated. Note first that 

(3.49) The contribution of the domain 

(3.50) {dist(x, A) < r, dist(x, Z) < 7, \Z X \ < p} 

to the remainder does not exceed Cfih~ 3, -fr 2 p 2 (while the main part of asymptotics is given 
by (3.1) with 7" > Ch\ log/?|); 
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therefore considering inner core given by (3.50) with r = /i 1 / 3 + 2<5 ' ; j = pL 2 / 3 + 45 ' one gets 

This remainder estimate is very rough and should be improved despite rather unam- 
bitious goal of this subsubsection. However let me note first that considering the strip 
{| x i| < /i _1 / 2+3<5 '} and thus taking r = 1, p = 1 and 7 = / u _1 / 2+3<5 ' one would get the 
remainder estimate 0(/i 1 / 2+3<5 h~ 3 ) which is 0(/i 2 /?" 3 ) as as /i > h~ 2 ^ 3 ~ 5 . Further, the con- 
tribution of this zone to the asymptotics is 0(h~ Ar /r 2 ) and thus taking r = 1, p = 1 and 
7 = p 35 '" 1 / 2 one would get 0(/i~ 1 / 2+35 '/?~ 4 ) which is 0(fi 2 h~ 2 ) as p > h 5 ~ 4 / 5 . Therefore in 
what follows until the end of subsubsection one can assume that 



(3.51) 

Note first that 
(3.52) 



h~ s <p< h 



s-i 



\F t ^ h -i T x T (t)V(uQl)\ < Ch s 



as Q is supported in Z( 7 ) = f {|xi| < 7} and T < T < 7~i with To = Ch\ log h\ and T"i 
this statement is empty in the case of the very strong magnetic field 



e/i" 



(3.53) 



e(/?|log/7|)- 1 < n< Or 1 



which is not the case now. 
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However looking at the propagation in direction K4 as \Z\\ x p one can easily notice 
that the propagation speed is x p and the shift for time T is x p T; this shift is observable 
as p 2 T > Ch\\ogh\. Therefore (3.51) holds as 7" < T < T x with T = Cp~ 2 h\ log h\ 
and Ti = e/i _1 7~ 1 . Under assumption 7 > CoyU -1 ^ 2 this short evolution is contained in 
6x7- vicinity of the original point 25 ). 

These two intervals [Ch\ log h, eji^ 1 ] and [Cp~ 2 h\ log h\, e/i _1 7 _1 ] overlap as 

(3.54) p>(Cph\\ogh\f' 2 
and therefore I arrive to 

(3.55) Under assumptions (3.51) and (3.54) estimate (3.52) holds as T E [T , Ti] with 
7~o = Ch\ \ogh\, Ti = e min( / Li~ 1 7 _1 , /x" 1 / 2 ). 

Therefore the contribution of the domain 

(3.56) {dist(x,A)< r, dist(x, T) < 7, \Z X \ > C{ph\ log h\) 1/2 } 

to the remainder does not exceed Cp,h~ 3 r 2 7(7 + fi^ 1 ^ 2 ). 
Combining with (3.49) I arrive to 

Proposition 3.9. Contribution of domain (3.50) with p x 1 to the remainder does not 
exceed 

(3.57) C///?~ 3 r 2 7(7 + /i~ 1/2 ) + C/i 2 /rV 2 | log h\ 

while the main part of asymptotics is given by (0.12) with (any) Tq > Cfih\ \ogh\. 

Remark 3.10. Estimate (3.57) is sufficient for my limited goal here. In particular it covers 
the inner core giving remainder estimate 0(fi~ 1 / 2 h~ 3 + p 2 ^ 3+5 h~ 2 ) which is not only better 
than 0(/i~ 1 / 2 /j- 3 + n 2 h~ 2 ) but than 0(pT x l 2 h^ + /i 3 / 2 /r 3 / 2 ). 

So, I need to cover the near outer and inner zones; in both 7 < r 2 and therefore I 
can use canonical form (2.38). Even without it, according to section 1 the drift speed is 
x iT x p 2 r^i~ 2 as p 2 r > Q7. 

Then the shift for time T is observable as [i~ 1 p 2 r r y~ 2 T x p > Ch\ \ogh\; plugging T = 
/i~ 1 7~ 1 one gets (3.58)2 below and therefore in (3.52) one can upgrade T\ = efi~ 1/ -f~ 1 to 
7~i = e/ry 2 provided 

(3.58) i, 2 P 2 r > cj, p 3 r > Cp 2 ^h\ \ogh\. 

25%> One can improve this statement and the estimate following from it. 
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However, the contribution to the remainder of the part of Z' out where 7~i = ep7 2 is 
0(p" 1/2 /r 3 ). Therefore one needs to consider only part of Z' out where (3.58)i2 is vio- 
lated. At this moment one can ignore subzone {p < C(ph\ log h]) 1 ^ 2 } and therefore one can 
upgrade 7~i = ep" 1 to 7~i = ep _1 7 _1 with no penalty. 

Moreover, using 0(p,h~ 3 r 2 ^ 2 p 2 ) estimate (with r = 1,7 = /j -1 / 2 +3<5 ) one can see easily 
that the contribution to the remainder of the zone where condition (3.58)2 is violated does 
not exceed Cp 1 / 3 /^ 7 / 3-15 which c^p -1 / 2 /? -3 + p 2 /?~ 2 ). 

Therefore one needs to consider only zone where condition (3.58)i is violated; again due 
to 0(p/?~ 3 7 2 p 2 r 2 ) estimate its contribution to the remainder does not exceed C p -1 / 2 h~ 3 ~ 5 
which is only marginally short of what I want (and only the case p < /j~ 2 / 5 ~ 5 needs to be 
considered) . 

The proper estimate of the zone where condition (3.58)i is violated can be done easily on 
the base of proposition 1.6. Really, the drift speed is x p _1 (p+ p')Ar^ 2 with p' = O^r^ 1 ), 
A = \p — p'\ and therefore one can upgrade 7~i = e/i _1 7 _1 to T\ = ep7 2 as 

p,~ 1 (p + p')Arj~ 2 x /i _1 7 _1 x A > C/)| \ogh\ 

with p' = 0(7 1 / 2 r~ 1 / 2 ). Therefore one can take 

A = C^r' 1 ' 2 ^ log/)|) 1 / 2 (p + p / )- 1 / 2 ; 

then as \p — p'\ > A one can take 7~i = ep7 2 . 

On the other hand, the contribution of the zone {(x,£) : \p — p'\ < A} R Z^ r ^ to the 
remainder does not exceed 

Cp/rVpV'A = C/i 2 /)- 3 r 2 7 5 / 2 r 3 / 2 (/ ) | log h\f' 2 p n ' 2 . 

Since p' < c{^/r) l l 2 the latter expression does not exceed Cp 2 /?~ 3 7 l:L / 4 r 5 / 4 (/?| log/?!) 1 / 2 . In 
particular, the contribution to the remainder of the part of Z^ where (3.58)i is violated 
does not exceed Cp 2 /?~ 5 / 2 7 n / 4 | log/?! 1 / 2 and therefore the contribution of the corresponding 
part of Z' out does not exceed Cp 5//8 /7~ 5 / 2 ~ 5 . 
So I arrive to 

Proposition 3.11. Under condition (2.1) the total contribution to the remainder of zone 
Z' out does not exceed C(p _1 / 2 /?~ 3 + p 2 /)" 2 ). 

3.2.2 Now I need to improve the previous estimate under non- degeneracy condition. How- 
ever, let us summarize the remainder estimate I actually derived: 

(i) The inner core intersected with {p > C(p/)| log h]) 1 ^ 2 } contributed 0(ph~ 3: y 2 ?Q) = 
o(p~ l l 2 h~ 3 ) where 70 = p~ 2 ' 3+AS \ ? = p" l l 3+25 ' ; 
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(ii) Z' out n{p > C{fih\ log h\yi 2 } contributed 0(/x/r 3 7 2 (/i/> 7l ) 2 / 3 ) = O^ 1 ^^ 1 '^ 5 ") where 



71 = a* " 1/2+4a '; 



(iii) {|xi| < 7i } n {p < C(/x/?| log /7|) 1/2 } contributed 0(/i 2 /r 2 7l | log = 0(^/2+5" /,- 2 ); 

each of estimates in (ii),(iii) are less than 0(/i -1 / 2 /? -3 ) as fi < h~ 1 / 2+s . Therefore in what 
follows one can assume that 



-l 



(3.59) h- 1/2+s <fi<Ch 

The better analysis is based on the canonical form (2.38). However, this requires the inner 
core to be treated separately: 

Proposition 3.12. Under condition (3.59) the contribution of 

Z\ r) = {dist(x, E) < r 2 , dist(x, A) < r} 

with r > / u~ 1 / 3+25 ' does not exceed C[ih~ 3 r 6 while the principal part of asymptotics is given 
by (0.12) with T=T= Ch\ log h\ . 

Proof. Looking at the canonical form (2.17) and scaling x" i— > x" r~ , x\ \— > xir -1 , /th->/i t = 
/zr 3 , h i— > = hr -1 one gets 2-dimensional magnetic Schrodinger operator with parameters 
/i*, /)*, which also is 1-dimensional /14-PDO with respect to X4 with ^,4 = \x~ x hr~ 2 \ for such 
Schrodinger operators one can prove easily that the main part of asymptotics is given by 
Weyl formula and is of magnitude h^ 1 h~ 2 = [xh~ 3 r 4 while the remainder estimate gains the 
factor fall* = fihr 2 and thus is of magnitude fi 2 h~ 2 r 6 . Really, if /x^ = 1 it would follow from 
the standard results (and no non-degeneracy condition is needed); as p* > 1 one needs just 
apply extra rescaling x' 1— > x'p* which is also rather standard. □ 

Corollary 3.13. Therefore even without nondegeneracy conditions contribution of the (ex- 
tended) inner core Z^ with r = /i -1 / 4 to the remainder does not exceed C ji" 1 / 2 h~ 3 . 

Thus in what follows I can apply canonical form (2.38). Let us fix some index n and some 
final 6(0, 1) element where in contrary to subsection 3.1 final means only after reduction 
to (2.38). This element is either regular, when on the original B(x' , er) element 

(3.60) |\// : 2 - 1 -(2n+l) y u/)| + |Vy(\/^ 1 -(2n + l)/i/))| 2 x L 2 > r 2 
or singular when on the original B(x', er) element 

(3.61) I Vff 1 - (In + l)fih\ + IVx^V/T 1 - (2n + l)(ih) \ 2 <r 2 . 
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Since analysis of the regular elements is similar to the analysis of the singular elements but 
simpler, I consider only singular elements, leaving regular elements to the reader. 
Let us introduce the scaling function with respect to (x3,x 4 ,£ 4 ) 

(3.62) £ = £ n = er- 1 ^,^ WT 1 ! + I 

where £ > C'~yr~ 1 \ \ogh\ will be chosen later. With respect to \Zf\ I use scaling function 
p = e(|Zi| + r 1//2 £) . Then ellipticity of A n is not broken unless ( Vf^ 1 — (2n + l)ph) x p 2 or 
p < cr 1 l 2 £. 

Then it follows from the standard theory that for u n holds as Q is supported in the 
corresponding ^-partition element, T e [To, T{\, r < Cr£ 2 with 



(3.63) T Q = Ch(p 2 + r£ 2 Y l \\ogh\, 

(3.64) 7"! = ep' 1 ^ 1 , 

(3.65) p + r 1/2 £ > Cp = C{fihj\ log h\) 1/2 
and moreover, under these conditions 

(3.66) \F t ^ h -i T x T {t)V{u n Qy ] )\ < Cph~ 2 r 3 £ 3 1 
Therefore 



(3.67) While the contribution of this {£; p)-subelement into asymptotics does not exceed 
Cph~ 3 r 3 (p 2 + r£ 2 )£ 3 7 its contribution to the remainder so far does not exceed 

Cph~ 3 r 3 (p 2 + rf)^ x T T-f 1 = Cp 2 h~ 2 r 3 £ 3 1 2 . 

Note that for each {£; p)-partition subelement in (x 3 , x 4 , £ 4 ; £i, £ 3 ) the number of indices n 
satisfying (3.62) and also violating ellipticity does not exceed 

(3.68) M = (Co^-V + rP) + l) 

Also note that under condition (0.8)^ the total number of £-subelements is 0(£ q ~ 3 ). Then 
the total contribution of all ^-elements to the asymptotics (as magnitudes of r, 7 are given) 
does not exceed 



(3.69) 



Cph- 3 1 r i - 2 ~ q)+ {p 2 + rt 2 ) x (c (/i/?)~V + r£ 2 ) + l)xf 
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which one can sum nicely to its value at maximal £ = 1, p = 1 and then to maximal r = 1, 
7 = resu iting OrV~ 1/2 . 

On the other hand, the total contribution of ^-subelements to the remainder so far does 
not exceed 

(3.70) C/i 2 /)- 2 7 2 r (2 ~ q)+ x {Co^hy^p + rif + l)xf 

which one can also sum nicely according to the same principles but the sum is C/x/? _3 7 2 
not estimated that well and £, r close to 1 are problematic. So I want to increase Ty. 

I remind that W n = cu^Vf^ 1 — (2/7 + l)ph) with uj = (o-ai)^ 1 f 2 and here (see proposition 
2.8) to ~ \{Zi, fi}|~ 2 / 3 rescaled; therefore for "small" p calculation of V1/V n and u;V( Vff 1 — 
(2n + l)p,h) are equivalent modulo a;p 2 r _1 . Note that 

(3.71) As (Vf^ 1 — (2/7+ l)ph) x p 2 and £ > Cp 2 ^ 1 the propagation speed with respect to 
x 4 does not exceed CyU _1 r~ 1 £ and is of this magnitude as |^ 4 l/l/ n | x l\ then the shift with 
respect to x 4 is observable as 

(3.72) p^rHT x£> Cji^r 2 ^ log h\\ 



(3.73) As (Vf^ 1 — (2/7+ l)ph) x p 2 and £ > Cp 2 ^ 1 the propagation speed with respect to 
p,~ 1 hD" does not exceed CpT 1 r~ 1 £ and is of this magnitude as |<9 X »I/1/„| x £; then the shift 
with respect to p,~ 1 hD" is observable under the same condition; 

(3.74) Finally, in the short-term propagation (as T < e/i -1 7 _1 ) the propagation speed with 
respect to (xi, X3) is of magnitude \Z\\ x p and the shift is observable as pT x p > Ch\ log h\. 

On the other hand, 

(3.75) Let p 2 > C7; then the drift speed (in the rescaled coordinates) with respect to x" is 
of magnitude C/i~ 1 7~ 2 p 2 and the shift for time T is observable as 

(3.76) ^VVT'x P > Chrl \ 
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Therefore one can increase T\ = ep lr y 1 provided one of the conditions (3. 72), (3. 76) 
holds with T = T\. Plugging T\ = e/i~ 1 7~ 1 into (3.72), (3.76) I get 

(3.77) p > p d = f C{p 2 h 1 'r 1 \ \ogh\) 1/3 + C 7 1/2 , 

(3.78) £ > £ = C^hjr' 1 ] log h\) 1/2 + CpV 1 

respectively. 

Now consider increased T\. I claim that 

(3.79) Under condition (3.78) one can take T = ep 7 2 £ 

Really, then magnitude of |V XliX4 ,f 4 W n | changes by no more than e£ as long as p < e\\ 
however, as p > e condition (3.77) is fulfilled and one can take T2 = ep7 2 ; 

Therefore the contribution to the remainder of the zone with the given magnitudes 
of p, £, 7, r satisfying (3.78) does not exceed expression (3.69) multiplied by TqT^ ; the 
resulting expression does not exceed 

(3.80) C/i/rV^+tp 2 + rf) x (Co^rV + r£) 2 + l) x £ q x 

h(p 2 + rf)- 1 x fji-^r 1 = 

C/i^/rVV 2 -^ ((p + ri) 2 + iih)^- 1 ; 

then summation with respect £, p results in Cp _1 /7 _3 7 r as q > 1 (then one should take q 
slightly larger than 1). As q = 1 one gets Cp _1//2 /? _3 (l + p/?| log h\) as provided summation 
is taken over p < c(£ + | log /jj^ 1 ); but the latter case is covered by arguments below. 

(3.81) Under conditions (3.77) and p > r l / 2 £ one can take T 2 = ep 7 2 pr^ 1 / 2 . 

Really, in this case due to bounds for propagation and drift speeds this guarantees that 
both I V XliX4i £ 4 W n \ and W n x \Zi\ 2 change by no more than r~ x ' 2 p, e x p 2 respectively during 
the evolution. 

Now instead of (3.80) one gets 

(3.82) Cp/7-V (2 - q)+ (p 2 + r£ 2 ) x (c {ph)-\p + r£) 2 + l\ x £ q x 

h{p 2 + rf)' 1 x n-^p- 1 < 

Cp^/rVV 2 "^ ((p + r£f + p/i^V 1 

and summation with respect to p, £, 7, r results in Cp" 1 / 2 /?^ 3 . Therefore 

(3.83) The total contribution to the remainder of elements in Z'^ut satisfying (3.65) and 
one of conditions (3.77), (3.78) does not exceed Cp -1 / 2 /? -3 . 
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3.2.3 So, I am left with two types of elements, which are not covered by the above 
arguments: 

(3.84) With p <P<P = C{p 2 h^r^\ log/?]) 1 / 3 + C7 1 / 2 and £ < fx r l pQ + p 2 r~ x \ 

(3.85) With p < p = C(/i/) 7 ) 1 / 2 and £ < £ x r~V 2 p Q . 

For (3.84)-type elements r£ 2 < p 2 and expression (3.69) becomes 

Q//rV (2-£? V x (V/?)-y + 1) x 

multiplying it by TqT^ 1 = Chp~ 2 \ log /j| ■ ^7 one gets 

Cp 2 h- 2 -f 2 r {2 - q)+ x ((/i/?)~V 2 + l) x ^1 log/ij 
which after summation over p becomes 

(3.86) Cp 2 h~ 2 7 2 r (2 - £?)+ x ((/i/>)~V 2 + l) X £"| log /»| . 
Meanwhile for (3.85)-type elements expression (3.69) becomes 

(3.87) Cp J h~ 3 jr^ + pl£ q . 
Let us plug p, po, £ and 7; one can take 

T = ^-V2 r l/2 ( ^ = ( ^ r l/2 ) l/2 i - = ^1/2^1/2)1/3 + ^l/4 r l/4 ^ 

£ = {fi^hr 1 ' 2 ) 1 / 2 + (p 1 / 2 /)) 2 ^" 2 / 3 + pT^r 1 ' 2 , 

and compensate this by multiplication the result by /? _<5 " with an arbitrarily small exponent 
<S">0. 

Take first £ = 4 = {p}l 2 hr x l 2 ) x l 2 \ note that then in (3.86), (3.87) factor r~ q ' A coming 
from £ q is more than compensated by the factor r coming from 7 2 or 7P0 respectively. 
Therefore summation with respect to 7, r results in the same expressions calculated as 
r = 1, 7 = / i-V2 ) £ = p = (^1/2^1/2 and p = (^ h y/3 + ^-1/4. 

(3.88) C/)- 3 ((^ 1 / 2 /)) 2 / 3 + ^ + phYp^hf' 2 x /T 5 ", 

(3.89) C^ 2 h-\p l l 2 h)^' 2 x /T 5 ". 

One can see easily that both these expressions do not exceed 
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(i) C(/i~ 1 / 2 /)- 3 + ^ 3 / 2 /T 3 / 2 ) as q = 1; 

(ii) C(/i- 1 /2/,-3 + ^1/2+^-2) as g = 2 ; 

(iii) C/i~ 1/2 /r 3 as g = 3. 

However one should consider (3.86) with £ = £ 2 = (/i 1 ^ 2 /?) 2 ^ 3 r~ 2 / 3 and £ = £ 3 d = /i" 1 / 2 /- -1 / 2 ; 
as g = 1 one can take again r = 1 obviously; as q = 2 it is not so but one can take then by 
<7 = 9/5 arriving to 0(yU~ 1 / 2 /7 -3 ). 

Then combining with (3.83) I arrive to 

Proposition 3.14. The total contribution to the remainder oj the zone Z' out does not exceed 
C^t -1 / 2 /?" 3 + C/i 3 / 2 /7~ 3 / 2 under condition (0.8)\, C/i _1 / 2 /)~ 3 + C/i 3 / 2 /)~ 1_<5 under condition 
(0.8)2 o,nd C/i -1 / 2 /) -3 under condition (0.8)3. 

Remark 3.15. The main part of this asymptotics is given by expression (0.12) with T = Tq. 

3.3 Estimates in the inner zone 

Now I want to consider the inner zone 

(3.90) Z- mn = {dist(x, Z) < 70 = Cfi~ 1/2 ■ dist(x, A) 1 / 2 }, 
in term consisting off the inner bulk zone 

(3.91) Z innJ i = Z- mn n {\Zx\ -dist(x,A) < e/i • dist(x, Z) 2 } 
which mimics the outer (bulk) zone Z out and the true inner zone 

(3.92) Z mnJ = Z- mn n {|Zi| ■ dist(x, A) > efi ■ dist(x, Z) 2 } 

which mimics the inner zone for 2-dimensional Schrodinger operator. One needs to remem- 
ber that trajectories can leave one of them and enter another one after a while. 

Note first that the inner core is already covered by proposition 3.12, so in the definition 
of the inner zone I can include condition dist(x, A) > /i" 1 / 4 and therefore use canonical form 
(3.28) there. Actually using the same method as in the proof of proposition 3.12, one can 
provide estimate 0(/i/?~ 3 r 2 7 2 ) = 0(h~ 3 r 3 ) thus covering even Z- mn n {dist(x,A) < /i~ 1//6 } 
but I don't need this. 
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3.3.1 I start from the inner bulk zone Z- Inn ji. Let \Zy\ x p, pr < ep^ 2 and condition 
(3.77) be fulfilled. Then exactly as in the outer zone the drift speed is x pT x p 2 r^~ 2 < 
cp and evolution speed with respect to /7D4 does not exceed c; then for time Ty = epr 
dynamics remains in Z mn \\. Therefore in the domain where condition (3.77) is fulfilled and 
P > C{p,h\ log/?!) 1 / 2 one can take Ty = epr and T = Ch\ \ogh\; however in the estimates 
this logarithmic factor is not needed. 

Therefore the contribution to the remainder of this domain intersected with Zr rilfP \ does 
not exceed Ch~ 3 ^fT^ 1 r 2 p 2 = Cp~ 1 ^ 2 h~ 3 r 2 p and the summation with respect to (r,7,p) 
results in 0(/i~ 1 / 2 /) -3 ). 

On the other hand, consider the domain where condition (3.77) is violated, but still 
P ^ Pi- The contribution to the remainder of of this domain intersected with Z( rrlp \ does 
not exceed Cph~ 3r y 2 r 2 p 2 where two factors 7 come as the width of the strip and the part 
of 7"-, -1 x ^7. The summation with respect to (r, 7, p): rp 2 < cy results in Cp^^bT 3 . 

Further, consider domain where condition (3.78) is violated but still p > p~y. In the sim- 
ilar manner one can estimate its contribution to the remainder by Cph~ 3 ^ 2 (p}l 2 h\ log h\) 2 ^ 3 
which is o(p 2 h~ 2 + p~ l l 2 h~ 3 ) for sure and is 0(/i~ 1//2 /)~ 3 ) as p < C(h\ log /?|)~ 4//5 . 

Finally, contribution to the remainder of the domain where p < pi does not exceed 
Cp 1 p 2 = Cp 3 ' 2 h~ 2 \\ogh\. 

Thus I arrive to 

Proposition 3.16. Contribution of Z- mn u to the remainder is 0(p 2 h~ 2 + p^^h^ 3 ) for 
sure. Furthermore, this contribution is is 0(p~ 1 ^ 2 h^ 3 ) as p < C{h\ \ogh\)- 1/2 . Again the 
main part of this asymptotics is given by expression (0.12). 

3.3.2 Again, the estimate achieved in proposition 3.16 is good in the general case but 
should be improved under non- degeneracy condition as 

(3.93) c{h\ log/7|)~ 1/2 < p < Ch- 1 . 

Again one should consider regular final elements and singular ones and I consider only 
latter. Further introducing the scaling function £ n by (3.62) I can define 7"o and T\ by 
(3. 63), (3. 64) as (3.65) holds; furthermore estimate (3.66) holds. 

Further, if condition (3.78) is fulfilled, I can increase T-y to 7~2 = el since on this time 
interval both |V Xl , X 4,£ 4 | change by no more than ey min(£, r). 

Then repeating the arguments of subsubsections 3.2.2 I find that the contribution of 
this part of Z-' nn u to the remainder does not exceed expression similar to (3.80) but with 
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p 7 2 replaced by 1: 

(3.94) Cfih-^r^+ip 2 + r£ 2 ) x (Co(/i/j) _1 (p + rl) 2 + l) x £« x r;(p 2 + rf 2 )" 1 x r 1 = 

C/rV (2 ~ ,)+ ((p + rl) 2 + (thy*- 1 ; 

Then the sum with respect to (I, p, r)-partition with p > £ does not exceed Ch~ 3/ y as q > 2 
and C/j~ 3 7| log /i| as q = 1 and one can get rid off logarithmic factor exactly as in the 
analysis in Z'out, II. 

Furthermore, condition (3.77) is fulfilled instead of (3.78) and p > R lj/2 £, I can increase 
7~i to 7~2 = ep and the contribution of this part of Z[ nn u to the remainder does not exceed 
expression similar to (3.82) but with p _1 7~ 2 replaced by 1: 

(3.95) C/i/7-V (2 ~ q)+ (p 2 + r£ 2 ) x (c {ph)-\p + rl) 2 + l) x £ q x h{p 2 + rf)" 1 x p' 1 = 

Ch~ 3 jr {2 - q)+ (Z + ^jiv 1 ; 

then summation with respect to (I, p, repartition results in Ch~ 3/ -f + C/i/? _2 7| log /j[ . There- 
fore the sum with respect to 7-partition with p > £ does not exceed C/! -1 / 2 /?^ 3 as q > 2 
and n~ 1/2 h~ 3 + C/j, 1/2 h~ 2 \ \ogh\ as q = 1. 

On the other hand, due to the same arguments as in subsubsection 3.2.3 the contribution 
to the remainder of the domain where both conditions (3. 77), (3. 78) are violated but 1 (3.65) 
still holds does not exceed (3.86). Then the summation with respect to (£, p, r, 7) partition 
results again in Cpih~ 2 {pi l l 2 h\ \ogh\) q ' 2 + 0(/i- 1/2 /r 3 ). 

Finally, contribution of the domain where (3.65) fails is tackled in the same way as in 
subsubsection in subsubsection 3.2.3 and estimated by (3.87). 

So, I arrive to the following copy-cat of proposition 3.14 and remark 3.15: 

Proposition 3.17. The total contribution to the remainder of the zone Z mn ji does not 
exceed CpT x l 2 ln~ 3 + Cp?l 2 h~ 3 l 2 under condition (0.8)i, CpT 1 ! 2 ^ 3 + Cp?l 2 r)^ x \ log/?| under 
condition (0.8)2 and CpT x l 2 h~~ 3 under condition (0.8)3. 

Again the main part of this asymptotics is given by expression (0.12). 

3.3.3 Now let us consider the true inner zone Z mn j defined by (3.94). The crucial differ- 
ence between this zone and Z out D Z mn ji is that in the classical evolution dist(x, Z) is not 
necessarily preserved 26 -* 

26 ) As |Zi| • dist(x, A) > C/i ■ dist(x, Z) 2 it is not preserved for sure; in the rest of this zone variation of 
dist(x, Z) is of the magnitude of dist(x, Z). 
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I remind that this subzone (as the whole inner zone) should be studied as r = dist(x, A) > 
c// -1 ' 4 . Also I remind that the contribution of subzone {\Zi\ < C(ph log h) 1 ^ 2 } to the 
remainder does not exceed C/i/) _3 7 x ph\ \ogh\ = Cp, 3 ^ 2 h~ 2 \ \ogh\) which is 0(fx 2 h~ 2 ) for 
sure and 0{p~ l / 2 h~ 3 ) as p < c(h\ log A?| ) x / 2 . 

Let us introduce 7 = /i _1,/2 (rp) 1//2 ; then in the classical evolution 7 is the magnitude 
of maxdist(x, T.) along the cyclotron movement. Note that the contribution of subzone 
{(/ih\ log/?!) 1 / 2 < \Zi\ < p} to the remainder does not exceed C\xh~ 3 ^ 2 p 2 = Ch~ 3 p 3 which 
sums to 0(p~ 1 / 2 h~ 3 ) in subzone {p < cp~ 1 ^\ log/7| -1 / 3 }. 

Consider first the general case of q = when our goal is 0(p~ 1 / 2 h~ 3 + ji 2 h~ 2 ) estimate. 
Then the total contribution of Z;„„j PI {\Z\\ > C(p,h\ log/?!) 1 / 2 } to the remainder is 0(/r 3 ). 
Therefore in this case one needs to analyze only p, < h~ 1 / 2 and only subzone {p > c(p~ 1 / 6 + 
p 2 ' 3 ^' 3 )}. 

I am going to prove that the contribution of Z- mn / fl {p > C(ph\ log Z?]) 1 / 2 } to the remain- 
der does not exceed Qi -1 / 2 /?" 3 + C(p}^ 2 h\ log h\y^ 2 h^ 3 (see proposition 3.19) but I prefer a 
bit less direct way to do this; I am concentrating on this zone due to the previous remarks 
and because I do not have non-degeneracy condition at the moment. Thus one can take 
7~o = Ch\ \ogh\ and T x = eyU -1 / 2 7 _1 = ep~ 1 (pr)~ 1 . Let us try to increase it. 

According to the series of propositions of subsubsection 1.4.3 the drift speed in the 
inner zone is ft/"p(£ 3 — k*p) + 0(/i -1 / 2 ). Then the shift during one cyclotron movement is 
x r P(G — k* p)T\ provided 

(3.96) rp|£ 3 - k*p\ > Cp- 1/2 
and it is observable as 

rp|£ 3 - k*p\ x p,- 1 / 2 x |£ 3 - k*p\ > Ch\ log h\ 

or equivalently 

(3.97) |£ 3 - k*p\ > C{rp)-W{jihnt\ \ogh\f/ 2 = C^p)-^^ log/))) 1 / 2 . 
Both conditions (3.96), (3.97) are fulfilled as 

(3.98) |6 - k*p\ > A d = f cirpyV^hl log h\fl 2 + cp'^rp)- 1 ' 2 . 

Then as (3.98) is fulfilled 27 ^ one can upgrade 7~i to T 2 = er|£ 3 — k*p\ because then not 
only magnitude of p but also of |£ 3 — k*p\ will be preserved during the evolution. However 

27 > This makes sense as A < p: i.e. 

(3.99) p > p 2 d ^ f r-^ 5 (p^ 2 h\ log /7|) 2 / 5 + ^V 1 / 3 ; 
otherwise one needs to take A = cp, thus redefining A r = ef min(A, p-i). 
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taking the evolution in the right time direction it could be increased to T 2 = er|^3 — /c*p| x 
exactly as in [Ivr6]. The standard calculations imply 

(3.100) The contribution of subzone of Z mn j where condition (3.98) is fulfilled intersected 
with {p > C(ph\ log/?!) 1 / 2 } to the remainder does not exceed Cp}l 2 br 3 . 

So, one needs to study near periodic zone 

(3.101) Z per = {|£ 3 — k*p\ < A}, 
intersected with {\Zi\ x p} fi {dist(x, A) x r} as 

P > c(fih\\o g h\y/ 2 + p 2 . 

Note however that as r < e then |Vu;| x r" 1 ^ and | Hesscj| x r~ 2 uj and therefore in the 
rescaled coordinates | V X3 , X4i ^ 4 W n \ x p 2 as r < e and p 2 > Cr and | Hess X3;X4 ^ 4 W n \ x p 2 
as r < e and p > Cr. In the first case periodicity is broken for sure and in the second 
case periodicity is not broken only as |Vl/l/ n | < C(p'jh) 1 ^ 2 which in intersection of Z per has 
measure not exceeding Cp}l 2 h^ instead of C(p}^ 2 h) 1 ! 2 ^ and since one can take 7~i = e/i~ 1 7~ 1 
here for sure the contribution of the reduced periodic zone to the remainder estimate does 
not exceed Ch~ A x (/i 1 ^ 2 /?) log/?|)7 x /j/77 < Cp}l 2 br 2 \ log /?| thus implying 

Proposition 3.18. The contribution of Z mn / fl {\Z-y\ > C(jih\ log h\) 1 / 2 + Cr} to the re- 
mainder does not exceed C/i" 1 / 2 /) -3 . 

Again the main part of this asymptotics is given by expression (0.12). 

So, only p < Cr and thus 7 < CpT x ' 2 r should be analyzed. 

Proposition 3.19. The contribution of Z- mn / fl {\Zi\ > C(ph\ log/?!) 1 / 2 } to the remainder 
does not exceed C/i~ 1/2 /r 3 + C(p 1,2 h\ log /7|) 1/2 /r 3 . 

Again the main part of this asymptotics is given by expression (0.12). 

Proof. Now I am concentrating only on even more reduced zone. In particular (3.99) implies 

r > (yLi 1 ^/,) | 0g ^1)1/3 + M -l/4 

One can take T\ = ep~ 1 / 2 (rp)~ 1 / 2 as condition (3.98) fails. Then the contribution of 
Z aer to the remainder does not exceed 



(3.102) C^/T^VpA = Ch~ 3 r 3 p 2 A = 

Ch^r 3 p 2 [c(rpY 1 l\p l l 2 h\ log/,)) 1 / 2 + c^ 2 (rp)^ 2 ) 
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Obviously when taking sum with respect to (r, p) partition the latter expression sums to 
itself with A = A\ r=l p=1/y=fl -i/2: 

(3.103) A = C(/i 1/2 /?| log h\) 1/2 + C/i~ 1/2 . 

Then (3.102) becomes O^" 1 / 2 /)- 3 + h~ 3 (p, 1 / 2 h\ log ft]) 1 / 2 ). 

One can check easily that the contribution of subzone {/?2 > p > C(p,h\ log/?|) 1//2 } is 
just 0{jih~ 3 ^ 2 p\) which is less than this. □ 

3.3.4 Again I would like to improve the remainder estimate achieved in proposition 3.19 
under non- degeneracy condition (0.8) q with q > 1 also getting rid off assumption \Z\\ > 
C(ph\ log/?|) 1//2 . In the analysis below instead of W n I consider its pullback u(V — (2n + 
l)?2/i/j) in the rescaled coordinates x' on Z; it is useful for Hessians but it transforms into 
properties of W n since I am interested only in the measures of sets |VtV„| x A. 

Proposition 3.20. As p, < ch" 1 ] \ogh\~ K the contribution of Z mn j PI {\Zi\ > Cr} to the 

remainder does not exceed C fi" 1 ^ 2 h~ 3 . 

Proof. Let us introduce 

(3.104) i = £ n = er- 1 |V X3 , X4 , e4 W n \ + I, 
and 

(3.105) Q = 0n = e\&-k*Wy 2 \+A. 

Then the contribution of (r, p, £)-elements to the asymptotics does not exceed (3.69)-like 
expression 

(3.106) Cp\T 3 1 r\p 2 + rf) x [C {ph)-\p+ rf) + l); 

if p 2 + ri 2 < Cpih^l \ogh\ then also r 2 < Cplry] \ogh\ and (3.106 does not exceed 
Cp 3 /j _1 7 3 | log/)| 2 < Cp 3 ^ 2 h" 1 \ \ogh\ 2 . So let us consider p 2 + ri 2 > Cp,rry\ \ogh\. 

As p 2 > Cr then periodicity is broken for sure (since Vl/1/„| > ep 2 — Cr and one can 
instantly upgrade 7" = Ch\ \og h\(p 2 + ri 2 )" 1 to T 2 = e. Otherwise it could be done as 
i > C(p,h\ log/?!) 1 / 2 . Then contribution of such elements to the remainder does not exceed 

(3.107) C^/rW + rl 2 ) x (c (/ih) ~\p + ri 2 ) + l) x h\ log h\{p 2 + ri 2 )- 1 = 

Cph^r 2 x (Co^rV + ri 2 ) + l) | log h\ 



3 ESTIMATES 



69 



because | Hess W n \ > ep 2 — Cr 2 and one can delete safely logarithmic factor because it is 
not in the estimate of the Fourier transform. 

Otherwise one can take T\ = ep,~ 1/ ~f~ 1 and then one gets 

(3.108) C/i/rVV + rf) x ^C (pih)' 1 (p + r£ 2 ) + l)xh\ \og h\(p 2 + rl 2 )" 1 ^ x £A = 

Cfi 2 h~ 2 ^ 2 r 2 x (C {fih)- 1 (p+ ri 2 ) + \ogh\ < Ch~ 3 x (fi 1/2 h) = 

Cfi 1/2 h~ 2 \ \ogh\ 2 

where one copy of (p}^ 2 h\ log h]) 1 / 2 comes as the width of zone with respect to £3 and another 
as the width with respect to |VI/I/„|. □ 

Proposition 3.21. Under condition (0. 5)i be fulfilled. Then 

(i) The total contribution of Z- mn / H {\Zi\ < er} to the remainder does not exceed CpT x l 2 h~ 3 
as pi < h' 1 ] log h\~ K . 

(ii) The total contribution of Z\ nn ,i H {|Zi| > er} to the remainder does not exceed 

(3.109) C/i" 1 / 2 /)" 3 + C^ A h~ 3 ' 2 \ log h^ 2 - 

in particular it is 0(p,~ 1 ^ 2 h~ 3 ) as p, < /?~ 6 / 7 | log/?|~ 2 / 7 28 ). 

Again the main part of this asymptotics is given by expression (0.12). 

Proof. Note that the most delicate is the case er < p < cr when even condition (0.8)2 does 
not yield (0.8)i for W n 2 ^ 

(a) Consider first r < ei, p x r, 7 x /i~ 1 / 2 r 1 / 2 p 1 / 2 x pr x l 2 r. 
Let us introduce a scaling function 

(3.110) e n = er 1 \VW n \ + ^l £ = C(/x 1/2 /7| log /7|) 1/2 

and the corresponding subpartition as p < cr. Note that actually I > I = C(p , yr~ 1 h\ log h]) 1 / 2 
is the condition that for time 7~i = e/i _1 7 _1 the shift is observable. Then exactly as before 

(3.111) The total contribution to the remainder of all ^-elements with i > £, p x r does 
not exceed CpT x l 2 h~ 3 . 

28 ) Actually I am going to prove later that one can drop logarithmic factors. In this case estimate cannot 
be improved without either additional assumptions or correction terms. 

29 > However (0.8)i with the first and the second order derivatives taken only with respect to A implies 
such condition for W n thus eliminating vicinity of A from the contributing more than C/i~/ 2 /?~ 3 to the 
remainder. 
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So only £ < £ subelements should be considered. Note that after rescaling |Va>£| x 1 
and therefore 

(3.112) For each I subelement there exists n such that |Vu;(V — (2n + l)f 2 pih\ > £ as 
\n— n\ > C£(/i/)) _1 and therefore over each ^-element there are no more than M = C£(fih)~ 1 
other subelements. 

Therefore the total contribution of them to the remainder does not exceed 

(3.113) C^VV + rf) x (c (lih)- 1 I+ l) x h(p 2 + ri 2 )" 1 x ^ 7 A x 

C^/rV x (P(nh)- 1 + e\ 

since A x i. After summation over r-partition one gets C/i 1//2 /) _3 | log h\ + Cph~ 2 £ which is 
exactly (3.109). 

(b) Consider now r x 1, p x 1, 7 x pr x l 2 . One needs to recover an extra factor fih 
in comparison with estimate of proposition 3.19 and only case p, > h s ~ 2 ^ 3 needs to be 
considered. 

Let us assume first that 

(3.114) I Hess(-^)| > e 

'2 

and let us fix 

(3.115) (>(= C(p 1/2 h\\ogh\y +enh 
and examine subzone 

(3.116) Q c = {3neZ + : | Hess W„\ x £}. 

Let us introduce (^-admissible partition in this zone; as r < ( one should replace it by 
r-admissible partition. Note that one can assume that 

(3.117) For each element of such partition no more than M = CoC(pih)- 1 magnetic numbers 
n satisfy condition | Hess W n \ x (. 

Really, if |Vo>£| x 1 then arguments of (a) work. If |Vcj^| x< e 1 but |\7/^ — 1 \/| x 1 then 
jVo'V'l x 1 and then |VI/1/„| x 1 for all n < C(fih)~ 1 and everything is easy. So, let both 
|Vu;£| x< ei and |V/2~ 1 V| < e±. Then as | Hessa;fe| x 1 then (3.117) is obviously true. 
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Finally, if | Hesso;^| < ei then (3.114) implies that | HesscL>\/| x 1 and then | Hess W n \ x 1 
for all n < C{fih)-\ 

Then for each index n described (3.117) all the above arguments could be repeated but 
with 

(3.118) £ = C x \VW n \ + ^i £ = C(C V>l bg />|) 1/2 
leading to the contribution to the remainder 

(3.119) C/i 2 /)~V x CC^r 1 x h\ log /i| x M x Qi 1/2 /T 2 | log h\ 2 
since (£ = A. 

So I am left with ( = ( and only with indices n | Hess W n \ < (. Here one needs to 
consider only subelements with £ < i i.e. with 

(3.120) |Vl/l/„<CC 2 . 

So, I am looking on ( partition elements satisfying this condition. 
Assume first that | Hessu;/^| x 1. Let us introduce scaling function 

(3.121) £=\\/ujf 2 \ + ^l 

then as £ > £ over each such subelement leave no more than C( 2 /(fih£) indices n satisfying 

(3.121) and then (3.119) is replaced by 

(3.122) C/i 2 /?- 3 7 2 x C^'V^r 1 x h\ \ogh\ x A x Ch^A^i'^ \ogh\; 
alternatively one can replace (3.119) by 

(3.123) C/i 2 /7~V x ({fih)- 1 x h\ \ogh\ x £A x Ch~ 3 A(£\ \ogh\ 

where £ is a width of of the strip where |Va>/^| < £. 

Comparing (3.122) and (3.123) one can see that the best choice of £ is £ = C 1 ^ 2 ; then 
both of them become C/?~ 3 A^ 3//2 | \ogh\ x C/x 1//2 /j~ 3 | \ogh\ 2 . 

On the other hand, as | Hesscu/^l < ei then due to the arguments proving (3.117) this 
analysis is not needed. 

On the other hand, if (3.114) is replaced by condition 

(3.124) | V (^)|>e 

'2 
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the analysis is essentially the same, but simpler, with (3.116)-(3.117) replaced by 



respectively and in (3.118) Hess W n \ x £ is replaced by |VW„| x (. 

(c) Finally as p < er condition (0.8)i translates for a similar condition to W n and then one 
can apply the arguments similar to those of proposition 3.20 which imply (i). I leave easy 
details to the reader. 

(d) Analysis of p 2 < c(ph\ log/?!) 1 / 2 is also easy. Let us introduce £ by (3.110). Then as 
p 2 + r£ 2 < Cp}! 2 h\ log/?| factor (p 2 + r£ 2 ) takes care of everything; otherwise factor A is 
retained but factor (p 2 + r£ 2 ) translates in the end into extra factor ph\ log/7. I leave easy 
details to the reader. □ 

3.3.5 Now the better remainder estimate should be pursued only as p 
Further, as q = 0, 1 the contribution of Z per to the remainder is better than the remainder 
estimates 0(p~ 1 ^ 2 h~ 3 + p 2 h~ 2 ), 0(p~ 1 ^ 2 h~ 3 + p 3 ^ 2 h~ 2 ~ 5 ) and therefore only q = 2, 3 should 
be considered. 

Moreover, it follows from the proof of proposition 3.21 that if summation is taken 
over Z per R {dist(x, A) < r} then the second term in (3.109) would get an extra factor r 4 
(from r 2 7 2 ) and therefore only zone Z per n {dist(x, A) >? = p~V 16 \ log/j]- 1 / 8 } should be 
considered. 

Moreover, then ( x ph (which is much larger than the first term in (3.114) and therefore 
for each ^-partition subelement only one magnetic number n should be considered. 

There is a minor problem as p > /j^ 1 ! \ogh\~ K because then contributions of {p < er} 
and {p > Cr} were estimated by p}l 2 br 2 \ log/)| K ; however in this case there also factor r 2 
at least coming if looking at Z per n {dist(x, A) < r} and therefore one needs to consider 
only Z per R {dist(x, A) > | log h\~ K }. Then in the first case condition (0.8)2 supplies factor 
(p}/ 2 h\ log/?) 1 / 2 ! \ogh\ K which is more than enough to take care of | \ogh\ K ; however I will 
just include both cases in the final analysis, 

Then repeating all arguments of subsections 2.8, 2.9 of [Ivr6] with the obvious modifi- 
cations, due to the presence of r (which is greater than 7 > r) and of (x2, p~ 1 hD2) in all the 
operators, one can prove easily the following 

Proposition 3.22. Under condition (0.8)2 the total contribution of Z mn j R {\Z\\ > e} to 

the remainder does not exceed Cp~ x l 2 h~ 3 while the main part of the asymptotics is given 



(3.125) 
(3.126) 



C>C = c(// 1/2 /j|l0g/7|)l +ph, 

ft c = {3neZ + : \VW n \ x(}, 
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by (0.12) plus a correction term 



(3.127) J £^(x',r)^(x')dx', 



where temporarily 



(3.128) 



£, 



corr 



MW def 



h- 1 




1 



and in this formula one can choose arbitrarily 



(3.129) T G [Ch\ log/iI.eV^V 1 ]. 

(3.130) T> Cr- 1 /2(/ 7 || g/ J |)V2( max( | 6 _ r | Zl || i/7 : 



/2 )) 



-l 



and pseudo-differential Q = I in Z per n {dist(x, A) > r} and supported in the same zone, 
slightly inflated. 

In the next section this correction term will be rewritten in more explicit form. 

3.3.6 Conclusion. So, with the remainder estimates described in theorems 0.1, 0.3 
asymptotics are derived where the main parts are given by rather implicit formula (4.1). 

The analysis in frames of theorem 0.5 is similar but much simpler. I leave it to the 
reader. 

4 Calculations 

So, according to the previous section, with the remainder estimates described there the 
main part of the asymptotics of J e(x, x, 0)^(x) dx is given by 



where Q L form an appropriate partition of unity and 7~ t > Tiaq a ls° described there. 

In this section I will make all the calculations transforming implicit expression (4.1) to 
a more explicit form, namely to 



There are many methods to make a reduction and I will employ them depending on the 
power of the magnetic field and non- degeneracy condition. 
Also I will calculate (3.122) in a more explicit form. 



(4.1) 




(4.2) 
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4.1 Moderate magnetic field. I 

I remind that I proved that for [i < h~ 5 the remainder estimate 0(/i -1 / 2 /7 -3 ) holds while the 

main part is given by (0.12) with T = T = Ch\ log h\. Now I want to prove the same under 
a larger upper bound for fi; further, increasing \x I want to recover a remainder estimate 
consistent with those in theorems 0.1, 0.3, 0.5 while keeping the main part given by (0.12) 
with T = T. To do this I need to estimate properly (1.22)-like correction term 

(4.3) r 1 f (F t ^ T (xAt) - XT»{t)) (TuQ}j) dr 

with T' = Tq derived in section 3 and T" = T 30 \ This will let me to rewrite (0.12) as a 
Weyl expression, may be with 0(/i 2 /?" 2 ) correction. Here Q = ip but in the analysis below 
I will need more general PDOs. 

This expression (4.3) is the sum of expressions of the type 

(4.4) h~ x J [F t ^ T ( X 2T{t)-XT{t))VuQl) dr = /T" 1 {F t ^ r XT{t)V uQ*)^ 

with x{t) = t _1 (x(|f) - and T = 2 k T", k running from 2 to [log 2 ( V / T")] . 

Proposition 4.1. Let assumptions of theorem 0.1 with 

(4.5) h~ 5 < fi < e(h\ \ogh\y 1 

be fulfilled and let Q be supported in {dist(x, T.) < 7}, 7 > /i" 1 / 2 . 

Then expression (4-3) with T' = e/i _1 7 _1 and T" = T does not exceed C\i~ x l 2 br z + 

C/2 2 h~ 2 31 ). 

Proof. Note first that expression (4.4) with T G [T, e/i _1 ] is negligible due to condition 
(0.2). So, until the end of the proof I can consider T" = e/i _1 . 

Further note that, as Q is supported in {|Zi| 2 > Cfih\ log h\} expression (4.4) with 
T G [T", T'\ is negligible as well. Then for a general PDO Q expression (4.4) does not 
exceed C^ 2 h~ 2 \ \ogh\. 

Then proposition is proven as fi < c(h\ log h| ) 2 / 5 and until the end of the proof I can 
assume that 

(4.6) {h\ log h\) 2/5 <fi<e{h\ log hi)' 1 . 

30 ) However in what follows I will consider some other parameters T > T" . 

31 ) Term C/i -1 / 2 /? -3 is probably of no need in some estimates here and below but it is in the remainder 
estimate anyway. 
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Further, as Q is supported in {dist(x, E) < 73 = c| log expression (4.3) would not 
exceed Cfi 2 h~ 2; y 3 \ \ogh\ x p?br 2 and as Q is supported in {(Z^ 2 < cfih} expression (4.3) 
would not exceed C/i 2 h~ 2 as well. 
So, the only zone in question is 

(4.7) {C Q ph<p 2 < Cph\ \ogh\} n {dist(x, E) x 7 , \Z X \ x p) 

with I log < 7 < e. In this zone one needs to consider expression (4.4). 

Using precanonical form and decomposition (2.16) one can rewrite (4.4) as the sum of 
the similar expressions but with u, Q, V replaced by u n , Q n and V respectively: 



(4.8) /T- 1 (F t ^ /j - lrX r(t)r / i/ n ^ y 

From the decomposed problem 

(4.9) (hD t - A n )u nn , = 0, 

(4.10) iw|t=o = S nn ,S{x' - y') 

one can prove not only that expression (4.8) is negligible as T > T p d = Chp~ 2 \ log h\ but also 
that its absolute value does not exceed Cph~ 2 'yp 2 (hp~ 2 7") s as t < T p . The summation of 

this expression with respect to T from e/i -1 to 00 results in C/i/?~ 2 7p 2 (/x/jp~ 2 ) s . 
Further, summation of the latter expression with respect to n results in 

CpLh- 2 1P 2 (php- 2 y 5 x p^fih)- 1 = Ch- 3 p 4 (php- 2 y 5 

since no more than p 2 (fih)~ 1 numbers n violate the ellipticity of A n . 

Finally, the summation with respect to p > c(fih) 1 ^ 2 results in Cp 2 h~ 2 'j. So, one can 
conclude that expression (4.3) with indicated T' and T" by absolute value does not exceed 
Cfi 2 h- 2 . □ 

One can improve the above error estimate under nondegeneracy condition. 

Proposition 4.2. Let assumptions of theorem 0.1 and condition (4-6) be fulfilled. Further, 
let Q be supported in {dist(x, E) < 7} with 7 > p?~ x l 2 . Finally, let condition (0.9) q with 
q = 1, 2 be fulfilled. 

Then expression (4-3) with T' = efi^ 1 ^ 1 , T" = T and Q supported in {dist(x, E) < 7} 

does not exceed 

(4.11) q Cp-^ 2 h- 3 + C{fih) 2+ ^ 2 \ log/>| K /r 4 7 . 
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Proof. There is no need to consider the inner core. Further, due to proposition 4.1 the case 
fi < Ch~ 2 / 5 is already covered. So, let us assume that p > /j~ 2 / 5 . 

Consider precanonical form (2.17) and define £ = |V*, p = \Zf\. Then even as p < 

p d = C(ph\ log/?!) 1 / 2 but £ > £ = f p one can still trade T' = ep~ 1/ ~f~ 1 to T" = e/x" 1 because 
then the shift for time T e [e/i _1 , T"\ is still observable. Here perturbation terms and their 
derivatives are at most C/i~ 2 <C £. 

On the other hand, the total contribution to the remainder of all partition elements 
with p < p, £ < £, dist(x, H) < 7 does not exceed Cp 2 £ 2+q ^h~ Ar which is exactly (4.11),. □ 

Remark 4.3. While the error estimate of proposition 4.1 is as good as I need in the general 
case, the error estimate of proposition 4.2 is as good as I need under condition (0.9), only 
under some restrictions to p and/or 7. 

4.2 Moderate magnetic field. II 

Now I want to trade T" = ep~ 1/ ~f~ 1 to a larger value To which provides the remainder 
estimate derived in section 3. Now I will need to use some canonical forms and thus 
different zones will be treated separately. 

4.2.1 I will start from the strictly outer zone. 

Proposition 4.4. Let assumptions of theorem 0.1 and condition (4-5) be fulfilled and let 
Q be supported in Z* ut fl {dist(x, T.) < 7} with 7 > CqP 5 " 1 ^ 2 . Then 

(i) In the general case expression (4-3) with T' = Tq and T" = T does not exceed 
(4.12), C/x- 1/2 /T 3 + C( Ai /)) (4+c ' )/2 /7- 4 7| log h\ K + C{ph)^' 2 h-^ 2 \ log h\ K 
with q = 0; 

(ii) Under condition (0.9), this expression does not exceed (4.12),. 

Proof. One needs to consider only Q supported in a zone in which in the proofs of propo- 
sitions 4.1, 4.2 term (4.3) with T' = ep^ 1 and T" = T was negligible. Then I am going to 
prove (4.12), without the second term. 

(i) Assume first that Q is supported in Z^ r 7 ) C Z* ut ,. 

Let us consider canonical form (2.54) and introduce function t by (3.38). If I > 1 (i.e. 
p > Cr~ 1 (ph)~ 1 ) I do not partition the "final" ball 6(0, 1) corresponding to the original 
(7; 7r _1 ) element; otherwise I make ^-subpartition. Then on any given £ subpartition 
element the propagation speed with respect to x 3 does not exceed /i~ 1 7~ 2 r^ p ; further the 
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propagation speed is of this magnitude as \d^ 3 H pn \ x l\ in this case the shift for time T is 
observable as it satisfies the logarithmic uncertainty principle 

yT 1 1 ~ 2 £ x T x £ > Cr 2 /i- 1 hj~ 3 \ \ogh\. 

Plugging T x // 7 one arrives to 

(4.13) £>£ l = C(/ih\\ogh\) 1/2 . 

Similarly, considering propagations with respect to other variables one arrives to the same 
conclusion as |Vl/l/ np | x £. 

Note that for given index p and for index n violating ellipticity I x £ p = (£* + r\p—p\fj,n) 
with some p < C (fih)~ 1 . 

Further, on each £* subpartition element condition (4.13) is fulfilled for all indices p 
as £* > l\ and condition (4.13) is violated for no more than Co£i(rfih)~ 1 (^> 1) indices 
p as f < £]_. Furthermore for each such index p ellipticity is violated by no more than 
Co(^i7(r/i/7) _1 + l) indices n. 

Therefore the total contribution to (4.3) of all singular subelements (i.e. subelements 
with £* < Cq£\) belonging to 2( r , 7 ) does not exceed 

(4.14) C/i 2 /)- 2 7 2 ^ 2 x (liirfihy 1 + l) x ^(r/i/?)" 1 + l) x £\; 

since the second and the third factors do not exceed C£i(rfih)~ 1 and C\ log h\ K respectively, 
I get C/)~ 4 (/i/7)' 3+C7 ^ 2 7 2 | \ogh\ K . This expression sums with respect to r, 7 to the same 
expression calculated at the maximal values of 7 and r = 1, and multiplied by | \ogh\, 
which is exactly expression (4.12) q . I increase K as needed. 

(ii) Zone Z* ut u is treated in the same way. □ 

While estimates (4.11)^ with q > and even (4.12) q with q > 1 are sufficient for my 
needs, (4.12) is not and needs to be improved to (4.12) q with q arbitrarily close to 1: 

Proposition 4.5. Let assumptions of theorem 0.1 and condition (4-5) be fulfilled and let 
Q be supported in Z* ut fl {dist(x, E) < 7} with 7 > Co//" 1 / 2 . Then expression (4-3) with 
T' = Tq and T" = T does not exceed 

(4.15) C/i- 1 / 2 /)" 3 + C/i 2 /j" 2 7 | log h\ + C(/i/j)( 3+£7 >/ 2 /r 4 7 2 | log h\ K 
with q < 1 arbitrarily close to 1. 
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Proof. Again, one needs to consider only Q supported in a zone {\Z\\ > C(fih\ log h]) 1 ^ 2 } in 
which in the proof of proposition 4.1 term (4.3) with T' = e/x -1 and T" = T was negligible. 
Then I am going to prove (4.15) without the second term. 

(i) Assume first that Q is supported in 2i r ,j) C Z* utl . One needs to consider only 
P < Co(rfjLh)- 1 . 

Let us introduce the scaling function and ^-admissible partition the in the "final" ball 
6(0, 1) corresponding to the original (7;7/" _1 ) element: 



(4.16) t k = 4 = min e^ 1 ]T | VU,p| (m+1)/(fe+W) ) + i k 



l/(m+l) 

j<k 

with k = m, ( = 7r _1 and 

iJcvi/(*+i) 



(4.17) £ k =(fih\\ogh\ K y 

As m = 1 this function coincides with £* given by (3.38). 

Consider first partition m-singular groups ((p, n), element) i.e. groups such that 

(4.18) \d a A np \ < cC^" |Q| Va : \a\ < k 

with k = m, £ = £ m and ( = Cm- 
Then contribution of all m-singular groups in 2( r , 7 ) could be estimated in the same 
manner as (4.14): 

(4.19) C> 2 /r 2 r 2 7 2 x (C£ k {"flih)- X + l) x (c^ifih)- 1 + 1 

with the same k, £ and ( as in (4.18). Here the second factor estimate the number of indices 
p involved and the third factor indicates the number of indices n violating ellipticity of A pn 
for given p. 

The groups which are not /(-singular are k-regular, i.e. they satisfy (4.18) with some 
I > 2£k and also 

(4.20) \d a A np \ x Ck£ k ~ H for some a : \a\ < k. 



Then one can rescale such element to 6(0, 1). After this rescaling conditions (4. 18), (4. 20) 

are fulfilled with £ re =l and ( re = f (£ k . So, one can apply the same rescaling and partition 
to it with (k — 1) instead of k and with crucial £ (separating (k — l)-singular from (k — 1)- 
regular) equal to £k-i/ ' £k- This means that if one returns to the original element, the radius 
of the /(-singular element would be £ k every time. 
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This continuation however has a new property: every time when regular or singular £'- 
subelements appear inside of some regular elements, their relative densities do not exceed 
£'/£ (all the time I refer to the same scale) and therefore the absolute density of all £'- 
subelements does not exceed Cfi' 1 . 

Therefore the total contribution of all /(-singular elements does not exceed (4.19) mul- 
tiplied by ik^m 1 

(4.21) «" m V 2 /rV7 2 x (c^ 1 ^)" 1 + e) x (V+V^r 1 + l) 

with £ = tk and ( = ^r~ l since I know that ( does not increase (even if it is different in 
the different elements). Then two last factors do not exceed | log h\ K for sure and the total 
contribution of all singular ((p, n),subelement) groups does not exceed £~ 1 /x 2 /7~ 2 7 2 r 2 | \ogh\ K 
and summation with respect to (r,-f) results in the same expression with the maximal 
values of 7 and r = 1. 

Microlocally this is sound as logarithmic uncertainty principle 

(4.22) P k > /i" 1 /?^" 3 ! \ogh\ K 

and to satisfy it for any 7 ^> p. -1 / 2 /" 1 / 2 one needs to take k = 1. 

This leaves me with 1-regular groups and there are three types of them: 

(a) With Q£ 2 < / -fr~ 1 fih\ \ogh\ K , which are covered by the same estimate (4.21); 

(b) With (I 2 > 7r" 1 / u/?| log h\ K ; for them the shift fi~ 1 ^y~ 2 r(£T with T = yU _1 7 _1 is 
observable since /x -1 r 2/ ~f~ 3 (£ x /i _1 7 _1 x £ > Cr 2 fi~ 1 h / -f~ 3 ; 

(c) 0-regular groups but then A pn is just elliptic on the corresponding elements. 

(ii) Zone Z* ut u is treated in the same way. □ 

4.2.2 In the general case my purpose is the remainder estimate 0(fi~ 1 ^ 2 h~ 3 + fi 2 h~ 2 ). 

Proposition 4.6. Let assumptions of theorem 0.1 be fulfilled. Then 
(1) As 

(4.23) h- 5 < n < Jio = = h 5 - 2 ' 5 

J e(x, x, r)ip{x) dx is given modulo 0([i~ 1 / 2 h~ 3 ) by expression (0.12) with T = T; 
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(ii) Moreover, as 

(4.24) h 5 ~ 2/5 < y, < h^ 1 

and tp(x) is a nice 32 ^ function supported in {\xi\ < 72 = t L ~ 5 '} J e(x, x, r)^(x) c/x given 
modulo 0(/i~ 1/2 /r 3 + /i 2 /r 3 ) fry expression (0.12) with T = T 

Proof. Combining propositions 3.2, 3.4 and 4.5 I conclude that this estimate holds in frames 
of (i) for if; supported in Z* ut and in frames of (ii) for ip supported in Z* ut D {|xi| < /i _<5 }. 

In propositions 3.11, 3.14, 3.16 and 3.19 it was essentially proven that this estimate 
holds for if) supported in the near outer zone Z' out , inner core, inner bulk zone Z mn ji and 
true inner zone Z lim j respectively. □ 



4.3 Intermediate magnetic field. I 

4.3.1 Now I want just to calculate (0.12) with T = Ch\ \ogh\. It is well-known (see for 
example [Ivr5]) that 

Proposition 4.7. Let n < h~ 1+s with arbitrarily small 5 > 0. Then under condition (0.2) 
(i) Asymptotics holds 



(4.25) h- 1 f (F t ^ h -i T XT(t)r x uQ;)dT~ «".m,o(x,A)/7- 4 + 2m+2 V 



-4+2m+2n, 2n 

as \\\ < e, T < T < e/i" 1 . 

(ii) Moreover, as Q = Q(x) K n ,o,ct(x, ■) = ^n,o,i(x, )Q(x) and 



/oo 
%((A - t) Th- 1 )S MW (x, r)dr~J2 ^,o,/(*, A)/T 4+2m+2 V 



-4+2m+2n,.2n 

with the Standard Weyl Expression 



(4-27) k ,o,/(x, A)/j = £ w (x, r) = -^(2r + V) 2 ^. 

Combining propositions 4.6(i) and 4.7(i) I arrive to 



Corollary 4.8. Under conditions (0.2) and (4-23) modified asymptotics (0.3) (with £ 
replaced by £ w ) holds. 



MW 



••!2 ) 



F.e. ip(x) = V'i(^i/72)^2(x2,x 3 ,x 4 ). 
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To prove theorem 0.1 under condition (4.23) it is sufficient to estimate properly 

(4.28) | / (£ MW (x, A) - £ w (x, A))^(x) dr\ 
which due to proposition 4.7(h) is equivalent to 

(4.29) | y^ MW (x,A)- Th~ l j°° ^{(\-T)Th- l )E M ' N (x,T)dT^4,(x)dT\. 

Proposition 4.9. (i) Under conditions (0.2) and (4-23) both expressions (4-28), (4- 
with A = do not exceed C /i _1//2 h~ 3 . 

(ii) Under conditions (0.2) and (4-24) both expressions (4-28), (4-29) with A = and a 
nice^^ junction ip supported in {|xi| < 72 = H~ S '} do not exceed Cii~ l ^ 2 h^ + Cfi 2 h~ 2 . 

Proof. Proof is standard, based on the same scaling functions and partition as in the proof 
of proposition 4.5 as ip is supported in Z* ut . These arguments work in the other zones 
as well since I do not need uncertainty principle anymore. I leave the easy details to the 
reader. □ 

Now propositions 4.6(i), 4.7 and 4.9(i) imply 

Corollary 4.10. Under conditions (0.2) and (4-^3) asymptotics (0.3) holds. 

4.3.2 In view of theorem 0.1 nondegeneracy condition should be used only as 

(4.30) ch 5 ~ 2/5 <fi< Ch- 1 

with 5 = 0. Then, under condition (0.9)i the remainder estimate derived in section 3 is 
0(/i~ 1/2 /)~ 3 + /i 3 / 2 /?" 3 / 2 " 5 ) i.e. it is 0(/i~ 1/2 /)~ 3 ) as fi < /?^ 3/4 and then there is no need to 
invoke (0.9) q with q > 2. However it makes life much easier and I will do it right now leaving 
more difficult analysis for the later; so 5 > is arbitrarily small in this subsubsection. 
First of all one needs to prove 

Proposition 4.11. Let conditions (0.2) and (0.9)^ be fulfilled and jj, < h 6 ^ 1 . Let iJj be a 

nice function supported in {|xi| < 7} with /i 1 / 2 " 5 < 7 < e. Then J e(x, x, r)^(x) dx is given 
by expression (0.12) with T = T with an error not exceeding (4.12) q . 

Proof. Combining propositions 3.7 and 4.4 I conclude that the statement holds for ip sup- 
ported in Z* ut D{|xi| < 7}. 

In propositions 3.14, 3.12, 3.17 and 3.22 it was essentially proven that this estimate 
holds for ip supported in the near outer zone Z' out , inner core, inner bulk zone Z- mn ii and 
true inner zone Z mn i respectively. □ 
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Proposition 4.9 and 4.7 imply 

Corollary 4.12. Let conditions (0.2) and (0.9) q be fulfilled and fi < h 5 ^ 1 . Then 

(i) One can define the main part of asymptotics by (0.12) with T = Ch\ log h\ (or equivalently 



under condition (0.9) q with q > 1; in particular fj,2 = H 1 ^ 2 \\ogh\ K and ^ = fi 4 / 7 |log/?| K 

I refer to the case of (4.31) q as a moderate magnetic field. While I can get rid off 
logarithmic factors, there is no point to do it right now. 

To prove theorem 0.3 one needs to estimate properly expressions (4.28) and (4.29). 



Proof. Proof is standard, based on the same scaling functions and partition as in the proof 
of proposition 4.4 as ip is supported in Z* ut . These arguments work in the other zones 
as well since I do not need uncertainty principle anymore. I leave the easy details to the 



Corollary 4.14. Theorem 0.3 is proven as fj, < fi q . 
4.4 Intermediate magnetic field. II 

4.4.1 Now I want to prove theorems 0.1, 0.3 as fi > Jx q defined by (4.23), (4.31)^ for q = 
and q > 1 respectively. To do this I need to consider expression (0.12) localized in this 
singular zone because in the regular zone one can always apply (0.12) with Q = I even if 
the expression is not as explicit as in the case Q = I . 

Proposition 4.15. One can rewrite with the same error as a remainder estimate derived 
in section 3 expression (0.12) as 




Proposition 4.13. Let conditions (0.2) and (0.9) q be fulfilled and \i < h s 1 . Let ip be a 

nice function supported in {|xi| < 7} with [i & ~ 1 ! 2 < 7 < e. Then 

(i) Both expressions (4-28) and (4-29) do not exceed (4.12) q . 

(it) In particular as 7 = 1 both expressions (4-28) and (4-29) do not exceed C/x -1 / 2 /?" 3 as 

/i < /J q . 



reader. 



□ 



(4.32) 
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with 



(4.33) ^ w (x,0) d ^const^(^(^ pn )Op n )| (x/; e;)=xu x Mx^x^/T 2 



".p 

MW 



exactly the same constant as in S 

Proof. Proof is rather standard: I consider canonical form of the operator and in this form 
I apply the method of successive approximations with operator in question A pn (x, / u _1 /7D x ) 
and unperturbed operator A p , n (y, fi~ 1 hD x ). I leave details repeating those in my multiple 
papers to the reader. □ 

Remark 4.16. Obviously, defining by (4.33) but with A pn replaced by A° pn one would 
get exactly 

(4.34) £$ w (x,0) = 

const ^0(v(x) - (2n + l)fihf 2 (x) - (2p + l)^i(x)) h{x)f 2 {x)Q pn ^ 2 h~ 2 

where A° pn is obtained from A pn by replacing perturbation terms B pn by 0. 

While one can take Q = I and then Q np will be the diagonal elements of ip transformed 
according to (2.55) the expression (4.33) is not very explicit either because of presence of 
B pn and the similar terms in (2.55). Because of this I want to take Q supported in as small 
zone as possible; the only restriction so far is that Q = I in the singular zone 

(4.35) Q sing = {£* < h = C(rfih\ log h\f' 2 , \p-p\< r^h)' 1 , 

\n-n\< C 7 r" 1 £ 2 ( / u/))~ 1 + l|. 

Therefore I will take Q supported in this zone (with increased C). 

Furthermore I actually need to consider not (4.32) but only its correction with respect 
to what is given by magnetic Weyl formula; namely 

(4.36) |(^ w (x,0)-^ w (x,0))^(x)c/x. 
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4.4.2 Now my goal is to estimate expression (4.36). I do it first under condition (0.8) q 
(including formally q = 0) and in the next subsubsection I improve case q = in the same 
manner as proposition 4.5 improves proposition 4.4. 

Proposition 4.17. Let norms of perturbation operators do not exceed ce, 

(4.37) e < fih 

and Q be supported in Z* ut . Then expression (4-36) does not exceed 

( C£ 7 2 (^) (1+q)/2 ^ 4 |log/)| K q>2, 

(4.38) , I C 1 2 (e{fih)h-' i + Ce 1/2 {fih) 2 h-^\\ogh\ K q = 1, 

{ C^C/^rVl log h\ K + C 7 2 (/i/>) 2 /r 4 c? = 0; 

Proof. Let us introduce a scaling function by (3.38). Further I replace £ by min(£, r) + £ 1 / 2 . 
Then in as £ = £ pn > C£i, the contribution of all £ pn groups to (4.36) does not exceed 
the left-hand expression of (4.14) with {£ instead of l\ = C(r/x/)| log /?|) 1//2 ), multiplied by 
Ce£~ 2 : 

(4.39) C^h' 2 j 2 r^+ x (^{rfih)- 1 + l) x (i^r/ih)- 1 + l) x ^ 2 e; 

one can prove it easily by considering such elements and integrating by parts if I > ce 1 / 2 . 
Here the third factor does not exceed | log h\ K for sure as I 2 < fih\ log h\ K and 7 < r 2 . 

Then as q > 2 (4.39) sums with respect to £ to its value as £ reaches its maximum 
min(£i, r); then it sums with respect to r s to the the first line in (4.38)^. 

As q — 1 one gets 

C{fih)h~ 4 n 2 (e + (/i/jjet 1 ) I log h\ K ; 

one can estimate £~ 1 e by e 1 ! 2 as I > e 1 ^ 2 and the summation results in the second line in 
(4.38),. 

The same approach works as q = as well but results in an unwanted now factor 
I \ogh\ K at C( / u/)) 2 /j~ 4 7 2 . However, as £ > e 1 ^ 2 ] \ogh\ Kl with the large enough exponent K% 
one can estimate the fourth factor in (4.39) by | \ogh\~ Kl and all summations will result 
in C(fih) 2 h~ 4 ~f 2 \ \ogh\ K °~ Kl with Kq independent on K\\ choosing K\ large enough one gets 
C y u 2 /)- 2 7 2 . 

On the other hand, one can estimate contribution of all groups with £ < e 1 ^ 2 ] log h\ Kl does 
not exceed the product of the first three factors in (4.39) with £ = e 1 ^ 2 ] log h\± (since one 
does not need to use an integration by parts then) which does not exceed after summation 
with respect to r, 7 the third line in (4.38),. □ 
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I remind that the "perturbation" terms in (2.50) are 
£ ^x((2p+l) A ifry)V-^ft l + 

2q+2/+2y>3 

E ^ >< ((2«+ l)^)*((2p + l)A*fry) 'f-n-^-'tfr-^-^v-Vh' 

2k+2q+2m+2s+2l>3 

and due to (2.51) their norms do not exceed C/i~ 2 r7~ 3 and thus one can take 
(4.40) £ = C min(/i" 2 7 - 3 , /i^V 1 ) 

as \xi\ x 7 since r < min(l, fj}~ s j 2 ) in 2* ut . Then condition (4.37) translates into 

def / i i/o * o u f /i" 1 /)^ 1 / 3 as u < h 5 '~ 2 / 3 , 

(4.4!) ^^Hfc.in^-Va,^)^^ ^ <* - ^ 

4.4.3 Consider g = first. Then plugging (4.40) into (4.38)o and taking the sum with 
respect to 7 one gets the same expression calculated as 7 = 1 i.e. Ch~ 3 \ \ogh\ K + Cfi 2 h~ 2 
which is 0(fi 2 h~ 2 ) as /1 > /) _1 / 2 | log/)| K . Therefore due to due to propositions 4.6, 4.11, 
4.15 and 4.17 I arrive to 

Corollary 4.18. Estimate {0.3) holds as fi > /r 1/2 | \ogh\ K . 

Remark 4.19. As h 5 ~ x < /i < cbr 1 the contribution of zone {|xi| < h s } to the asymptotics 
is 0(/?~ 4+<5 ) = 0(fi 2 h~ 2 ) and one does not need to use 4.6 and 4.11. 

So estimate (0.3) remains to be proven as h 5 ^ 2 ^ 5 < fi < /? _1 / 2 | log h\ K ; I will do it in the 
next subsubsection. 

Consider case q > 1 now. Plugging (4.41) into (the first term of) (4.38)^ one gets 
C(fihy i+q ^ 2 h~ A min(yU _2 7 _1 , yU 2_5 7 ); after summation with respect to 7 > 72 one the same 
expression (4.38) q calculated at 

, s def , x i/„ f _1 /T 1/3 as/i< /? 5 '- 2/3 , 

(4.42, , = ,3*^,^)^^ 

which is C(/i/7)(^ 1 )/ 2 /?- 8 / 3 )|log/7| /< = 0(/i~ 1 / 2 ) as /i < //'~ 2 / 3 and C^)^ 2 / 2 /?- 5 / 2 /!" 5 as 
fi > h s ~ 2 / 3 ; the latter expression is 0(/i -1//2 /)~ 3 ) for q >2. 
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So, as q = 2 one gets 0(fi 1//2 /7 3 ) as an estimate of the contribution Z* ut D {|xi| > 72} 
into 

(4.43) I J (e(x, x, 0) - £ MW (x, 0))^(x) dx\ 

On the other hand, due to propositions 4.9, 4.11 and 4.13 and their corollaries and also 
propositions 3.20, 3.21, 3.22 contribution of zone Z* ut D {|xi| < 72} U Z' out U Z- mn to (4.42) is 
estimated by (4.12) q calculated at 7 = 73. One can see easily that the result is 0(ji~ l l 2 h~ 3 ) 
as /i < h 5 ~ 2 / 3 , q > 2; as /i > h 5 ~ 2 / 3 one gets 

c^ 2 h- 3 + c(^) (2+q)/2 ^y +1/2 + c(/i/)) (1+ ' )/2 r 3 / 

where the third term is far less than the second one. 

Then for q = 2 this result does not exceed C/i _1 / 2 /7 -3 as [i < h s ~ 2 / 3 and statement (i) 
of the proposition below is proven. On the other hand, as ji > h s ~ 4 / 5 72 is below of the 
bottom of Z* ut (see 2) and therefore statement (ii) is proven: 

Proposition 4.20. Let conditions (0.2) and (0.8)2 be fulfilled. Then 

(i) As /i < /j 5 ' -2 / 3 expression (4-42) does not exceed C^i^^h' 3 and therefore estimate (0.6) 
holds with ^ = 0; 

(ii) As /? 5 '~ 2 / 3 < jjL < ch^ 1 contribution of Z* ut into the (4-43) does not exceed C/i _1 / 2 /?~ 3 
while contribution of Z' out U Z- mn (including the inner core) does not exceed 

(AAA) C^ 1/2 h- 3 +Cfi 5/2+8 h~\ 

Remark 4.21. Plugging (4.40) into the second term of (4.38)i one gets after summation its 

value as 7 = 1 i.e. Cjih~ 2 \ log h\ K . Because this and other ugly terms and because estimate 

under condition (0.8)i is not a part of my core theorems, I am no more considering q = 1, 

leaving to the reader either to derive some estimate of (4.43) or to derive 0(fi~ 1 ^ 2 h~ 3 + 
^3/2^-5-3/2) estimate f 

some more complicated expression. 
4.4.4 To finish the proof of theorem 0.1 I need to prove 

Proposition 4.22. Let condition (0.2) be fulfilled, h 5 ~ 2 ^ 5 < /i < Ch" 1 and ip be a nice 
function supported in {72 = < |xi| < e} with small enough exponents 5 > 0, 5' > 0. 
Then 



(4.45) I ( (S(x, 0) - £ MW (x, 0))^M dx\ < C^h' 2 . 
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Proof. One needs to consider only zone r > [i~ 5 ' . Here and in the assumption 5' > could 
be taken arbitrarily small. 

Let me introduce scaling function £ m as in the proof of proposition 4.5 and the corre- 
sponding partition. Then due to the same analysis as in proposition 4.5 the contribution 
of all the m-regular groups with to the remainder is less than C/i _1 / 2 /7 -3 + Cfj>~ 2 ~ 5 h~ 2 . 

So one needs to consider only m-singular ((p, /7),subelement) groups. Corresponding 
subelements are of the size £ = (fj,h\ log h\ K Y^ m+1 > and on each such group |<9 2 *4. pf) | < 
C(/i/?| log /)|)( m_1 )/( m + 1 ). Consider e/i/)-subpartition. Then for each element of it and for each 
index p 

\V{x)-f 1 {x){2p + l)ph-f 2 {x)(2n + l)ph\ > e ph\n-n(p)\ 

and the same is true for this expression perturbed by B pn = 0(e), e = yU _2+4<5 ', and also for 
each element there exists index p 

\V(V{x)-f 1 {x){2p+l) f ih-f 2 {x){2n + l)ph)\ > e ^h\p - p| 

as n = n(p). 

Then repeating arguments of the proof of proposition 4.20 one can see easily that the left 
hand expression of (4.45) does not exceed C/i~ 2 /?~ 2 (l + e Xl/ ( > 1 (/u/)/c) _1 ) < Cfi~ 2 h~ 2 since 
e iih so perturbation does not violate ell ipticity of elliptic A pn . □ 

Corollary 4.23. Theorem 0.1 is proven completely. 
4.5 Strong magnetic field 

To prove theorem 0.3 or slightly worse estimate under condition (0.8)2 one needs to improve 
estimate (4.44) of the contribution of Z' out U Z mri (including the inner core) into remainder 
as /i > /? 5 ~ 2 / 3 and to calculate correction term (3.127)-(3.128) as fi > /7~ 6 / 7 | log/j|~ 2 / 7 . 

4.5.1 To improve estimate (4.44) of the contribution of Z' out U Z- mn (including the in- 
ner core) into remainder without calculating the correction term one can notice that this 
expression is given by (0.12) with T = T Q = Ch\ log h\(p 2 + £ 2 ) -1 matching one for 2- 
dimensional magnetic Schrodinger operator (2.26) and therefore the asymptotics with the 
Weyl expression for (2.26) which one can rewrite easily as 

(4.46) {2n)- 2 ph- 3 9(2T+ V - {2n + l)[ihf 2 - (2p + l)p,hf^\ hh^fg 

where I skipped perturbation B' n = 0{n~ 2 ) which under condition (O.82) leads to the relative 
error 0(/i~ 2 ). 
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While each expression (4.46) is of magnitude fih~ 3 and after integration over (Z' out UZ- inn ) 
it acquires factor 7 = / u~ 1 / 2 r 1 / 2 and after summation over n it acquires factor (fih)" 1 ; so I 
get a required expression 

J £ MW (x, r)^(x) dx 

of magnitude /Z -1 / 2 /) -4 as ^ is supported in (^ ut U Z mn ), an error has magnitude 

^-l/2^-4 f 5/ 2 x ^-2 = ^-5/2^-4 = 

Therefore I arrive to (almost) final result: 

Theorem 4.24. Let F be of Martinet-Roussarie type and condition (0.2)be fulfilled. Let ijj 
be supported in the small vicinity ofH. Then 

(i) Under condition (0.8)3 estimate (0.6) holds; 

(ii) Under condition (0.8)2 the left hand expression of estimate (0.6) does not exceed Cfx~ 1 ^ 2 h~ 3 + 
Ch s ~ 5 / 2 with an arbitrarily small exponent 5 > 0/ 

So far ££w w defined by (3.127) - (3.128) and is 0(^ 4 h- 3 / 2 \ log/)] 1 / 2 ); a more explicit 
representation is given by (4-51). 

4.5.2 First of all note that formula (3.57) [Ivr6] one can rewrite in the form which does 
not require g^lr = Sj^. 

Lemma 4.25. One can rewrite formula (3.57) [Ivr6] as 

/ <;„l/( !/ + 1 )/( 2l/ U- 1 A / \ 

(4.47) = pny^h-'h^K-^V^^^G ^ M )VS'\ T 

where h~ = fj}' v h, g' = g^i = g u g so ^fg' 6x2 is a Riemannian density on E, 

(4.48) <P= (f-dist(x,r) 1 - ly )| L , 

dist(., .) is calculated according to the Riemannian metrics (gjk) and function G is given by 
(3.53) [Ivr6]. 

Really, in this formula V, <fi and ^fg 1 dx2 are invariant with respect to change of the 
coordinates while multiplication of operator by uj 2 is equivalent to substitution g^ k 1— > u 2 g k 
and V 1— > lu 2 V which leads to f h cu 2 f, g' \— > uj~ 2 g', dist(x, T.) \— > cu" 2 d\st(x,T.) (up 
to a factor equal 1 at E) and <j> ^ LU u+1 (f); so yH" 1 )/^)^ 1 /" an d V^-^cj) 1 ^^ are 
invariants. 
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As v = 2 and fj, < Ch~ x (3.47)-(3.48) become 

(4.49) = (27r)- 3 / 2 /7- 1 ^ 1 /2 ft .-i/2 V i/8 i/4 G 

with an error 0(h~ 1 h) = 0(/i -1 / 2 /? -1 ) and 

(4.50) </> = \\Vf\\ = {Y,g Jk djfd k f) 



2-kH 



V2| 



where norm |.| corresponds to the Riemannian metrics (gjk)- 

Similarly, repeating arguments of [Ivr6] leading to calculation of the correction term 
(minimal modifications I leave to the reader) one can calculate (3.127) — (3.128) modulo 
0(fi~ 1 ^ 2 h~ 3 ) deriving the following modification of (4.49)-(4.50): 



(4.5i) crv) 



(2,)- e (y-(2 B+ i)6^)"y/«c( *( v ~ (2 " + 2 ^ ft)3/v ' /2 )6v? 



where = g^V and y/g 7 dx' is a Riemannian density on Z, is now defined by (4.50) with 
derivatives taken only along K 1; f — f\ 



(4.52) 



I Vk^iIki I 



where |.| corresponds to the Riemannian metrics (gjk) restricted to K4. One can see easily 
that only zone {dist(x, A) x 1} contributes after integration since otherwise 
|Vz((V - (2n+ l)f 2 fih)(f)- 2 / 3 ) \ is disjoint from 0. 

Theorem 4.26. Statement of theorem 4.24 holds with £j££ = 0( / u 5/4 /?~ 3/2 ) defined by 
(4.51) -(4.52). 



A Additional Results 
A.l Proof of Theorem 0.5 

Proof of Theorem 0.5 cannot be obtained by the simple rescaling since fe/i and fi/i/|xi| scale 
differently. However, arguments of sections 1, 3,4 work with little or no modifications: 
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A. 1.1 In section 1 propositions 1.5-1.9 do not require condition (0.2) while proposition 
1.10 holds under condition (0.10) instead of (0.2) where T = Ch\ \ogh\/£ with 

(A.l) £ = emax(|xi|, \V\) +£, £ = C/ih\ log h\, 

instead of T = Ch\ log h\. 

Further, propositions 1.11, 1.12(i) do not require assumption (0.2) and proposition 
1.12(h) holds with modified T under condition (0.10). Therefore proposition 1.13 also 
remains true: one needs to estimate \F t ^ h -i T XT(t)f~(uip y Qy)\ as T = T and this is done 
easily as £ > 3£ by rescaling and standard Schrodinger operator analysis; as £ > 3£ it can 
be done easily by more crude approach as well. 

Furthermore, propositions 1.14, 1.15 do not require assumption (0.2) and due to the 
previous modifications arguments of proposition 1.16 lead to the same estimate 0(/i~ 1 / 2 /? -3 ) 
for the contribution of Z out to the remainder estimate under condition (0.10). 

Finally, propositions 1.17-1.18 also do not require assumption (0.2) and proposition 1.19 
holds under condition (0.10) instead of (0.2); proposition 1.20 also does not require (0.2) 
and propositions 1.21-1.22 hold under condition (0.10). 

A. 1.2 Section 2 does not require condition (0.2) at all. 

A. 1.3 In section 3 analysis in the strictly outer zone Z* ut (subsection 3.1) leading to 
propositions 3.2, 3.4 does not require condition (0.2); the analysis and condition (0.10) 
leads to proposition 3.6 which is the worthy substitution for proposition 3.7. 

In the arguments subsections 3.2 and 3.3 in the absence of condition (0.2) one again 
should pick up 7= Ch\ log h\/£ with £ defined by (A.l); in comparison with the arguments 
of these subsections the factor £ _1 appears but it is compensated by a factor (| W\ + fih) 
which appears due the change of range of factors p 2 or (2/7 + l)/i/j. This would add an extra 
term Cp}^ 2 h~ 2 \ \ogh\ to the remainder estimate which is subordinate as p < C(h\ log /? | ) 1 . 

On the other hand, as p > (h\ log/?!)" 1 I can take e£-admissible partition with respect 
to x 1 and note that as 1/1/ < e 1 ph at some partition element, it will be classically forbidden 
and therefore its contribution to the asymptotics would be 0; moreover, with the rescaling 
method I can in the right-hand expression of estimates and using condition (0.10) I can 
take T = Ch there thus getting rid off the logarithmic factor: 

Proposition A.l. The total contribution of zones Z' QUt and Z- mn \ Z per to the remainder 
does not exceed C[x~ x l 2 h~ z under condition (0.10) while the main part of asymptotics is 
given by (0.12) with T = Ch\ \ogh\. 
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Moreover in the vicinity of {1/1/ = 0} condition (0.10) also reads as 
(A.2) |V Z (V<T 3/2 )| >eo 

which kills periodicity for all n not just for all n but one as it would be the case far from 
{ 1/1/ = 0} thus leading to 

Proposition A.2. The total contribution of zone Z per to the remainder does not exceed 
C^ -1 / 2 /?" 3 under condition (0.10) while the main part of asymptotics is given by (0.12) with 
T = T modified as above and no correction term is needed. 

A. 1.4 In section 4 condition (0.10) allows me to trade easily T = Ch\ \ogh\ to T = efi' 1 
to T = e/i _1 7 _1 to T = efi'-f 2 as each next expression is larger than the previous one 
thus leading to theorem 0.5 as fi < Z?" 5-1 covering everything by a moderate magnetic field 
approach. 

Alternatively as fi > the strong magnetic field approach works well: in Z* ut one can 
just skip 0(/i~ 2 7~ 3 ) perturbation thus replacing by £ MW as in proposition 4.17 with 
an error well below /i _1//2 /7~ 3 ; similarly in Z' out U Z- mn one can just skip 0(/z~ 2 ) perturbation 
thus getting magnetic Weyl expression again. 

I leave details to the reader. 
A.2 Special case 

Consider the special case of operator A defined as A\ + An where A/, An are operators of 
type (0.1) in variables x/ = (xi,X2) and x// = (x3,X4) respectively with magnetic intensities 
A(xi) = |xi| and f 2 (xu) = 1: 

(A.3) A, = /j 2 D 2 + (hD 2 - /ix 2 /2) 2 - 1, 

(A.4) A,, = h 2 D\ + (hD 4 A - fix 3 ) 2 . 

One can separate variables and prove spectral asymptotics with the remainder estimate 
h~ 2 Ri where /?i is the remainder estimate for operator A\ + r for the worst possible r > 0. 
According to [Ivr6] R\ = C^/i" 1 / 2 /) -1 ). 

Then the main part of asymptotics is exactly as in estimate (0.6) with the correction 
term is given by (4.51) and this term is 0(h~ 3 h 1 ^ 2 ) 33 \ 

33 ' So one loses factor /ih in comparison with theorem 4.26 in the estimate of the correction term. 
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The leading part of the correction term is rather irregular: it is (2n) 2 h 3 h 1 / 2 K 2 l n with 

(A.5) /(*.») S £ (y - (2n + l)fe)' / V^ c( 5 °( v/ - P" ^^^V ^l/ 

with ft = /i 1/2 /7 < £ = /i/7 < 1. 

Obviously that /„ = 0(1) and but it is not clear that no better estimate like l n = 0(e a ) 
is possible; one can see easily that sup e / xe ^ xS \l(s', h')\ > c~ x e. 

I hope that some readers will be able make a numerical experiments. 

A. 3 About term C[i 2 h~ 2 in the remainder estimate 

Is term 0(fi~ 2 h~ 2 ) really needed in the general case? The answer most likely is positive. 
Consider operator A defined as A\ + An where 

(A.6) A, = h 2 D 2 + (hD 2 - fix 2 /2) 2 -l-kxt 

and An is given by (A. 4). Then 

(A.7) W pn = V-{2n+ l)f 2 /ih - (2p + l)f lfi h = 1 + kx x - (2n + l)fih - (2p + l^/ih 

and as k = (2p + l)fxh, 1 = (2n + l)p/? then {W pn = 0} = {x x > 0}. 

While it does not mean that {A pn = 0} = {xi > 0} because of perturbation, I believe 
that perturbing slightly A\ one get achieve the latter result. I leave it to the reader. 
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